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A Method for Predicting Interchannel 


Modulation due to Multipath Propagation 
in FM and PM Tropospheric 
Radio Systems* 


By C. D. BEACH and J. M. TRECKER 
(Manuscript received August 7, 1962) 


This paper describes a method for predicting the magnitude of interchan- 
nel modulation due to multipath propagation on angle-modulated tropo- 
spheric scatter radio systems. Values of signal-to-intermodulation ratio, 
S/I, are calculated for various pairs of signal reflecttons in the troposphere, 
taking into account the base bandwidth and frequency deviation of the 
system, the antenna patterns, the path geometry, and climatic conditions 
during the “worst month” of propagation. The lowest value of S/I (worst 
intermodulation) computed for such pairs of signal reflections ts then cor- 
rected envpirically to account for multiple reflections. The result represents 
the median value of S/I expected during the worst month of transmission 
on a specified path. 

The method yields results that, when compared to measured results from 
four widely different paths normalized to worst-month conditions, have a 
standard error of estimate of about 2.6 db. 
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I. INTRODUCTION 


Various theories have been proposed to explain the mechanism by 
which energy radiated from a transmitting antenna is reflected or re- 
fracted in the atmosphere and ultimately arrives at a receiving antenna 
beyond the horizon.!?:? Apart from the merits of any particular theory 
of propagation, the geometry of a typical over-the-horizon, or “tro- 
pospheric scatter,’ path suggests an inherent multipath transmission 
problem. As illustrated by Fig. 1, components of signal reflected to the 
receiving antenna from successively higher points in the atmosphere 
are delayed by increasing amounts behind the earliest arriving com- 
ponent. 

When such time-delayed RI signals — called echoes henceforth — 
are demodulated by an FM or PM receiver, they are converted into 
energy at undesired regions of the baseband. In a frequency-division 
multiplexed, multi-channel voice system, this undesired energy is com- 
monly called “interchannel modulation,” or simply “intermodulation.” 
It adds to other distortion energy generated by equipment nonlinearities 
and to the ever-present white noise, to produce the total background 
of noise faced by the desired speech signal. Whether or not ‘“‘path inter- 
modulation”? noise is objectionable in a given system depends on its 
magnitude relative to these other components of ‘noise. 

The general theory of interchannel modulation due to transmission 
echoes in angle-modulated systems has been treated by Bennett, Curtis, 
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Fig. 1 — Symmetrical tropospheric scatter path. 


and Rice.* The objective of the present paper is to outline a method for 
predicting the magnitude of intermodulation due to echoes in tropospheric 
scatter radio systems, making use of the intermodulation theory devel- 
oped in Ref. 4 and the reflection theory of tropospheric propagation 
given in Ref. 3. 

The first step in the method, following specification of a suitable model, 
involves calculation of the relative amplitudes of signal echoes received 
from various angles of elevation above the horizon. The corresponding 
time delays of these echoes are then calculated, relative to the signal 
arriving via the antenna beam axes, or “‘centerline”’ path. 

Next, signal-to-intermodulation ratios are determined for selected 
combinations of echo amplitude and corresponding time delay, radio 
system base bandwidth, and frequency deviation. 

These initial, primarily analytical, steps provide an indication of the 
magnitude of intermodulation which would be generated by single echoes 
in an angle-modulated system without diversity. The results postulate 
a so-called “‘worst echo”; that is, an echo having a combination of 
amplitude and time delay that produces more interchannel interference 
than any other combination. 

In calculating the values of echo amplitude and time delay, the 
meteorological conditions of the path are assumed to be the median 
values during the worst month of propagation. Thus, the calculated in- 
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termodulation represents the median value for a single worst echo during 
the worst month of transmission on the given path. 

From this point, the method becomes empirical in accounting for vari- 
ations in the calculated median value of intermodulation as propagation 
conditions change, and in accounting for the effects of multiple echoes 
and diversity reception. Satisfactory analytical methods of handling 
these factors have been elusive, to say the least. However, experimental 
results have provided encouraging means of estimating correction factors 
for their effects, and these corrections are discussed. 


II. MODEL OF TROPOSPHERIC SCATTER PATH 


The typical geometry of a tropospheric scatter path is illustrated in 
Fig. 1, showing two narrow antenna beams intersecting to form a 
“common volume.” The lower edges of the pattern represent the beam 
centerlines; each is pointed toward the horizon on the azimuth of the 
other antenna. The upper edges of the patterns are arbitrarily chosen to 
correspond to some value of antenna gain — e.g., 3 db — less than the 
centerline maximum. Throughout this paper, the antenna pattern above 
the centerline is assumed to be the same as half the free-space pattern. 


2.1 Reflection Theory 


The reflection theory of beyond-the-horizon transmission proposed by 
Iriis, Crawford, and Hogg? has been found to yield estimates of re- 
ceived power which are in good agreement with reported experimental 
data. The expression for received power proposed by this theory will be 
adopted in the present model for use later in deriving an expression for 
echo amplitude. 

Assuming identical transmitting and receiving antennas, and “‘inter- 
mediate-size” reflecting layers,* the received power in a system with 
symmetrical geometry as in Jig. 1 is given by: 


M»> 1 1 a 
roe) 
3 a’ 6a (2+) 6 (1) 


where 
P, = transmitted power, in watts, 
M = a constant of the path proportional to the number of reflecting 
* Ref. 3 proposes somewhat different expressions for the cases of large, inter- 
mediate, and small reflecting layers, where the dimensions are defined in terms of 


wavelength and Fresnel zones. All three may be present on a given path at vari- 
ous times, but intermediate-sized layers are probably most prevalent. 
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layers and the squares of path length, reflecting layer length, 
and change in gradient of refractive index at the bottom of the 
common volume, 

\ = wavelength, meters, 

antenna beamwidth to the half-power points, radians, 

6 = angle between a chord joining the ends of the path and the axis 
of the antenna beam, radians, 

a = half the length of the chord between the ends of the path, in 
meters. For practical paths, this can also be taken as half the 
ground distance between terminals, 


R 
I 


and 


a 

1 1 6 
0 al’ 8 Qa 
ji+5 | I+5 





Equation (1) indicates that the received power is critically dependent 
on the angle @, and it will become apparent later that path intermodula- 
tion also depends heavily on this parameter. Therefore, it is of interest 
to consider those factors which influence the value of 6 on a path of 
given length. These factors — atmospheric refractivity, unequal eleva- 
tion of antenna sites, and nonhorizontal antenna take-off angles — will 
now be discussed. 


2.2 Effect of Refractivity on Angle 6 


As is customary in many radio problems where the transmitted signal 
is refracted by the atmosphere, the lines of propagation in Tig. 1 have 
been drawn straight over an earth having an effective radius greater than 
the true earth’s radius. The symmetrical path is then defined completely 
by the path length LZ (or alternatively, the chord length 2a) and the 
angle @ between the chord and the earth-tangent ray at either end of the 
path. 

The value of @ which gives the correct altitude of intersection for 
horizontally directed antenna beams, allowing for refractive bending, is 


_f_ iL 
~ OR 2RUK 


where F# is the effective earth’s radius, R, is the true earth’s radius (3960 
statute miles), and K is a constant. 

The value K = 4 is used in many radio applications. It has been found 
in the present work that better results are obtained if the value of K is 
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Fig. 2— Minimum monthly mean value of N, (Taken from Ref. 6). 
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determined for the specific path under study. We need, then, to deter- 
mine the relation between K and refractivity. 

A method for determining the approximate value of K for a given path 
is developed in Appendix A. The initial step in the procedure requires 
the selection of appropriate values of minimum monthly mean refrac- 
tivity at sea level, N., for the geographical coordinates of the stations 
at each end of the path. Fig. 2 gives such data, collected from 306 
weather stations over a period of years. 

Next, values of refractivity, N, , are calculated for the altitude at each 
station. This is done using equation (22), repeated here for convenience 
as: 


N, = N. exp (—0.0322h) (3) 


where N, and N, are in ‘“N-units,” and h is the station altitude, in 
thousands of feet. 

Finally, the values of K corresponding to the surface refractivity at 
each station are read from Fig. 3, and a single average for the path is 
estimated from the values at the end points. Unless the path is very long 
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Fig. 3 — Effective earth-radius factor vs refractivity at surface. 
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or there is a large difference between the altitudes of the two stations, 
the values of K determined for the two stations are likely to be almost 
the same. 


2.3 Effect of a Difference in Antenna Site Altitude 


The angle @ for a symmetrical path with no obstructions can be de- 
termined by consideration of refractivity alone, as outlined above. When 
the altitudes of the antenna sites above sea level are significantly dif- 
ferent, or when the antenna centerlines are directed above or below the 
horizontal, corrections are necessary to determine the angles of the beam 
centerlines above the chord. The correction factors may be (and gener- 
ally will be) different at the two ends of the path. The effect of a differ- 
ence in altitude of antenna sites will be considered first. 

Fig. 4(a) illustrates the geometry of an unsymmetrical path. The 
asymmetry here is due to the high altitude above sea level of the antenna 
at the left. For paths of practical interest, the difference in altitude, A, 
between antenna sites is not likely to exceed a few thousand feet. (The 





6,=O+Ph+%, 62 =9-B+6, 
(4, IS NEGATIVE) 


(b) 


Fig. 4— (a) Antenna 1 at altitude A above Antenna 2; (b) negative take-off 
angle at raised site 1, positive take-off angle at site 2. 
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greatest difference in existing systems known to the authors is about 
8000 feet, on a path in the Arctic.) Thus, since A is very small relative 
to the path length, Z, the correction angle 8 to account for the elevation 
difference is given with sufficient accuracy by 


A 
Te (4) 


2.4 Effect of Nonhorizontal Antenna Take-Off Angle 


When an antenna site is elevated and there is no foreground obstruc- 
tion the beam centerline axis is generally depressed from the horizontal 
and aimed at the optical horizon to increase the received signal. Con- 
versely, when there is an obstruction between an antenna and the normal 
horizon it may be necessary to elevate the antenna axis above the hori- 
zontal. 

In the case of Fig. 4(a), for example, the beam from the elevated 
antenna at the left would normally be depressed by an angle 6, in aiming 
at the horizon. This is illustrated in Fig. 4(b). Also in Fig. 4(b), an 
obstruction has been added in front of the antenna at the right, requiring 
that its beam axis be elevated by angle 6. 6; at either end will be con- 
sidered positive when measured upward from the horizontal. 

By combining the effects of antenna altitude difference and nonhori- 
zontal take-off angles, the following general expressions are obtained for 
the angles 6, and 4, between the chord joining the antenna sites and the 
centerlines of the respective antenna beams: 


A; Oa Be Oa (5) 
62 = O6—B+ be (6) 


where @ is defined by (2), 8 is defined by (4), and @ is at the end with 
greater altitude. The factor 6 generally can be ignored if the antenna 
sites are within about a thousand feet of the same altitude above sea level. 


2.5 Concept of ‘Echoes’ 


The conventional diagram for a tropospheric scatter path, as illus- 
trated in Fig. 1, shows well-defined antenna beams projecting from each 
terminal. This is a convenient model for many purposes, particularly for 
determining received signal power. It can be misleading when one is 
concerned about interference effects, however. We will therefore discard 
the notion of a bounded antenna beam, and turn instead to the geometry 
illustrated in Fig. 5. 

In this illustration, the earliest-received — and generally the strongest 
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——— ECHOES 


—— MAIN SIGNAL PATH 
(BEAM AXES) 





Fig. 5 — Illustration of tropospheric transmission in terms of ‘‘main signal” 
and ‘‘echo”’ paths. 


— signal component travels along the axes of the antenna beams and is 
reflected at their intersection in the troposphere. This will be called the 
“main signal.’? Other components of signal are assumed to be smaller in 
amplitude because they are somewhat off the maximum-gain centerline 
of the antenna beam, and also because the reflected power decreases with 
increasing altitude of the reflection point. These reflections will be called 
“echoes” because of their delay behind the main signal. 

Although there must be an infinitely large number of echoes from a 
continuously illuminated common volume, some sections of the volume 
probably produce higher power density than others as seen by the re- 
celving antenna. It is helpful in this argument to consider the total 
received power as some equivalent result of a finite number of signal 
paths. 

In the remaining analytical sections of this paper, it is assumed that 
reflections occur independently, at random, and only in the plane of the 
antenna beam centerlines. Among the entire ensemble of such echoes, 
one is to be designated as the “‘worst’’ echo — the echo calculated to 
produce most distortion in the output of the system. Other echoes — 
including those on either side of the beam centerlines — will contribute 
to the total distortion, but their effects are included later as part of the 
empirical correction for multiple echoes. 


2.6 Summary of Medel 


The results up to this point can be summarized as follows: A path 
structure has been proposed which replaces the conventional 3-db an- 
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tenna beamwidth by a series of signal paths. The shortest of these paths 
is assumed to follow the antenna beam centerlines, and is designated the 
main signal path; the rest are echo paths. Critical parameters in the 
model are the angles 6, and 62 between the antenna beam centerlines and 
the chord joining the two antenna sites. The received power, given by 
(1), depends critically on these angles. When the path geometry is sym- 
metrical, (which would be the case with both antennas at sea level and 
no foreground obstructions), the angles depend only on path length and 
atmospheric refractivity (2). When there is a significant difference in the 
altitudes of the antenna sites or non-horizontal antenna take-off angles, 
6, and 6. must be modified accordingly (5, 6). 

The foregoing model will be used in deriving expressions for echo 
amplitude, echo delay, and finally, intermodulation. 


III. ECHO AMPLITUDE 


The magnitude of an echo relative to the main signal over a tropo- 
spheric scatter path is derived in Appendix B. The derivation makes use 
of (1) to express the power received in vanishingly small beamwidths 
a and ew, representing main signal and echo paths at angles @ and 6. , 
respectively (See Fig. 6). The “reflection loss,” r. , of the received echo 
power relative to the received main signal power is given by (28); thus: 


7 
re (in db) = 10 logy (*) = 10 logo p’. (7) 


It can be shown that (7) is valid for the unsymmetrical path of Fig. 7 





Fig. 6 — Geometry used to calculate ratio of main signal power to echo power 
in a symmetrical path (7). 
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Fig. 7 — General unsymmetrical path. 


as well, where the main signal and echo angles are related by: 
~~ = = p (8) 


Equation (7) also holds for the “small layer’ version of the reflection 
theory. For the “large layer” theory, which may be more applicable 
under high-signal conditions, the reflection loss r, can be shown to vary 
as (p)°. 

The derivation of (7) assumed that the antenna gains in the directions 
of the echo path were equal to the gains along the main signal path. In 
many cases, particularly those involving half-power beamwidths less 
than about one degree, the difference in antenna gain between the angles 
of the echo and main signal paths will add a significant factor to the dif- 
ference between the received echo and main signal powers. Thus, if the 
combined loss due to the decrease in the radiation patterns of the two 
antennas at the angles 6@,, and 62. is given by ra, in db, the resultant 
ratio of main signal to echo power is: 


r(indb) =rn+ ra. (9) 


The reflection loss, r, , is related to p by (7). The antenna loss, ra , can 
be determined with sufficient accuracy from antenna patterns such as 
those illustrated in Fig. 8. The beam axes correspond to angles @, and 
62. Thus, the antenna losses at the echo angles 6, and 62, can be read 
from the appropriate antenna patterns at angles from the axes given by 


Ate == A = (p = 1); (10) 
at one end of the path and 
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HALF-PATTERNS 
ONLY SHOWN aa 
“| 


RELATIVE POWER IN DECIBELS 





ANGLE IN DEGREES FROM BEAM AXIS 


Fig. 8 — Typical free-space patterns of parabolic antennas (750 me). 


O22 — 02 = (p — 1) 6 (11) 
at the other. 


IV. ECHO TIME DELAY 


The time delay of an echo traveling the path 6,., 62. relative to the 
main signal traveling the path 6, , 42 (Fig. 7) is: 


T = 2.68 L 6,02(p’ — 1) microseconds. (12) 


Here the path length, LZ, is in statute miles and p is as defined by (8). 
The derivation of (12) is given in Appendix C. 


V. INTERMODULATION, SINGLE ECHO 


A method for calculating intermodulation distortion caused by a single 
echo in the transmission medium of an angle-modulated system is given 
in Ref. 4. The method assumes a baseband loaded with random noise 
simulating many frequency-division multiplexed voice channels. The 
signal-to-intermodulation ratio, S/J, in a narrow frequency slot any- 
where in the baseband output can be computed for given values of signal- 
to-echo power ratio (r), echo delay (T), base bandwidth (f,)* and rms 
frequency deviation («). The technique is applicable for signal-to-echo 
power ratios greater than about 5 or 6 db. 

*-The theory presented in Ref. 4 assumes that the baseband extends from zero 
to a top frequency f, . In practical systems, the lower edge of the band starts at 


some nonzero frequency, say 12 ke. For systems of many channels, the assump- 
tion of a band extending to zero introduces negligible error. 
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Figure 5.7 in Ref. 4 shows a set of curves from which one can determine 
the S/J ratio in the top channel of an FM system without pre-emphasis. 
Vig. 9 is a similar set of curves computed for a pre-emphasized I'M sys- 
tem. The preemphasis characteristic, Fig. 10, is typical of those used on 
many tropospheric scatter systems. The curves of Fig. 9 are probably 
also sufficiently accurate for phase modulation. 

Fig. 9 is used as follows: 

Given values of base bandwidth, frequency deviation, and an assumed 
value of time delay, compute the rms deviation ratio, o/f,, and the 
product of top baseband frequency and echo delay, f,7 . Note the value 
of the contour line on Fig. 9 at the intersection of the lines for the com- 
puted values of o/f, and f,T’. The value of the contour at the intersection 
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Fig. 9 — Contours of constant interference. (Value in db by which S/I ratio 
exceeds signal-to-echo power ratio in multichannel pre-emphasized FM system. 
Value is for top channel only.) 
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Fig. 10 — Pre-emphasis typical of tropospheric scatter systems. 


represents the amount by which the S/J ratio exceeds the received 
signal-to-echo power ratio, r, in db. Thus, the values of the contour and 
the signal-to-echo power ratio must be added to arrive at the S/Z caused 
by the echo. 

The manner in which S/TJ varies as a function of echo time delay can 
be determined from Fig. 9, for a system of given parameters and a single 
echo of known magnitude. A family of curves of S/J vs. time delay can 
be developed, where each curve in the family applies to an echo of fixed 
amplitude. Fig. 11 shows such a family. The particular case chosen for 
illustration is that of a one-megacycle baseband, pre-emphasized and 
frequency modulated with an rms deviation ratio of one. Each line of 
this family of curves corresponds to a cross-section of the contours of 
Tig. 9 on the vertical line o/f, = 1.0. To the cross-section from Fig. 9 has 
been added, in the case of each curve on Fig. 11, the indicated signal-to- 
echo power ratio 7. 

Note that the family of curves in Fig. 11 applies generally to any sys- 
tem having the selected base bandwidth and deviation. It has not been 
tied to any particular path. The problem now is to find the combination 
of echo amplitude and delay time for a selected path which produces the 
most distortion in the system. To illustrate the method of doing this, 
the following path parameters have been arbitrarily chosen: 

Path length, Z = 200 miles 

Refractivity factor, K = 3 

Antennas, both ends = 60-foot parabolas, with patterns as in Fig. 8 

A symmetrical path 
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Fig. 11 — S/I vs echo delay (single echo) for pre-emphasized FM; f, = 1.0 
mc, o/fe = 1.0. 


Equation (2) is used to calculate 6: = 6 = 0.019 radians. This is the 
angle between the chord separating stations 200 miles apart, and the 
antenna beam axis for a horizontal take-off angle. 

Next, choose an arbitrary value of reflection loss, r.. Suppose 5 db is 
chosen, corresponding to the power ratio 3.16. This is the ratio of the 
reflection loss along some echo path at angle 6,, to the reflection loss 
along the main signal path at angle 6; (or 62). To this value of reflection 
loss must be added the reduction in antenna gain at the angle of the 
echo. 

To find the reduction in antenna gain at the angle 0, corresponding to 
5 db reflection loss, first compute p from (7). Then (10) or (11) can be 
used to compute the angular difference (0. — 61) or (8. — 62). The differ- 
ence in this example turns out to be 0.19°. At this angle from the beam 
axis of a 60-foot antenna, Fig. 8 shows that the antenna loss, r. , for two 
antennas at the angle of the echo path is about one db. The total loss 
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along the echo path, relative to the main signal path, is then r, (5 db) 
plus r. (1 db), or 6 db. 

The time delay of this echo, given by (12), is about 0.075 microseconds. 

Turning now to Fig. 11, the S/J ratio corresponding to the echo delay 
T = 0.075 microseconds and the echo amplitude ratio r = 6 db is about 
35 db. 

Similar computations are now performed for other arbitrarily-selected 
values of reflection loss, say 10, 15, and 20 db. Each case is associated 
with a separate echo path, and the total echo loss and delay along that 
path is computed. The results of such calculations are summarized in 
Table I, together with the S/I ratios that each echo individually would 


TasBLE I — EXAMPLE OF SINGLE-EcHOo S/I CALCULATION 








Reflection Loss. re | Two-Ant. Loss, ra Total Loss, r Delay T S/I 
(db) (db) (db) (usec) (db) 

5 1 6 0.075 35 

10 3 13 0.18 30 

15 8 23 0.32 35 

20 18 38 0.52 50 


produce if it were the only contributor. These values of S/J are also 
plotted on Fig. 11, resulting in the overlay curve showing how S/J varies 
with echo delay and amplitude for the specific path used in the example. 

The important point to note in Fig. 11 is that there is a peak in the 
S/TI intersection curve, occurring at time delays near 0.16 microsecond 
for the 200 mile example path and the assumed baseband and deviation. 
This tendency toward a ‘‘worst echo” effect has been found to be char- 
acteristic of all paths studied. The location of the peak on the time delay 
scale varies depending on the path length, base bandwidth, deviation, 
and other parameters involved. In particular, it has been observed that, 
as path length increases, the peak of the intersection curve for wide 
frequency deviations tends to shift to the region of echo amplitude r less 
than 5 or 6 db. It is in this region that the intermodulation theory loses 
accuracy. This has not appeared to be a serious limitation in work to 
date, however. 


VI. LIMITATIONS OF THE ANALYTICAL METHOD 


The analytical method as outlined predicts the intermodulation in the 
output of a receiver due to a single echo, when the received power is 
assumed given by (1). There are three significant limitations to the 
method at this point: 
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(z) It predicts a single fixed value of intermodulation, whereas meas- 
urements show that there are both long-term changes in the average 
value as well as more rapid variations about an average during a short 
interval. 

(iz) The analysis predicts the effects of single echoes only; practically, 
there must be a continuum of echoes spread over something vaguely 
defined as a “common volume,” and the total intermodulation gener- 
ated by these echoes is bound to be greater than that of a single ‘‘worst”’ 
echo. 

(477) The theory of intermodulation due to transmission echoes as 
presented in Ref. 4 is valid only for echo amplitudes smaller than about 
5 or 6 db relative to the main signal. On tropospheric paths, theory sug- 
gests that the strongest echoes are normally associated with the shortest 
time delays. Having short time delays, they would not contribute greatly 
to interchannel modulation, although they would be responsible for signal 
fading and the associated increases in fluctuation noise during fades. 
This is not the entire story, however. While strong echoes may normally 
be associated with short time delays, there are surely some meteorological 
conditions which produce strong reflecting or scattering regions at ab- 
normally high altitudes. Occasionally, then, an echo might be compar- 
able in magnitude to the main signal and be delayed enough to cause 
severe distortion. It should be emphasized that this apparently is not 
the average situation, which is of primary concern in this paper. 

Kefforts to develop satisfactory analytical techniques to remove the 
limitations noted above have thus far been unsuccessful. Empirical cor- 
rections have therefore been developed to account for them, based on the 
comparison of calculated values and results of tests presented in the 
following section. 


VII. COMPARISON OF CALCULATED AND MEASURED INTERMODULATION 


Measurements of intermodulation have been made on four tropo- 
spheric scatter paths with diverse climates, path lengths, and siting con- 
ditions. Significant parameters of the four paths are given in Table II. 

The four systems were Gaussian-noise loaded except for a slot near the 
top of the baseband. Median values of interference power were measured 
in this slot at the receiver during approximately 14 to 3-minute samples. 
The median of each sample was adjusted to account for fluctuation 
noise and equipment intermodulation, leaving the value of intermodula- 
tion noise chargeable to the path. A more detailed account of test pro- 
cedures and methods of data analysis is given in Ref. 7. 

Table III shows a comparison of the raw measured non-diversity 
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TasBLe IT — PARAMETERS OF Four Test PATHS 








Path 1 Path 2 Path 3 Path 4 





Trans. | Rec. Trans. Rec. Trans. Rec. Trans, Rec. 











Location Caribbean Arctic Arctic Arctic 
(over (over ice) (over water) (over water) 
water) 


Length, Stat- 185 194 440 340 
ute Miles 





Radio Freq., 725 900 900 800 
me (Nomi- 
nal) 











Ant. Dia., ft. 60 60 30 60 120 120 120 120 














Site Elev., ft. | 200 10 8000 1100 70 1100 700 1100 
(Above Sea 
Level) 





Eff. Earth 1.55 1.26 1.37 1.38 
Radius 
Factor, K 


Ant. Takeoff 0 =0 /—0.0035}|—0.004 |—0.010 |—0.0026|—0.0005)—0.0014 
x; O:, Ta- 
dians 


X%’s Between |0.0151/0.0151| 0.0227] 0.0087} 0.0375] 0.0309} 0.0295} 0.0308 
Ant. Axes 
and Chord, 
6, and 62, 
Radians 














results from the four test paths and the results calculated for the single 
“‘worst”’ echo. The same data are presented as a scatter plot in Fig. 12. 
With the exception of the Path 2 data, each point represents an average 
of the median values from at least 5 to as many as 20 test samples, for 
one combination of base bandwidth and frequency deviation. Only 3 
samples were available for each combination from the Path 2 tests. 


VIII. ADJUSTMENT OF MEASURED INTERMODULATION TO WORST-MONTH 
PROPAGATION CONDITIONS 


Results from the four tropospheric scatter test paths have shown that 
there is a decided correlation between average intermodulation and 
average path loss (or received power). As path loss increases, intermodu- 
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TABLE III — ComPaARISON OF SINGLE-ECHO CALCULATIONS AND 
Raw Measurep S/J, DB 





Path 1 Path 2 Path 3 Path 4 
Top Baseband RMS Dev., 
Freq., kc kc* 





Calc. | Meas. Calc. | Meas. | Calc. | Meas. | Calc. | Meas. 























108 62 — _ ae pated 51 38 — — 
125 — — 68 48 45 35 51 47 

250 — — 62 43 40 31 45 43 

500 “— oo 57 37 36 30 41 38 

800 = _ 54 33 36 32 — — 

1000 = = — — — ae: 40 33 

300 62 = _ — — 43 32 —_ — 
125 — — 59 43 37 28 42 41 

250 55 48 53 38 32 27 37 37 

500 49 41 48 33 27 27 32 32 

800 = a 45 30 26 27 — — 

1000 45 36 = _ — oe 30 28 

2000 43 34 — = — — — — 

552 62 — — _ = 37 31 — — 
125 — = 53 36 31 27 37 35 

250 50 42 48 33 26 24 32 32 

500 44 36 43 29 22 24 27 27 

800 —— = 40 27 21 23 — — 

1000 40 29 _ — — = 25 23 

2000 37 28 — = — — — — 

1052 125 = — 49 30 — ms _— _— 
250 44 37 43 27 = _ 33 27 

500 39 32 37 23 — = 27 25 

800 = cam 34 21 = — 22 22 

1000 35 26 — _ = _ — — 

2000 32 24 — = = = 20 18 


* The deviations used on Path 3 were slightly smaller than listed here. In each 
case, however, calculated values of S/J are given for the actual deviations used 
in tests. 


lation also increases. While the relationship between path loss and 
intermodulation may not be highly obvious from the various equations 
presented to this point, the fact that intermodulation does vary and 
that a correlation with path loss shows up in measured data leads to the 
following questions: 

(t) To what value of path loss does the computed value of intermodu- 
lation correspond? In other words, what is a useful point of reference for 
the observed dependence? 

(it) What is the rate of variation of intermodulation with path loss? 

The point of reference for the relation between intermodulation and 
path loss has already been suggested in the earlier discussion of refrac- 
tivity. It was indicated there that in calculating the effective earth radius 
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Fig. 12 — Calculated single-echo S/Z vs average measured (unadjusted) S/T. 


factor, K, the minimum monthly mean value of refractivity indicated in 
Fig. 2 should be used. Minimum refractivity corresponds to minimum 
bending of the radio waves, maximum values of antenna centerline 
angles 6, and 62, maximum time delay for an echo of given amplitude, 
and therefore maximum intermodulation. Correspondingly, (1) indicates 
that minimum refractivity (maximum 6) tends to minimize the average 
received power. Thus, the observed direction of variation of intermodula- 
tion and path loss is borne out by our model. It may be deduced that the 
value of intermodulation computed on the basis of minimum monthly 
mean refractivity corresponds directly to what is commonly called 
“worst-month” path loss —the monthly median path loss during the 
worst, month-long period of propagation during the year. (There should 
be some probability attached to this value, but it is usually approxi- 
mated without this refinement.) Estimates of worst-month loss on a 
new path are generally derived from data from other paths, with subse- 
quent adjustment based on measurements over the new path. 

It can be noted that the choice of worst-month conditions as the point 
of reference is somewhat arbitrary. We might have chosen instead to use 
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an average annual value of refractivity, compute an average annual value 
of intermodulation, and relate the latter to the average annual path loss. 

Apart from the point of reference chosen, there remains the second 
question posed above: what is the rate of variation of intermodulation 
with path loss, for propagation conditions other than the reference condi- 
tion? If data were available on the time variations of refractivity for 
any path of interest, it might be possible to compute an approximate 
relationship. Unfortunately, such detailed refractivity data are not 
generally available. As a substitute, we are forced to use empirical ad- 
justments based on results of measurements of intermodulation and 
received power. 

Fig. 13 is typical of the data which have been obtained. Shown are 
two plots of intermodulation (in terms of S/J ratio) versus received 
signal measured on Path 3 — a 440-mile over-water path in the Arctic. 
Each point represents the median values of the two variables during a 
test sample of about 90 seconds. The total test period covered about ten 
days, during which time a wide range of propagation conditions was 
encountered. The rather wide scatter of points is an indication that 
factors other than average path loss can affect the value of intermodula- 
tion during any specific test sample. For example, a localized region of 
high reflecting ability could produce a strong echo and unusually high 
intermodulation. 

The data of Fig. 13 apply for a base bandwidth of 300 ke, rms fre- 
quency deviations of 125 and 250 ke, and no receiver diversity. The 
slope of the curve drawn through the 125 ke deviation scatter plot is 
slightly greater than one; that is, intermodulation and path loss varied 
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Fig. 13 — S/I vs path loss — 440-mi arctic path, 300-ke baseband, no diversity. 
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approximately db for db. The slope of the 250 ke plot is slightly more 
than 0.8. The least squares method was used to calculate these slopes. 

Similar scatter diagrams plotted for other combinations of these trans- 
mission parameters and the other three paths indicate that the change 
in intermodulation ranged from about 0.4 db to 1.2 db per db of path 
loss change. The lowest values were noted on Path 1 — a 185-mile path 
between Florida and Cuba. There was no apparent pattern to the slopes 
of the scatter plots as a function of baseband or deviation on any of the 
four paths tested, except that the slope was generally lowest for the 
highest deviation ratio. There is a technical justification for the latter 
effect. contained in the contour curves of Fig. 9. Note that a variation 
in time delay of an echo of given amplitude represents a change in the 
ordinate f,7T on these curves. It is clear that for a given change of time 
delay (say, corresponding to a change in f,7' from 0.1 to 0.2), a system 
operating at a high deviation ratio (toward the right on the o/f, scale) 
will in general exhibit a smaller change in S/J than will a system operat- 
ing at a low deviation ratio. Except for this trend, however, the effects 
of different path lengths, antenna sizes, take-off angles and climates 
confound any attempt to find a precise dependence of S/J on path loss. 

The data referred to above, together with similar data from three 
other paths, indicate that an adjustment to the computed value of S/J 
of about 0.7 db per db of path loss change is reasonable when one is try- 
ing to predict the probable average value of S/J for a path at times other 
than the worst month. If one is trying to compare a computed worst- 
month value of S/J with results of measurements on a specific path, it 
will be more accurate to determine the average slope of scatter plots 
similar to those of Fig. 13, and use that slope to place the computed and 
measured results on a comparable path loss basis. Thus, the correction to 
be added to the median measured S/J during a test period to adjust it to 
worst-month conditions can be expressed as: 


D = S(TPL — WPL) (13) 


where S is the average slope, 7,PL is median path loss during the test 
period, and WPL is the worst-month median path loss. 

During the test periods for which the results in Table III and Fig. 12 
apply, the median path losses on Paths 1, 3 and 4 were four to eight db 
less than worst-month estimates. During the tests on Path 2, which were 
conducted over a period of 10 days, the median path loss was about 8 db 
greater than the worst-month median estimate. Thus, adjustments to the 
measured data are necessary to convert them to propagation conditions 
comparable to those assumed in calculations. This was done as follows: 
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For each path, the average slope was determined of the least-squares 
lines drawn through scatter plots of S/J versus path loss (similar to those 
illustrated in Fig. 18). The average slope was then used to adjust the 
measured S/J results to the worst-month path loss condition, using 
(13). The estimated worst-month median path losses, test period median 
path losses, and the slopes used for the S/J adjustments for each of the 
four test paths are given in Table IV. 








TaBLE IV 
Path 1 Path 2 Path 3 Path 4 
Worst Month Median Path 197 db 213 db 243 db 240 db 
Loss (WPL) 
Test Period Median Path 193 db 221 db 235 db 236 db 
Loss (TPL) 
Avg. Slope of S/I vs Path 0.5 0.7 0.9 0.6 
Loss 
Total Adjustment to Meas- —2 db +6 db —7 db —2 db 


ured S/I Data 


When the adjustments indicated in the last line of Table IV are made 
to the measured data plotted in Fig. 12, the result is as given in Fig. 14. 
With the exception of the Path 3 data, the adjusted points exhibit a 
noticeably narrower spread than before the path loss adjustments were 
made. However, the group of points as a whole still differs significantly 
from the calculated results for single echoes. 


IX. ADJUSTMENT OF CALCULATED INTERMODULATION TO ACCOUNT FOR 
MULTIPLE ECHOES 


It has been emphasized that a calculated value of intermodulation 
represents only the contribution of the single worst echo among echoes 
which have at least 5 or 6 db lower amplitude than the first-arriving 
signal element. In fact, however, there are probably very many echoes 
present at any given time on a tropospheric scatter path. The combined 
effects of a multiplicity of echoes is to produce more intermodulation 
than the maximum attributed to any single echo. Thus, it seems qualita- 
tively reasonable to attribute the remaining differences between calcu- 
lated and measured results in Fig. 14 to the effects of multiple echoes. 

The results for Paths 1 and 2 in Fig. 14 appear to be closely intermin- 
gled, and they differ from the calculated values for single echoes by an 
average of about 10 db. The Path 4 results form a group somewhat closer 
to calculated values; their average difference is only about 4 db. As al- 
ready noted, the differences for Path 3 vary more widely. This is thought 
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Fig. 14 — Calculated single-echo S/Z vs average measured S/I normalized to 
worst-month path loss conditions. 


to be at least partially the result of a propagation anomaly observed on 
this path, the effects of which will be discussed shortly. There are also a 
number of other recognized sources of possible errors in both the calcu- 
lated and measured data, and these will also be discussed. Since there is 
no way at this time to resolve the quantitative effects of these errors, we 
are left with the undesirable but rather inescapable alternative of trying 
to estimate an adjustment for multiple echoes from the scatter of points 
as they stand in Fig. 14. 

Fig. 15 shows the result of subtracting nine db from each point in Fig. 
14. This value gives a little extra weight to the data from Paths 1, 2, and 
3 for the following reasons: (a) there is relatively good agreement among 
the adjusted data for these three paths, and (b) up to the time of this 
writing, there has been considerable uncertainty as to the path loss data 
which should be assumed for Path 4. 

While 9 db appears to be a reasonable average adjustment, it should 
not be expected that the intermodulation measured during any single 
test sample would be 9 db more than the calculated figure. The number 
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Fig. 15 — Calculated S/Z adjusted for multiple echoes vs average measured 
S/I normalized to worst-month propagation. 


of significant echoes is bound to vary from test sample to test sample as 
atmospheric conditions change. We are aware, too, of various sources of 
error in calculated or measured results. These errors are discussed in a 
later section, and it will be apparent that some of them can have an ef- 
fect on the values assumed for path loss and multiple echo adjustments. 

With the caution that some change in the number can be expected as 
better understanding and more and better data become available, it is 
felt that nine db is a reasonable number to account for the effects of 
multiple echoes in calculations of path intermodulation. 


X. ANOMALY ON PATH 3 


Fig. 15 shows that there is a wider variation of differences between 
measured and calculated results for Path 3 than for the other paths. 
This might be explained by the fact that the extreme length of the path 
makes it possible for strong reflections from the tropopause to cause more 
interference, at times, than echoes in the troposphere alone. The altitude 
of the tropopause over Path 3 was reported at 25,000 to 35,000 feet from 
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day to day during the test period. When the tropopause is in the lower 
half of this range, according to the method of calculation described earlier 
in this paper, reflections could cause distortion that is significant with 
respect to the maximum from a single echo in the troposphere. The effect 
is greatest for low FM deviations, as indicated by the calculated results 
in Table V. 

Fig. 16 shows the same comparison as Fig. 15 except that on Path 3, 
a reflection is assumed to occur at 27,000 feet with no more reflection 
loss than the main signal. Thus, the echo at the receiver would be reduced 
below the main signal only by the antenna losses due to the elevated 
angles of the path to the tropopause. In plotting Fig. 16, the interference 
from this echo was added to the single echo maximum before making the 
9 db correction for multiple echoes. 

To make the calculated S/J for Path 3 correspond even better with 
the measured values, it would be necessary to weigh the effect of re- 
flections from the tropopause more heavily than the effects of the tropo- 
spheric reflections. One way to increase the effect of the tropopause 
reflections would be to assume its altitude to be 25,000 or 26,000 feet. 
Another possibility is to assume that the tropopause provides a better 
reflecting surface than do the layers in the troposphere, so that, except 
for the antenna loss, the echo from the tropopause would be stronger 
than the main signal. 

Still another possibility is contained in the following argument. In 


TasBLe V — Errect or TropoprausE on S/J Ratios, Pato 3 








Calculated S/I, db Measured S/I 
Top Baseband RMS Dev. Tro 
? po- Tropo- Sum of 9-db : 
Frequency, ke ke sphere pause First Cor: Adj, to 


r 

ra Raw Worst 
cho, Two Multiple 

Echo 27,000’ Echoes Echoes Month 











108 62 51 49 47 38 38 31 
125 45 44 41 32 35 28 
250 40 41 37 28 31 24 
500 36 44 30 26 30 23 
700 36 50 36 27 32 25 
300 62 43 42 39 30 32 25 
125 37 36 33 24 28 21 
250 32 33 29 20 27 20 
500 27 36 27 18 27 20 
800 26 41 26 17 27 20 
552 62 37 38 34 25 31 24 
125 31 33 29 20 27 20 
250 26 29 24 15 24 17 
500 22 31 22 13 24 17 
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Fig. 16 — Calculated vs measured S/J when tropopause reflection is added to 
calculated results for path 3. 


applying the model of the tropospheric scatter path described in this 
article, the power in the lower half of the antenna beams has been neg- 
lected as though it were blocked by the horizon. Suppose that power in 
the lower half of the antenna beam grazing the surface of the sea were 
reflected (and spread somewhat by the convexity of the surface) with 
negligible loss. This power, with almost no antenna loss, may be reflected 
again from the tropopause and cause high distortion because of its long 
delay and low loss. This argument is particularly suited to Path 3 where 
the antennas at one end are located on a cliff 1100 feet above the sea. 

There seems little point in speculating further about the reasons for 
the relatively peculiar behavior of the results for this path. While each of 
the arguments above has some plausibility, the fact remains that infor- 
mation is not now available to prove any of them. 


XI, SOURCES OF ERRORS IN CALCULATED AND MEASURED DATA 


The foregoing comparison of calculated and measured results indicates 
that the analytical-empirical method of predicting intermodulation has a 
standard error of estimate of about 2.6 db relative to measured results. 
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Part of this error is due to the limited accuracy of the empirical adjust- 
ments for path loss and multiple echoes. There are in addition a number 
of other possible sources of error in both the calculated and measured 
data, and these will now be discussed. 


11.1 Path Loss Data 


Estimates of worst-month path loss are difficult to make, and errors 
of 2 or 3 db are not unlikely. In fact, different methods of estimating 
path loss show that much discrepancy. 

Determination of the median path loss from tests is also subject to 
errors. Knowledge of the values of transmitter output power, received 
signal power, RF plumbing losses, and antenna gains is required, and in- 
evitably some error creeps into measurements of these parameters. Any 
error in such measurements should be approximately constant (assuming 
negligible changes between test periods) for a given path. Hence, it 
should not affect the slope of an S/J versus path loss scatter plot for 
that path. However, a constant error would affect the value of the ad- 
justment used for multiple echoes. 


11.2 Antenna Orientation and Pattern 


Ideal free-space antenna patterns and centerlines directed at the 
horizon are assumed in calculating intermodulation by the method out- 
lined in this paper. In practice, local siting conditions may have some 
effect on an antenna pattern and modify the assumed amplitude rela- 
tionship between echoes at different angles from the horizon. The errors 
due to these factors should be a constant for any one path, but when 
the results for several paths are brought together as in Figs. 13 to 16, 
they will contribute to the spread of data. 

The antenna patterns assumed in all the calculated results in this 
paper were those shown in Fig. 8. These patterns — determined from 
scale-model measurements — apply at a frequency of about 750 me. As 
Table II indicates, the frequencies on the test paths varied from about 
725 me to 900 me. In view of the uncertainty noted above relative to 
effects of local siting conditions, it has not appeared worthwhile to adjust 
the nominal patterns for 750 mc to patterns for specific test frequencies; 
however, it is worth pointing out this factor as a possible contributor to 
the spread of data. 


11.3 Extrapolation of S/I vs Time Delay Curves 


In the example used to illustrate the calculation of single-echo inter- 
modulation, the peak of the S/J versus time delay curve (Fig. 11) was 
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quite pronounced. This example assumed a path length of 200 miles and 
an rms deviation ratio of unity. lor longer paths and higher deviation 
ratios, it is sometimes found that the highest value of S/Z which can be 
calculated falls on the 6-db echo amplitude curve. This indicates that 
the peak of the curve may be further to the left of the 6-db contour, or in 
the region where the intermodulation theory of Ref. 4 is not valid. Al- 
though we have taken the liberty of extrapolating into this region to 
obtain certain of the values used in Figs. 13 to 16, it is clear that this is 
another source of possible error. 

It is worth pointing out that in those cases where such extrapolation 
has been necessary to arrive at an estimate of intermodulation, the 
resulting value of S/J is poor enough to render the system rather useless. 
In a sense, then, the fact that there may be an error due to the extrapola- 
tion is academic. 


XII. EFFECT OF DIVERSITY ON INTERMODULATION 


The calculated results have been compared here to the S/J ratios 
measured without diversity. No approach has been found to date which 
permits analytical prediction of the diversity improvement; however, 
experimental results give some indication of the advantage to be ex- 
pected. 

Measurements of S/J were made with dual and quadruple diversity 
reception as well as without diversity on Paths 2, 3, and 4, and with 
dual diversity and nondiversity on Path 1. Generally, the difference in 
S/I between nondiversity and dual, and between dual diversity and 
quadruple, was 2 to 5 db. On the three Arctic paths, though, it is ques- 
tionable whether the full advantage of diversity was being realized be- 
cause the combiners were not working as well as expected, and there was 
some unexpected interference from testing equipment. 

Dual diversity results for Path 1 showed an improvement in inter- 
modulation over nondiversity during high deviations, but the non- 
diversity results were about equal to dual diversity results for low 
deviations. 

To generalize, on the basis of measured data, one might expect to 
realize 3 or 4 db improvement in median S/J ratio with dual diversity 
over nondiversity reception, and another 3 or 4 db improvement for 
quadruple diversity over dual diversity. 


XIII. SUMMARY OF ANALYTICAL-EMPIRICAL METHOD 


The relatively good agreement between calculated and measured re- 
sults indicates that the approach described in this paper is some progress 
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in the right direction. Further measurements and more accurate meth- 
ods of measurement may show the way to refinements in the analysis 
which will improve its accuracy and minimize the present reliance on 
empirical adjustments. 

There is given below a step-by-step summary of the method suggested 
for predicting the average intermodulation due to the propagation me- 
dium in a tropospheric scatter system. While the sequence is rather long 
and involved, many of the steps can be formulated in graphical or chart 
form for general use in a variety of cases. In fact, a computer program 
has been developed to accomplish steps 2 through 11, using as inputs the 
parameters listed in step 1. 

The steps involved in the method are as follows: 

1. Determine these parameters of the system and the path: 

a) Path length, statute miles 

b) Antenna beam pattern for the appropriate antenna size and 
operating frequency 

c) Antenna site elevation, feet, at each end 

d) Antenna take-off angles, 6 and 6 

e) Minimum monthly median values of radio refractivity at sea 
level, N., in the areas of the two antenna terminals. (From 
Fig. 2.) 

2. Calculate the surface values of refractivity, NV; , using (3) and the 
values of N. found in step 1(e). 

3. From Fig. 3, determine the values of the earth radius factor K 

’ corresponding to the values of N, found in step 2. Estimate an 
average value of K for the path. 

4. Compute the base angle, 6, defined by Fig. 1 and (2). 

5. Compute the correction angle 8, (4), to account for the difference 
in antenna site altitudes. 

6. Compute the total angles 6, and 6, between the chord joining the 
antenna sites and the centerlines of the respective antenna beams, 
using (5) and (6). 

7. Choose arbitrarily four or five values of reflection loss, r. , spaced 
through the range of approximately 5 to 30 db. 

8. Compute the value of p corresponding to each value of reflection 
loss selected in step 7, using (7). 

9. Determine the additional loss of the echo path due to the antenna 
patterns, for each echo selected in step 7. Use (10) and (11) to 
determine the echo angles, and the antenna patterns of Fig. 8 to 
determine the reduction of gain at those angles. 

10. Calculate the total loss for each echo, using (9). 
11. Compute the time delay, 7’, for each echo, using (12). 
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12. For given base bandwidth and frequency deviation, use Fig. 9 to 
determine the value of S/J for each echo. (For a straight-FM sys- 
tem, use Figure 5.7 in Ref. 4.) 

13. Plot S/Z versus time delay for each echo, and determine the mini- 
mum value of S/J. Alternatively, the minimum can usually be 
estimated with sufficient accuracy from inspection of the calcu- 
lated S/I values for the several echoes. 

14. To account for multiple echoes, subtract 9 db from the minimum 
value of S/I found in Step 13. 

The result of step 14 is the predicted value of the monthly mean S/I 
ratio during worst-month propagation conditions. The result applies to 
the top channel of a pre-emphasized I'M system operating without 
diversity. For dual or quadruple diversity systems, add 3 db or 6 db, 
respectively. For other than worst-month propagation, add 0.7 db to the 
S/T ratio for each one-db decrease in path loss. (The correction factor 
may be slightly higher or lower than 0.7 db for a specific path, and should 
of course be modified appropriately if the S/Z versus path loss function 
has been measured.) 

The geometry of some paths may present unusual propagation char- 
acteristics calling for special treatment beyond the steps outlined above. 
The case of the 440-mile path discussed earlier is an example; on this 
path, it appeared that reflections from the tropopause may have con- 
tributed significantly to the intermodulation in the system. The implica- 
tion of this case is that propagation anomalies can be expected on paths 
several hundred miles long, where the common volume extends to very 
high altitudes. 
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APPENDIX A 


Effect of Refractivity on Effective Earth Radius 


We shall outline here a method for determining the approximate 
value of the effective earth’s radius factor, K, applicable to a selected 
tropospheric scatter path. 

K is expressed with sufficient accuracy for our purpose by® 

R 1 
| le ee 
R, dn (14) 
1+ dh R. 
where dn/dh is the vertical gradient of refractive index, n. In turn, define 
refractivity: 


N = (n — 1)10¢ (15) 
Within a kilometer of the earth’s surface, 

dn As AN —6 

TY 0”? (16) 
where 


N, is the refractivity at the surface, and N;, is the value at altitude h. 
When h = 1 km, it has been found from extensive data’ that 


ANixm = —7.32 exp (0.005577 N,). (18) 


One useful model of the atmosphere assumes an exponential decrease 
in N with altitude. Then: 


N = N, exp (—Ch) (19) 


where h is altitude in kilometers above the ground. To fit the data quoted 
above, at h = 1 km, 


N; 


a ANA 


(20) 

Considering dn/dh to be constant with altitude (a good approximation 
only for the first 4000 or 5000 feet, but used here throughout because it 
introduces no appreciable error in calculating K for even the longest 
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path considered) and approximating it by the change in refractivity for 
the first 1000 feet only, 


AN _ —N, 
1000 ft. 1000 


Substituting (16), (18), (20) and (21) in (14), K is found to be a 
function of NV, . This is plotted in Fig. 3. 

Since the quantity to be calculated later is intermodulation during 
worst-month conditions, the K of interest is the value for the worst 
month. Maximum intermodulation generally corresponds to the mini- 
mum value of K, since it will produce least bending and hence longest 
echo delays. Minimum K corresponds to minimum N,. N, may be 
calculated from N,, the sea level value of N, by 


N,; = N. exp (—0.0322 h) (22) 
where h here is the station altitude in thousands of feet. Minimum values 


of monthly mean JN, are given in Fig. 2, based on data collected from 
306 weather stations over a period of several years.® 


[1 — exp (—0.3805 C)). (21) 


APPENDIX B 


Calculation of Echo Amplitude Relative to Main Signal Amplitude 


The technique involved here to calculate the magnitude of trans- 
mission echoes in a tropospheric scatter path makes use of (1) for re- 
ceived power. It will be convenient and helpful in this derivation to refer 
also to lig. 6, which shows symmetrical main signal and echo paths 
having equal beamwidths a and a, , respectively, and directed at angles 
6 and @, , respectively. 

Equation (1) can be used to calculate the relative power received from 
the two paths, assuming the antenna gains in the directions @ and 6, to 
be equal. (The effect of differences in antenna gain at 6 and 6, is discussed 
in the main text.) Expressing the ratio of the power of the main signal 
to the power of the echo by 7’, we have, from equation (1), 


oe © 2+) = 


we +a)s(@) 


Setting a, = a and letting a — 0, the limit, r, , will be the ratio between 
the power of a signal component traveling the centerline path and an 
echo traveling some longer path defined by @, . 
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Let f(a/@) be written as 





O° (26 + *] 
= — , 2 
oe E *@+a 80+ a) ve 
Then 
: ’ 0," : g(0 a) 
= —— nes : . 2 
Pe Re gee) 7 
The above limit is undefined, and so is 
lim / (Gye) (26) 





a->0 g' (Be ’ a) ; 
where g’ = dg/da. Using L’ Hospital’s rule with the second derivative, 


lini go" (Byo) () (27) 


a0 g” (8c, a) = 8 
or, 
0.\" 
n= (%) = oy (28) 


r, will be called the “reflection loss” of an echo relative to the main 
signal. 


APPENDIX C 


Calculation of Echo Time Delay 


The time delay of an echo which travels the path defined by 1. , 42 
instead of the path 4, , 6. (Fig. 7) is 
AL 
= 300) seconds, (29) 
where 
ALS (Lie + Line) — (Ly aa Iz) meters. (30) 
By the law of sines, 
2a (sin A. + sin 0) 2a (sin 6, + sin 6) 
sin (Ore + Boe) sin (61 -+- 42) 
22 2a (sin O16 + sin Boe) (31) 
SIN O1¢ COS Axe -+ SIN Oe COS Are 
x 2a (sin 6, + sin 62) 
sin 6, cos 62 + sin 02 cos 01 


AL = 





36 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


For angles 6 less than 0.1 radian* the following approximations are valid, 
applied to the above equation: 





2 
cos@= 1 — 2 
sin 0 = 0 
1 9 
eee ss 2 
2 
2a=L 
Then 
AL = 5 (elie — 0402) meters. (32) 
Finally, for Z in statute miles, 
T = 2.68L 6,0:(p' — 1) microseconds (33) 
where 
Bie One 
oar an 
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A Wideband Transistor IF Amplifier for 
Space and Terrestrial Repeaters Using 
Grounded-Base Transformer- 
Coupled Stages 
By W. F. BODTMANN and C. L. RUTHROFF 


(Manuscript received April 27, 1962) 


A wideband transistor IF amplifier is described. A measured model has 
a gain of 41 db and is flat to +0.1 db from 50 to 100 mc. Data on noise 
figure, effect of temperature changes and power output are presented. An 
electronically variable attenuator is described which is suitable for use in 
such an amplifier as part of the AGC circutt. 


I. INTRODUCTION 


The general objectives of the work discussed in this paper are as 
follows: 

(i) To design a wideband IF amplifier suitable for use on heavy-route, 
long-haul microwave radio relay systems.!? 

(ii) To design a wideband IF amplifier with low power consumption 
and good temperature stability for use in microwave radio relay sys- 
tems? which may be battery powered and which must operate at out- 
door ambient temperatures. 

In accordance with these objectives it was decided to design an ampli- 
fier to be substantially flat over a 40-me band centered at 70 mc. Tran- 
sistors were chosen as the active elements because of the low bias power 
required, and the common-base circuit configuration was selected be- 
cause it appeared likely to result in a stable band shape with respect to 
temperature. 

A part of the initial philosophy was to start with a stable wideband 
amplifier. The desired frequency selective circuits could then be added 
as linear lumped circuits whose gain and phase properties are in good 
control. Such circuits could include any desired gain and phase equaliza- 
tion. This approach results in an IF amplifier with an excellent delay 
characteristic. 
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WHERE: 
Zin © Yet Zpli-e) 

2. ZyparptJwl, 
Z.(I-@& )>> 2, 
Z.>>Z, + Zp 


Rin 


—* 
00Q 


(d) 


Fig. 1 — (a) Typical interstage circuit, (b) equivalent circuit for the inter- 
stage, (c) equivalent circuit for computation of input impedance, (d) equivalent 
interstage circuit including losses. 
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Il. THE IF crrcurirT®6 


The common-base configuration was chosen for several reasons. The 
transistors available have alpha cutoff frequencies in the region 500-700 
me. If these transistors were used in the common-emitter configuration 
in this application, the IF band would extend well beyond the beta cut- 
off frequency and considerable equalization would be necessary to shape 
the IF characteristic. Furthermore, changes in ambient temperature 
result in a shift in both alpha and beta cutoff frequencies. The alpha 
cutoff is well beyond the IF band but the beta cutoff is not. Hence, for 
operation over a reasonably large temperature range, the common-base 
stage has better stability with respect to band shape and midband gain. 

Simple broadband transformers’ are available which can be used in 
the common-base configuration and result in an interstage with at most 
one tuning element per stage. 

The interstage circuit is shown in Fig. 1(a). Z is the stray inductance, 
Cr the total interstage capacitance, ZL, is an inductance added in the 
base of transistor Q2, Ti is a wideband transformer and RA; has been 
added to help shape the interstage characteristic. This circuit has been 
redrawn in Fig. 1(b), where the assumptions have been made that the 
output of transistor Q; is a constant current and 7 is an ideal trans- 
former with an impedance ratio of 4:1. Transistor Q2 has been replaced 
by its input impedance, which is considered to be Riy + jwLry where 
both Rr and Lry are functions of frequency. The blocking capacitors 
and bias circuits have been omitted. From Fig. 1(b) we see that the out- 
put capacitance of transistor Q,, the input impedance of transistor Qe , 
stray inductance L, and resistance R; form a peaking circuit. The equiv- 
alent circuit used to compute the input impedance of Q. is shown in 
Fig. 1(c). It has been assumed that Z.(1 — a) > Zrand Z. >> Zi + Zp. 
In the Appendix it is shown that the frequency response of Fig. 1(b) is 








Io _ (Rr — nw)" + (uy + Xt — Xe)" (1) 
IT, XX? 
where Rp = A tret+r 
w= (re + yols)/(L +7) 
Y= S/fa 
re = emitter resistance 
7‘, = high-frequency base resistance. 


The high-frequency response of this circuit can be controlled by ad- 
justing R, or L,. Amplifiers have been built using both schemes. The 
effect of changing R; is to change the Q of the circuit. The effect of adding 
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a suitable inductance, LZ, in the base is a little more complicated but 
results in a better frequency response. 

An inductance in series with the base of a transistor changes the 
transistor input impedance by decreasing the resistive part and increas- 
ing the reactive part. (See Appendix.) Thus, instead of adjusting R, , 
a variable inductance in series with the base can be used to tune the 
interstage. Even though this method of tuning adds inductance to the 
emitter circuit it does not decrease the bandwidth of that portion of the 
response which is flat to a few tenths of a decibel. In fact if R, and Ly 
are properly chosen, the amplifier will be compensated over a wide 
frequency range. This is demonstrated in Figs. 2, 3, 4. Fig. 2(a) shows 
the measured responses of the peaking circuit with and without added 
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Fig. 2 — (a) Measured responses of series peaking circuit, (b) circuit used to 
measure interstage series peaking response. 
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Fig. 3 — Amplitude-frequency response with no added inductance in base of 
transistor. 


inductance in the base circuit. The circuit measured is shown in Tig. 
2(b). The computed curves are identical to the measured ones when the 
following parameters are used in (1): 


fa = 600 mc Cr = 4Supf Ves = —10v 
L = 0.02uh Ly = 0.05uh Iz = 10 ma 
ry = 50 ohms re = 3 ohms 

Rr = 115 ohms when Ly, = 0.05uh 


Rr = 110 ohms when = 0. 


Fig. 3 shows the computed response of a single stage and computed 
and measured responses for eight identical interstages, with no added 
inductance in the base lead. The transistor parameters are as listed 
above. I, = 0 and R, = 54 ohms. For the computed responses a has 
been assumed equal to one and the transformers ideal. Since these 
assumptions are not exactly true, the measured response can be expected 
to be poorer than the computed response. The midband insertion gain 
between 75-ohm impedances is 41 db. 

Fig. 4 shows the computed responses of a single interstage with two 
bias conditions and the measured response of eight interstages with 
inductance added in the base lead. The bias conditions for curve 1 are 
Ves = —10v and Jz = 10 ma; for curve 2, Vcg = —6v and Iz = 6 ma. 
The bias conditions for the measured response are the same as for curve 
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Fig. 4— Amplitude-frequency response with inductance added in base of 
transistor. 


2 in all except the last two stages which are discussed later. As suggested 
by Fig. 4, reducing the bias power results in some sacrifice in bandwidth. 
The final selection of bias conditions is a compromise between a flat 
response over the desired bandwidth and low power consumption for 
the amplifier. For the computation of curve 2 the circuit parameters 
are: 


fa = 400 me Cr = 4.5upf R;, = 59 ohms 
r, = 50 ohms In = 0.05ph 
re = 4.3 ohms L = 0.02uh 


The measured midband insertion gain between 75-ohm impedances is 
41 db for eight interstages. The de power consumption is 0.7 watt. 


lil. INPUT AND OUTPUT CIRCUITS 


A complete schematic diagram of the amplifier with a matched input 
circuit and two output circuits, is shown in Fig. 5. The input circuit 
shown is similar to the interstages in that a resistor is added in series 
with the emitter to provide the proper input resistance of 75 ohms, and 
a capacitor is shunted across the input to peak the circuit at higher 
frequencies. The return loss of the amplifier input as a function of fre- 
quency is shown in Jig. 6. 
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Fig. 5 — IF amplifier schematic. 
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Fig. 6 — Return loss vs frequency for input circuit of IF amplifier of Fig. 5 
measured at 75 ohms. 


The input conditions for operation with a low noise figure are discussed 
in Section VIII. 

Two output circuits are shown in Fig. 5. One output circuit consists of 
a broadband transformer and a shunt peaking circuit. The other has a 
double-tuned transformer and is shown as the alternate output circuit 
in Fig. 5. The return losses of the two output circuits, measured at 75 
ohms, are given in Fig. 7. The double-tuned transformer adds a measured 
variation of about one-tenth of a db to the frequency response. All data 
discussed in this paper were taken with the output circuit having the 
broadband transformer. 

The AGC detectors shown in Fig. 5 have been designed so as to have 
no effect upon the frequency response. AGC problems are discussed in 
Section IX. 


IV. POWER OUTPUT AND FREQUENCY RESPONSE 


The power output of the amplifier depends upon the bias conditions 
in the output transistors, which in turn depend upon the maximum 
ambient temperature expected. The last two transistors in Fig. 5 have 
increased bias currents to allow for increased power output. The bias is 
arranged so that the rated junction temperature is not exceeded at the 
high end of the ambient temperature range, 130°F’. 

The compression characteristic of this amplifier for a single frequency 
sine wave input at 70 mc is given in Fig. 8. In addition to this information 
it is desirable to know the frequency response as a function of output 
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Fig. 7 — Return loss vs frequency for output circuits of IF amplifier of Fig. 5 
measured at 75 ohms. 


level. This is shown in Fig. 9. Fig. 9(a) shows the response for an output 
power of +6 dbm and a negligible compression. Fig. 9(b) shows responses 
for several other levels where it is to be noted that while the details of 
the responses have been preserved, the level differences between the 
curves have not. 


V. TEMPERATURE-FREQUENCY RESPONSE 


The frequency response as a function of temperature has been ob- 
tained for the amplifier that does not have added inductance in the base 
circuit. The data are presented in Fig. 10, plotted relative to the gain at 
50 mc. The 50-me gain decreases at high temperatures and increases at 
low temperatures. The total change in 50-mc gainin therange —40°F 
to +130°F is under 2 db. 


VI. TRANSFORMERS 


The interstage transformers used in this amplifier have been described 
elsewhere’ and only a brief description is given here. These transformers 
are of a transmission line type and consist of a short bifilar winding on a 
small nickel-zine ferrite toroidal core. The winding is formed by twisting 
a pair of wires and then winding the pair on the core. Fig. 11 is a sche- 
matic of the transformer and a plot of the measured frequency response 
into resistive terminations. The turns ratio is 2:1. 

This type of transformer may be regarded as an ideal transformer in 
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Fig. 8 — Compression characteristic for amplifier of Fig. 5. 


series with a transmission line whose electrical length equals that of the 
bifilar winding. The analysis in this paper has assumed only the ideal 
transformer, which is reasonable as long as the winding is short compared 
to a wavelength at the highest frequency of interest. 

The permeability and Q of this ferrite has been optimized in this fre- 
quency range to give minimum loss. With this ferrite and a winding 1.6 
inches long, the transformer has the performance shown in Fig. 11. 


VII. GAIN PER STAGE 


All of the current gain is obtained from the interstage transformer so 
the maximum current gain per stage using a transformer with an im- 
pedance ratio of 4:1 would be 6 db. There are three sources of current 
loss in each interstage which reduce this gain. As seen in Fig. 1(d), they 
are: 

(i) The transformer inductance Ly and resistance Rr which shunts 
the emitter-base circuit of Q2. 

(11) The emitter bias resistor R: which shunts the input circuit of Qo . 
This loss can be reduced at the expense of greater bias power. 

(iii) The current gain of Q@. = a < 1. 

The total loss from these sources is 0.8 db so that the actual measured 
gain per stage is 5.2 db. 
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Fig. 9 — (a) Frequency response of amplifier and base line at output power 
level of +6 dbm and insertion gain of 41 db, (b) frequency response of IF amplifier 
as a function of the output level in dbm. 


VIII. NOISE FIGURES 


The amplifier of Fig. 5 has its input impedance adjusted to 75 ohms 
as indicated in Fig. 6. With this input circuit it is suitable for use as a 
main II amplifier. or use as a preamplifier, however, the input circuit 
can be modified to improve the noise figure. Such a circuit is shown in 
Fig. 12. Note that the turns ratio of the interstage transformer is 8:3. 
The measurements reported below were made on the amplifier of Fig. 5 
with the input circuit of Fig. 12 and using a diode noise generator. The 
bias on the first stage was J; = 3 ma, Vcg = —7V. 

Some amplifier noise figure measurements with W.E. 2N1195 transis- 
tors in the input circuit are shown in Jig. 13. The curve is the noise 
figure for the best of 5 transistors, and the amplifier noise figure is plotted 
versus generator resistance. Also shown are the amplifier noise figures 
for the 5 transistors with a generator resistance of 70 ohms. These meas- 
urements were made in a 0.7-mc band centered at 71 mc. At 50 me the 
noise figure is typically 0.5 db lower and at 100 me, 0.75 db higher. 

It should be emphasized that these are amplifier noise figures and not 
transistor noise figures. In such a broadband amplifier, care must be 
taken to minimize the noise contribution of the second transistor. This 
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Fig. 10 — Frequency response vs temperature of IF amplifier. 


accounts for the higher turns ratio of the interstage transformer. Use of 
the transformer described in Section VI results in an increase of 0.5 db 
in noise figure. 

The noise contribution of the second transistor depends upon the 
available gain of the preceding circuit,’ and the available gain depends 
in part upon the generator impedance and the transformer turns ratio. 
It is to be expected, then, that for narrower bandwidths the noise contri- 
bution of the second transistor can be further reduced, since higher turns 
ratios can be used with the corresponding increase in available gain of 
the first stage. 

Yor the broadband common-base amplifier the best noise figures are 
obtained with transistors having low optimum source resistances. This 
comes about because the available gain increases as the optimum source 
resistance decreases, and the greater the available gain, the lower the 
noise contribution of the second transistor. It follows that for best results 
a low-noise transistor having a low optimum source resistance should be 
used. 

The effect of the input circuit of Fig. 12 on the frequency response of 
the amplifier is quite small. With this input circuit the response was 
measured using a 75-ohm source impedance and the response equivalent 
to that of Fig. 9 was obtained by slight readjustment of the base lead 
inductance in the second stage. 


IX. AUTOMATIC GAIN CONTROL CIRCUIT 


Wideband transistor IF amplifiers present special problems with re- 
spect to AGC in that it is not always feasible to merely adjust the bias 
to obtain gain control. The method presented here consists of an elec- 
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Fig. 11 — Transmission-line type transformer with impedance transformation 
of 4 to 1. (a) Transmission line form, (b) low-frequency equivalent, (c) amplitude 
vs frequency of transformer measured between 75 ohms and 18.75 ohms. 
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Fig. 12 — Input circuit for low noise figure. 
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Fig. 13 — Noise figure measurements. 
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Fig. 14— AGC attenuator. 





+3.5 


WIDEBAND IF AMPLIFIER 51 











AMPLITUDE IN DECIBELS RELATIVE 
TO 75MC VALUE 


50 55 60 65 70 75 80 85 90 95 100 
FREQUENCY IN MEGACYCLES PER SECOND 


Fig. 15 — Frequency response of variable IF attenuator of Fig. 14. Output 
level held constant at —37 dbm. 


tronically variable attenuator which can be inserted in the amplifier at 
any suitable point. The attenuator is described schematically in Fig. 14 
and is basically a T-pad formed by three diodes. The three diodes are 
supplied with current at point A through the 750-ohm resistance. Tran- 
sistor Qs controls the current through diode 3 which in turn controls the 
current through diodes 1 and 2. The resistance of the diodes is a func- 
tion of the current through them so the attenuation of the T-pad can be 
controlled by adjusting the current in diode 3. This is done by applying 
the AGC detector output to the emitter of transistor Q3 . 

If Q; is cut off, the current at point A splits between diodes 1 and 2, 
making them low resistances; diode 3 becomes a high resistance. This is 
the minimum-loss condition. When Q; conducts heavily, all of the current 





Fig. 16 — Wideband IF amplifier. 
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at point A goes to diode 3, lowering its resistance and cutting off diodes 
1 and 2. This is the maximum-loss condition. The reference diodes (TI- 
650) provide a 4-volt bias that assures that diodes 1 and 2 can be cut 
off and also provides collector voltage for Q3. 

The frequency response of the attenuator is given in Fig. 15. These 
curves show the variation in response relative to the 75-mc value as the 
power input to transistor Q, (Fig. 14) is varied from —1 dbm to —31 
dbm. For inputs greater than +3 dbm the response changes shape rather 
than slope. 

The minimum measured insertion loss of the circuit of Fig. 14 is 6 db. 
The power consumption of the attenuator is 0.35 watt. 


X. MECHANICAL ARRANGEMENT 


The photograph of Fig. 16 shows the mechanical arrangement of one 
section of the amplifier which has 41 db of gain. The transistors are 
mounted inside the chassis to minimize stray inductances and to prevent 
accidental grounding of the collector, which is connected to the transistor 
header. 


XI. DISCUSSION 


It has been demonstrated that high quality wideband IF amplifiers 
can be made using transistors. Such amplifiers are capable of operation 
over large temperature ranges and small input power. For instance, an 
amplifier with a flat bandwidth of 50 mc, 90 db of gain, +12 dbm of 
output power, 34 db AGC range, a noise figure of 6 db, and capable of 
operation over a temperature range of —40°F to +130°F is feasible. 
Such an amplifier would require about 2 watts of bias power and would 
contain 23 transistors and 6 diodes. 

It is possible, of course, to trade one attribute for another. For in- 
stance, a flat band extending to nearly 200 mc can be achieved at the 
expense of increased bias power. If a narrower bandwidth is desired, 
then the gain per stage can be increased. This can be done, for example, 
by using two transformers per stage or using another type of transformer. 


APPENDIX 


The input impedance of the equivalent circuit of Fig. 1(c) is approxi- 
mately 


Zw ret Al — a) (1) 


where 
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Z(1 ee a) > Zr 


Le = LL + Zp . 
Now let 
ao 
= ——— and wl 
am ay ™ 
oa f/f 
LA = rp + aX. 


Xy = oly = YHalsn = V¥Xba 
After these substitutions (1) becomes 


(te te) +) = (re) + YP Xt0) + jr’ + Xba) 
oy eo ee See ge ms ( 


Thus, by adding inductance J in the base circuit the input impedance 
is changed by the addition of the real term 


[—7’Xse/(1 + y’)] and the reactive term j[Xsay’/(1 + y’)). 


Zin 2) 


When the input impedance of the transistor becomes part of the peak- 
ing circuit of Fig. 1(b), it becomes part of the resonant circuit formed by 
Cr, Ri, L, and Zry . The total inductance of this circuit will increase 
and the total resistance will decrease with increasing frequency due to 
the variation of Z;y , thus changing the shape of the resonance curve. 
The greatest effects will be observed at the higher frequencies. 

The peaking circuit to be analyzed is shown in Fig. 1(b). The total 
interstage capacitance is C7 , stray inductance L, and inductance added 
at base L,. R; is used to adjust circuit Q. The output current of tran- 
sistor Q, is Jo . The transformer has been assumed ideal. 

Io is the emitter current of Q. . Solving this circuit we get 


Io _ Ba + Rw + J Kty — Xe) 


: 3 
i 4K (3) 
and substituting from (1) and (2) 

ry + Xo, .[ (re + Xs) _ 
| cata “ie © ee a, Xe| 
Is =x (4) 


= Rr—-—pt juny + Xz — Xe) 
Xs 
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Inf’ _ (Rr — »)” + (uy + Xt — Xe)? (5) 
Ts Xo ' 
where 
Rr=Ritre +r 
- te + yuls 
Z dag? 
Y= ti fa . 
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Monitoring the Percussive Welding Process 
for Attaching Wires to Terminals 


By J. C. COYNE 
(Manuscript received September 5, 1962) 


At present, no satisfactory nondestructive test is known which can be 
practically applied in the shop to detect a weak percussive-welded connection. 
Also, the breaking strengths within a population of percussive welds are 
distributed in such a way that large sample sizes are needed to determine, 
with reasonable confidence, if the population contains an excessive number 
of weak welds. 

Momitoring the duration of the welding arc and the approach speed of the 
wire, for each weld as it 1s made, has been found to provide an effective 
control of the process. An analysis of the process yields necessary conditions 
which the are duration and wire speed must satisfy. Test data are presented 
which confirm the necessity of these conditions, and show that when they 
are satisfied, percussive-welded connections can meet quality objectives. 


I. INTRODUCTION 


The low-voltage percussive welding process as applied to wired con- 
nections has been described and its merits pointed out in previous 
papers.'-§ Therefore, only a brief review is included here. 

Percussive welding is a form of capacitor discharge arc welding. The 
parts being welded, in this case a wire and a terminal, are themselves 
the only necessary electrodes (see Fig. 1). Also, the only advance prepa- 
ration of the electrodes needed is to cut a tip on the wire end. The weld 
is made by propelling the wire across a 50-volt gap toward the terminal. 
When the two are nearly touching, an arc is established between the wire 
tip and the terminal. Initially, the power surge in the arc develops 
sufficient heat to melt back the wire tip faster than the electrodes are 
closing. But, with time, the power decays and melting slows down until 
eventually the electrodes impact, extinguishing the are and completing 
the weld. A welding voltage of 50 volts is commonly referred to as “low 
voltage” as opposed to “high-voltage percussive welding”’ used to attach 
the fixed contacts to the wires of the wire spring relay. 
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Fig. 1— Welding circuit. 


The welded joint made by this process is usually as strong as the wire 
itself. Nevertheless, it appears to be characteristic of the process to 
produce also, if adequate controls are not present, a small scattering of 
weaker welds, some ranging nearly to zero strength. Naturally, such a 
distribution of weld strengths is unacceptable from a connection relia- 
bility standpoint. Moreover, the usual statistical sampling techniques 
for assuring satisfactory quality in production have limited value here, 
as shall be seen later. For this reason, other more effective measures were 
sought. 

Extensive investigations were first conducted to find a nondestructive 
screening test which might be applied to each completed weld. However, 
at the present time no satisfactory test has been devised which can be 
practically applied in the shop. More recently, attention has centered 
on control of the process by monitoring certain important welding 
parameters as each weld is made. The analysis and tests pertaining to 
monitoring are the subjects of this paper. 

In the preliminary investigations, the time histories of wire velocity, 
current, and voltage were examined for aberrations which always coin- 
cided with the occurrence of a bad weld. Although this study revealed 
no easily measurable, consistent predictor of bad welds, it was observed 
that welds were consistently good if the time duration of the arc was 
within a certain range. Conversely, the percentage of defective welds in- 
creased for either longer or shorter times. Also, this optimum range 
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shifted, compressed, or expanded depending on the welding circuit con- 
stants, the wire velocity, and the wire tip shape. 

Further testing and analysis were done to explore the effectiveness as 
well as the practical possibilities of monitoring arc duration. Two factors 
seemed to be in favor of this. First, the are duration is a sensitive indi- 
cator of most possible changes in welding conditions. For instance, varia- 
tions in the electrical characteristics of the welding circuit, wire tip 
geometry, or approach speed of the wire will all affect the duration of 
the arc. Second, it is a relatively simple measurement to make. 

The areas of fractures for welds occurring at either very short or very 
long are durations were observed to have certain distinguishing char- 
acteristics. In the case of very short arc durations, the welded area is 
usually noncircular and has less area than the cross-sectional area of the 
wire, as shown in Fig. 2(a). This results from the fact that the chisel- 
shaped tip on the wire does not become completely melted back during 
the short time duration of the arc. For the same reason, the crater burned 
into the terminal is usually undersized. The loss of weld strength in this 
case depends on the reduction of weld area. 

In the case of very long are durations, the smooth break along the 
interface of the two metals indicates an absence of fusion and suggests 
that the molten metal film over the electrodes had prematurely solidified. 
Fig. 2(c) shows an extreme example of this. This type of weld is invari- 
ably quite weak. As the arc duration is shortened, a corresponding in- 
crease in the average area of fusion in the weld is observed, but with 
considerable variability. A comparable improvement in the average 
weld strength also occurs. 

If the arc duration is further shortened until it is in the optimum 





(a) (b) (c) 


Fig. 2— Appearance of weld breaks: (a) Excessively short arc durations do not 
provide sufficient melting energy and result in a weak weld due to the reduced 
area of fusion. (b) Arc durations within the optimum range result in a strong weld. 
(c) Excessively long are durations allow excessive heat conduction losses and re- 
sult in a weak weld due to premature solidification. 
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range (assuming that certain other conditions are also satisfied), then a 
full-area weld of wire to terminal will result as shown in Fig. 2(b), and 
the weld strength will equal or exceed the wire strength. The other con- 
ditions referred to above include: (z) No serious geometrical disparities 
in the surfaces of the arcing electrodes which would prevent a uniform 
full-area mating when contact is made. (i7) Adequate melting to remove 
the entire wire tip. (277) No spreading or other disturbance of the arc 
due to surrounding fixture metal at electrode potential being too close. 
(iv) Development of sufficient force between electrodes at contact to 
effect good mating and fusion. 

The foregoing observations suggest that premature solidification is a 
major cause of partial-area fusion. Excessive heat losses will occur at 
long arc durations, since the arc power decays at an approximately 
exponential rate, while the rate of heat conduction losses decreases in- 
versely as the square root of time. Consequently, heat conduction losses 
grow with respect to arc power until eventually heat is being conducted 
away faster than it is being generated. Although radiation and convec- 
tion losses are also present, these are small by comparison. 

A second major cause for partial fusion in the weld zone is that the 
contours of the arcing surfaces can interfere with full area mating. This 
is particularly true at very short arc durations for reasons previously 
mentioned. However, to varying degrees, there will always be some 
differences in contour, because the two electrode shapes are vastly differ- 
ent at the start of the process. A proper choice of wire tip shape will 
minimize the disparity in the contours of the arcing surfaces. Any small 
differences that might still be present will not be harmful if the elec- 
trodes impact with sufficient force. 


Il. ANALYSIS OF WELDING PROCESS 


2.1 General 


Having seen in a descriptive way why either extreme of arc duration 
is undesirable, the welding process will be examined analytically in 
order to derive, for the general case, maximum and minimum monitoring 
limits of are duration. 

The major independent variables of the process to be used in this 
derivation are: (7) the wire and terminal metals; (27) the welding 
voltage, V., which will be taken as 50 volts; (cz) the welding circuit 
effective resistance, R, and inductance, LZ, which will be treated as 
‘Jumped”’ constants; (zv) the power supply capacitance, C, also a 
lumped constant; (v) the wire velocity, U, which will be taken as the 
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average wire velocity during the time interval of the are; (vz) the shape 
of the wire tip, especially its length, /, and volume, ». 

The more important dependent variables to be used in the derivation 
are: (¢) the welding current, J; (27) the arc voltage, V; (2i7) power ex- 
pended in the arc, P; (i) heat conduction losses, Q; (v) the burnback 
of the wire tip, b; and (vz) the are duration, f, . 


2.2 Some Relationships and Approximations Useful in the Derivation 


The first step in the derivation will be to establish useful relationships 
for are voltage, current, power and heat losses. Wherever possible, 
approximations will be made to simplify the formulation. 

The are voltage (V) will be approximated by a straight line function, 
as shown in Figure 3: 


(V1 7 V2)t (1) 


where V; and V2 are the voltages across the arc just after arc initiation 
and just prior to are extinction respectively. For copper wire and nickel- 
silver terminal metals, V2 has been found to be approximately a 10- 
volt constant, independent of welding voltage and current, wire velocity, 
tip geometry, or wire gauge. Investigation of a few other metals also 
showed no appreciable variation from this 10-volt constant. 

V, showed a slight dependence on the parameters previously men- 
tioned, but in each case it was noted that an increase in V; corresponded 
to conditions producing a faster melting rate along the wire relative to 
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Fig. 3— Typical traces of current and arc voltage. 
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the wire speed. For instance, V; would be increased if either the wire 
speed decreased or a longer and thinner tip were cut on the wire. In 
other words, it appears that the arc voltage is mainly a function of the 
instantaneous separation of the welding electrodes. 

Although V; varied in the manner described, the magnitudes of the 
variations were small — in the order of 4 volts even when the peak weld- 
ing current ranged from 300 to 900 amperes and the wire velocity ranged 
from 50 to 100 inches per second. A reasonable value for V; , which will 
be used in the remainder of the analysis, was found to be 20 volts. 

The differential equation describing the current flow for the welding 
circuit of Fig. 1 is: 

dl 
To 


with initial conditions at zero current and Vo volts across C’. Substitut- 
ing (1) into (2), the welding current may be determined. The general 
solution for welding current consists of two parts: the standard “homoge- 
neous solution” for a capacitor discharge through an RL series circuit, 
and an additive “particular solution” arising from the time dependent 
are voltage function. The damping ratio of the welding circuit will be 
defined in the usual manner as the dimensionless ratio: 


p= 24/5 (3) 


and confined to a range of one to infinity. For D equal to one (critical 
damping), the welding current becomes: 


us A Ve= Vile CV = Ve) (27 + 1) —2r Vi — Vo 
eg ee ee 


where + = ¢/RC and ta = ta/RC. 
For D approaching infinity (zero inductance) the welding current 
becomes: 


dl I dV 
+R +547, =0 (2) 


_jJVo— Vi) (Vi — Ve)| = , Vi Ve 
= ets eo eee aaa "a ~ 
For intermediate values of D, the current magnitude will be between 


(4) and (5). 

When 7 equals 7, , either equation gives the current at arc extinction 
(I,). If 7a is sufficiently large, say greater than 4.0, then J, from both 
(4) and (5) becomes 
Vi — Vo 


Rrra (6) 


since in either case the exponential term becomes small by comparison. 


Iq 


IR 
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If 7, is not sufficiently large to use this approximation, then the con- 
tribution to J, made by the exponential term should be included. Ac- 
count is taken of this contribution by rewriting (6) with N substituted 
for Vi — Ve: 


_N _wNe 
 Rrq te 


where N is a number larger than Vi; — V2 and dependent on 7,. N may 
be calculated by equating (7) to (4) or (5), with r = 7,4. Some repre- 
sentative values are shown in Table I, calculated for V, = 50, Vi = 20 
and V2 = 10 volts. 

The instantaneous power, P, expended in the arc is the product of V, 
given by (1) and J, given by (4) or (5). At are extinction, using (7) for J, 
and with ¢ = i., 


Le (7) 








NCV2 


Pos 
ta 


(8) 

The electrical energy, Hy , dissipated as heat in the arc, as a function 
of time, is obtained from a time integral of the power. To simplify the 
formulation without any great loss of accuracy, it will be assumed that 
the are voltage remains constant (equal to V1) and therefore the second 
terms of both (4) and (5) vanish. Then for D = 1 


(Vo — Vi)te" 
R 


pay, (9) 


and for D — « 


(Vv. — Vide“ 


Using this approximation, the errors in arc voltage and current are of 
opposite sign, and therefore these partially cancel in the product. The 
remaining error is minimum at arc initiation and maximum at arc ex- 
tinction. Table II has been prepared to show the relative magnitude of 


Pe VV; (10) 


TasLe I — Vauues or N In (7) 
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TasLe II —Arc Power VERSUS r 


10 
(5) times (1) leppioumatiogs 


0 600/R 600/R 
1 218/R 221/R 
2 92/R 81/R 
3 = Ta 47/R 30/R 


the error made by using the approximation. The calculations have been 
made for r, = 3.0 and for the same values of V, , V1, and V2 used previ- 
ously. Since the contribution to the total arc energy is many times 
greater at arc initiation than at are extinction, it is evident that the 
percentage error is further reduced by the integration. 

Therefore, for D = 1 


Eg & CVV. — Vi) CL — €* — 2re7) (11) 
and for D—> 
w= CVi(V. — Vi) (1 — €7). (12) 


The rate of heat conduction losses (dQ/dt) at any time after arc ini- 
tiation can be taken as 


dt ~ Vrat Vi oe 


where A is heat conducting area; e.g., wire area, 
k is thermal conductivity, 
a is thermal diffusivity, and 
T is molten metal temperature. 

Substituting, for instance, the appropriate constants for 24-gauge 
copper wire, K becomes approximately 7.0 watt (sec)?. This number is 
based on using an estimated average value for 7’. There is evidence that 
some very small amount of copper vaporizes during arcing.®: ° Conse- 
quently, the molten copper must range between its melting temperature 
and its boiling temperature. 

The total heat conduction loss between arc initiation and time ¢ is 


Q = 2K Vt. (14) 


2.3 Maximum Limit on Are Duration 


A maximum limit on arc duration, tamaxy, may now be specified in 
order to prevent premature solidification. It is required that the power 
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dissipated in the are just prior to are extinction always exceed heat 
conduction losses. Assuming half of the are energy is expended at each 
electrode and that radiation and convection losses are negligible, 

P. — dQ 

a ee 15 

a a (15) 
Implicitly, it is postulated here that if the above inequality is satisfied 
at arc extinction, then it must necessarily be satisfied for all earlier 
times. 

Substituting (8) and (13) into (15) and solving for ta 


NCV’ 
ta < ( oR ‘ = latmax) . (16) 





Equation (16) is our first monitoring inequality and it is plotted in Fig. 
4 for different RC products, and for the two extreme damping condi- 
tions, versus capacitance. 

This formulation of tagnax) includes an inherent safety margin with 
regard to monitoring the process since the time to accomplish solidifi- 
cation is not taken into account. 


2.4 Minimum Limits on Arc Duration 


Burnback (0) is the length of wire melted back during arcing. With 
negligible error, it can be assumed that the electrodes are touching at arc 
initiation. (A potential of 50 volts will not normally initiate an are 
across the air gap until it has closed to approximately 0.00005 inch.)4 
Consequently, the distance advanced by the wire during arcing must be 
equal to the average wire velocity, U, multiplied by the arc duration. 
This product will be very nearly equal to the total melted wire length, : 
differing only by the thin layer of molten copper on the wire and by the 
crater depression in the terminal. These are both relatively small and 
of opposite sign. Therefore, 


b= Ut. (17) 
In order to melt the entire wire tip length, J, 
b > 1+ safety margin. (18) 


If the wire velocity is measured by the time, t, , for the wire to tra- 
verse contacts separated by 1, then (18) reduces to the monitoring in- 
equality: 


la > ly. (19) 
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Fig. 4— Computed curves for ta(max): 24-gauge copper wire, nickel-silver 
terminal, welding voltage 50v, Vi = 20v, Ve = 10v, K = 7 watt(see. )i. 


Preliminary tests indicated the need for another minimum limit on 
arc duration, ta(min) , independent of t, , since a scattering of weak welds 
sometimes occurred at slightly longer arc durations than predicted by 
(18). In general, fain) must be at least long enough to provide for 
sufficient melting energy; however, other reasons may make a longer 
time necessary. This is discussed further in Section V. On this basis, 


la > ta(min) (20) 
where facmin) 18 determined empirically. For the tests reported in Sec- 
tion IV, teminy was found to be 200 microseconds. 

2.5 Necessary Power Supply Capacitance 


The melting energy (HZ) necessary to remove the entire wire tip, 
assuming the molten metal is at its melting temperature, is 


PERCUSSIVE WELD MONITORING 65 


Ey = (Tyce + A)pv. (21) 


where 7’, is the metal melting temperature, c is the specific heat, \ is the 
latent heat of fusion, p is the metal density and v is the wire tip volume. 

The necessary power supply capacitance, C’, can be calculated from an 
energy balance 


1Eyn —-Q= Ex. (22) 


Substituting in (12), (14), and (21) and evaluating for average damp- 
ing conditions, for copper wire and for f4 = taimaxy) = 3RC, 


C & 600av (23) 


where a is a Safety factor at tamax). In general, a must include allowances 
for the reduction in input energy resulting from shorter arcing times as 
well as the greater demand for energy resulting from the variational 
increases in tip volume. The safety factor at tagniny, (6), will be, of course, 
smaller than a and may be evaluated from the following approximate 
relationship 


rae 1—- exp (— oni? CO) 


PT — exp (hi mn/RO) © 


(24) 
An independent measurement of the melted wire volume versus capaci- 
tance agreed closely with the volume predicted by (23). 


2.6 Monitoring Inequalities 


To recapitulate, the inequalities in (16), (19), and (20) restrict the 
allowable range of are durations for percussive welding 


. < ta < latex) (25) 


ta (min) 


where faimin) aNd taimax) are fixed limits evaluated either from the con- 
stants of the process or by empirical methods, and ¢, and ¢, are measure- 
ments made for each weld. The measurements of f, and #, can be made 
in a number of ways. Fig. 5 shows one method, in which contacts sepa- 
rated by the minimum burnback distance provide step voltages which 
are indicative of ¢, , and a coil around one welding lead provides induced 
voltages at the beginning and the end of the are which are indicative of 
t. . These voltage indications are fed to a timer and alarm — any suitable 
apparatus capable of determining whether the monitoring inequalities 
in (25) are satisfied. Various other measurement methods might be em- 
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Fig. 5 — Wiring schematic of monitoring circuit. 
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ployed, such as a light beam and photocells in the case of ¢, and the are 
voltage change in the case of é,. 


IiI. OTHER CONSIDERATIONS 


The described monitoring method naturally does not preclude the 
possibility of a defective weld. Other causes still exist. 


3.1 Bad Wire Tip 


Inconsistently or improperly formed tips are a potential hazard. 
Usually, monitoring will detect this condition since arc duration is sen- 
sitive to variations in either tip length or volume. Nevertheless, certain 
tip variations — for instance, variations in symmetry — are undetected 
by the measurement of are duration. 
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3.2 Arc Dispersion 


The wire as it approaches the terminal must ‘‘see’’ the terminal to the 
exclusion of any other metal at terminal potential. Otherwise the arc 
will spread and be weakened, resulting in poor welds despite monitoring. 


3.3 Mechanical Effects at Closure 


The most critical time in the weld process is, of course, when the two 
weld metals come together and fuse. The possible influences on weld 
strength at this time are manifold; the following are representative. At 
closure, the arcing surfaces should match reasonably well in order to 
achieve a uniform fusion zone. Also, any large disparity in contour will 
cause a time lag between initial electrode contact (arc extinguishes) and 
final full-area contact. The time lag is increased if the wire is appreciably 
decelerated. During such a time lag, premature solidification might be a 
problem in those areas that mate last. The shape of the original wire tip, 
especially its symmetry, plays an important role in determining how 
well the electrode surfaces will mate, since arcing generally is concen- 
trated between the nearest points on the two electrodes. 

With regard to impact forces at the weld interface, it has been found 
that up to a point these forces improve both connection strength and 
fatigue life. This might be attributed to various effects, some of which 
could be, (7) overcoming small contour differences on the arcing sur- 
faces, (77) heading, or increasing the area near the weld zone, (22) strain 
hardening the copper above the weld, (2v) promoting better mixture in 
the fusion zone. 


IV. DESCRIPTION OF TESTS 


4.1 General 


Welds were made for test under controlled laboratory conditions and 
monitored in accordance with the analysis. In order to test only the 
effects of the parameters in the monitoring inequalities, all other welding 
conditions were optimized to the best of our knowledge. 

Wire tips were cut with a hand-held pair of diagonal wire cutters. 
Despite its seeming crudeness, this tip has been found to be sufficiently 
reproducible and has given good results in the past. The tip length on 
the wire averaged 0.019 inch with approximately half the volume of an 
equal length of wire. The wire was 24-gauge tinned copper. 

The terminal was held flush with or slightly projecting from the ter- 
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minal clamp in order to minimize arc dispersion. The terminal was 
0.025 inch x 0.062 inch nickel-silver. 

The welding power supply contained 2300 microfarads capacitance. 
This value produces an adequate energy margin. No delay line or other 
type of network was used to supplement the capacitor bank. . 

The wire was accelerated with an adjustable tension spring and was 
permitted to fully impact into the terminal so that the weld interface 
experienced the complete impact forces. Weld strengths were measured 
by the “combined test”? whereby the wire is bent over parallel to the 
terminal and pulled until either the wire or the weld fails. This test has 
proven to be more sensitive and more severe than a straight tension test. 
The break was classified as a “wire break” if the failure occurred in a 
portion of the wire completely above the terminal surface. Actually, 
these “wire breaks” occurred randomly along a span of wire of about 
three inches, so that usually there was very little difficulty in distin- 
guishing a ‘“‘wire break” from a “weld break.” Furthermore, if a break 
occurred close to the weld-affected zone, and its classification was doubt- 
ful, it was usually taken to be a ‘‘weld break” since the probability is 
small that a “‘wire break” would occur randomly this close to the weld. 


4.2 Phase I 


The testing was divided into two phases. In Phase I, the burnback 
inequality at all times was satisfied, whereas the arc duration was varied 
in order to demonstrate that (16) and (20) are necessary conditions for 
optimizing the welding process. The arc duration was varied by means 
of changing the wire velocity. At least 600 welds were made and tested 
within each 50-microsecond interval between 150 and 500 microseconds, 
making a total of 4730 welds. The testing was not extended below 150 
microseconds because of practical limitations of high wire velocity. 

The following steps, based on the analysis, show the calculations made 
to determine the monitoring limits applicable to this first phase of test- 
ing. 

(a) The welding circuit constants were determined by substituting 
measurements from an oscilloscope current trace into the theoretical 
current expression. Accordingly, the effective total series circuit re- 
sistance was calculated to be 0.040 ohm, the inductance 0.75 micro- 
henry and the power supply capacitance 2300 microfarads. The capaci- 
tance figure agrees closely with an independent measured value. From 
(8) the damping ratio, D, becomes 1.1. The RC product becomes 92 
microseconds. 

(b) From Fig. 4, using the critical damping curves, taanax) Is 320 mi- 
croseconds. 
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(c) In accordance with previous tests at the same welding conditions, 
tatminy Was taken at 200 microseconds. 

(d) The average wire tip volume was measured to be approximately 
3 X 10-§ cubic inches. Using (23), the safety factor (a) at 320 micro- 
seconds is calculated to be 1.25 for a capacitance of 2300 microfarads. 
The safety factor is reduced to about 1.15 at 200 microseconds using (24). 


4.3 Phase II 


In Phase II, the burnback was varied in order to demonstrate the 
necessity of adequate burnback (18). The burnback was varied from 
0.012 to 0.022 inch to span the average tip length of 0.019 inch. At 
least 200 welds were made and tested within each 0.002-inch increment, 
making a total of 1126 welds. The test was planned to be within the 200 
to 350-microsecond interval of arc duration at a constant wire velocity of 
63 inches per second. The only welding parameter independently varied 
was circuit resistance. This method was chosen not only to facilitate a 
sensitive control of burnback but also to avoid the possibility of exceed- 
ing tagnax) , Since it was desired in this phase to test only the effect of re- 
ducing burnback. 

The computations of monitoring limits for Phase II follow. Since 
resistance was varied in this test, no one value of tagmax) applies through- 
out. But to indicate a typical case, at the longest range of arc durations 
in this test (820-350 microseconds), the incremental resistance was ap- 
proximately 0.060 ohm, making the total resistance 0.1 ohm. The RC 
product then is 230 microseconds and fa(max) 1s calculated from (16) and 
Table I to be approximately 600 microseconds (N = 15; ta = 2.5). 

From (24), the safety factor for energy at 320 microseconds is approx- 
imately 1.0, which means that this is the minimum allowable are dura- 
tion on a “sufficient energy basis.”” The burnback at 320 microseconds 
and 63 inches per second is 0.020 inch, in good agreement with the aver- 
age tip length (0.019 inch) and giving us a check on the calculations. 
Consequently, in this test we can expect weld quality to degrade for 
are durations less than 320 microseconds. 


V. DISCUSSION OF TEST RESULTS 


5.1 General 


The test results of Phase I, shown in Fig. 6, clearly demonstrate that, 
for the conditions of this test, the welding process is optimized when the 
are duration is between 200 and 300 microseconds. Within this range, 
99.2 per cent of the samples failed in the wire and 99.85 per cent of the 
samples exceeded 90 per cent of the wire strength. 
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Fig. 6 — Results of test No. 1, weld strength vs arc duration: 24-gauge copper 
wire, nickel-silver terminal, capacitance 2300 uf, resistance 0.04 ohm, diagonal- 
cut wire tip, variable wire velocity. 


The estimated limits of the optimum range, 300 and 200 microseconds, 
agree closely with both the results of previous tests and the analysis. At 
arc durations less than 200 microseconds there is a relatively fast drop- 
off in weld quality. At arc durations longer than 300 microseconds there 
is also a dropoff in quality, although not as rapid as in the first case. As 
the arc duration deviates further from optimum toward longer times, the 
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Fig. 7 — Results of test No. 2, weld strength vs burnback. Conditions same as 
test No. 1, except that wire velocity is constant 63 in/sec and resistance is varied. 


points lie along the 63 inches per second velocity line. Table III sum- 
marizes the test results for the regions of Fig. 8 and indicates the moni- 
toring inequalities which were violated in each case. The “‘order of merit”’ 
of weld quality is based on the test results in Figs. 6 and 7. Table III, 
in conjunction with Fig. 8, demonstrates the necessity of each of the 
monitoring inequalities. 

In Fig. 8 it is seen that the constant burnback line of 0.020 inch 
forms the boundary between regions D and E. The fact that this bound- 
ary, only 0.001 inch more than the average tip length, effectively 
separates the good from the poor welds implies that burnback is the 
controlling factor in weld quality in Test Phase II. An independent 
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percentage of weaker welds begins to level off. Although this was unex- 
pected, it may be explained plausibly on the basis of varying solidifica- 
tion time. 

The time interval that can safely elapse after tajmax) depends, of course, 
on the heat contained in an unknown thickness of molten film over the 
electrodes. This thickness is probably quite variable. At normal welding 
velocities, it is quite thin; however, there is evidence that it generally 
increases on the wire for longer arc durations (slower wire speed).® There- 
fore, the time to solidify the molten film will also be variable and will, 
in general, become longer at longer arc durations. 

The results of the second phase of testing, shown in Fig. 7, indicate 
a rapid loss of weld quality as the burnback is decreased. In agreement 
with the preliminary calculations, the weld quality was relatively good 
for arc durations greater than 320 microseconds or burnback in excess 
of 0.020 inch. 


5.2 Mapping Over-all Test on Burnback vs Arc Duration Axes 


The results of the over-all test program can best be discussed with the 
aid of Fig. 8. On this figure, the over-all test range is mapped into five 
regions and the monitoring limits visually displayed with respect to 
these regions on a burnback versus arc duration axis. Some further 
explanation will be helpful in interpreting this figure. 

For any weld, the average wire velocity during arcing and the dura- 
tion of the are may be measured. These measurements determine a point 
on Fig. 8, where the product of the measurements (burnback) is the 
ordinate and arc duration is the abscissa. As wire velocity is varied and 
a number of points are plotted, a curve takes shape. This curve relates 
the melted wire distance to the time after arc initiation for a particular 
set of circuit conditions and wire tip. Accordingly, the upper curve of 
Fig. 8 is the sketched-in melting curve for the welding conditions in the 
first phase of the test program. The lower set of curves is for Phase II, 
where each curve in the set would be determined by a corresponding 
value of circuit resistance. The intersection of a particular melting curve 
with a line of constant velocity locates the arc duration and burnback 
of a weld for the welding conditions corresponding to that melting curve. 
This point will be referred to as an “operating point” of the welding 
process. Naturally, there will be a statistical scattering of the data about 
the operating point due to normal test variations. 

The first test phase extends over regions A, B, and C while the second 
test extends over regions D and E. Note that in Phase I the operating 
points coincide with a melting curve, while in Phase II the operating 
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Fig. 8 — Graphical representation of over-all test program. 


ta(min) limit therefore is obviated in this case and the entries in the table 
agree with the burnback column. On the other hand, in region A, weld 
quality is poorer, despite an average burnback of 0.022 inch. Here 


ta(min) 18 the controlling limit. 


It is significant that there is this distinction in the controlling lower 
limit of arc duration, and the question is asked, Why does weld quality 


Weld Quality 


Test Phase Region Order of 


TABLE III 


ta < tatmax) 


Merit 
I A 3 
I B 1 
I C 4 
II D 5 
II iE 2 





* See Section 5.2. 
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drop off in region A? A comparison of welding conditions leads to some 
speculations as to the reasons. In region A, the electrodes close at the 
highest current density and with the fastest wire speed. It is reasonable 
that extremes in these variables could be undesirable. For instance, dur- 
ing the high current portion of the arcing period, the melting process 
and temperatures are undergoing rapid transitions. Probably some 
vaporization is taking place. If the wire speed is fast, closure occurs dur- 
ing this high current portion of the arcing period. Under such volatile 
conditions, the weld process is apt to be erratic. A relatively slower 
speed would lengthen the arc duration and allow conditions to become 
more stable. 

In general, poor weld quality will result if the wire speed is allowed to 
become too slow. This is seen clearly in Fig. 8 by letting the operating 
point move to the right along the upper melting curve. By the time the 
wire speed has decreased to 63 inches per second, the ‘operating point”’ 
is well into the poor area in region C. However, if we now hold the wire 
speed constant and add 0.060 ohm to the welding circuit, the operating 
point shifts into region E, and to a different melting curve, resulting in an 
improvement in weld quality. Hence, a change of circuit resistance in 
this case will restore the good weld quality which had been lost by slow- 
ing the wire speed. This suggests that the slow wire speed is responsible 
for the poor weld quality, not through any mechanical effects, but rather 
by lengthening the duration of the arc past a maximum safe limit for a 
particular melting curve. 


5.3 Statistical Interpretation of Data 


The data of Fig. 6 are replotted in Fig. 9. There are seven distribu- 
tions shown, one for each 50-microsecond increment of are duration. 
In general, all the distributions consist of two portions. The first portion 
is the distribution of welds which failed near the weld at less than wire 
strength, while the other is the distribution of welds which failed in the 
wire. In each case the portions representing weld failures are fairly 
linear, with a steep slope; however, in some cases the lines curve as they 
join the wire break portions. Because of the choice of ordinate, the por- 
tions representing wire failures are horizontal straight lines at 100 per 
cent of wire strength. 

To aid in interpreting these curves, the wire strength is thought of as 
an upper cutoff boundary on the weld strength. The following example 
is an analogous situation: suppose pennies were being pitched to a line. 
The density of pennies would be normally distributed on either side of 
the line. Now suppose that an adhesive-coated wall were placed at the 
line, all other things remaining the same. All the pennies which would 
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Fig. 9 — Cumulative distribution of data from test No. 1. 


ordinarily pass over the line now hit the wall and stick to it. A portion 
of a normal distribution is left to one side of the wall and all the remain- 
ing pennies stuck to the wall. Plotting the distribution of pennies around 
the wall would result in a distribution similar to that found in the case 
of percussive-welded connections. The pennies which stick to the wall 
represent wire breaks, and the truncated normal distribution of pennies 
to one side of the wall represents weld breaks. 

It appears, therefore, that the distributions of weld strengths in Fig. 
9 are normal below some value of weld strength since the lower strength 
portions of all the curves are essentially straight lines. Therefore, the 
probabilities associated with the lower tails of the distribution can be 
read directly. The two curves which represent the welds satisfying the 
monitoring inequalities approach the 0.01 percentage point of the dis- 
tribution at approximately 76 per cent and 80 per cent of wire strength. 
All the other curves are 45 per cent of wire strength or less at this equiv- 
alent point on the distributions. Thus these data indicate that, with 
monitoring, there is a much smaller probability of finding weak welds 
than in an equivalent size sample of unmonitored welds. 
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Based on these data, it appears that percussive welding is capable of 
meeting present quality objectives of a maximum defect rate of one in 
ten thousand, a defect being defined as a connection which is weaker 
than 60 per cent of wire strength. 


5.4 Generalizations Regarding Weld Strength Data 


The following are some typical characteristics of the strength dis- 
tribution of low-voltage percussive welds. These generalizations are not 
drawn solely from the present data, but have been observed repeatedly 
in various sets of data. 

(a) When the process is optimized, the distribution of weld strength 
is bunched tightly around the wire strength. 

(b) When the process is off optimum, a long tail develops in the dis- 
tribution of weld strength. This tail is made up of relatively weak welds 
which usually comprise a small percentage of the total number of welds 
in the distribution. The tail is an incomplete normal distribution, 
bounded on the upper side by the wire strength. 

(c) The over-all quality of a population of percussive welds is defined 
completely in terms of the tail of the distribution. 

The following factors all tend to limit an effective assurance of quality 
based on small sample size destructive testing: 

(a) Only a very small percentage of the data yields useful information 
regarding the tail of the distribution. 

(b) The large standard deviation of the tail makes large sample sizes 
necessary if confidence in the sampling is not to be sacrificed. 

(c) The allowable percentage of defects is small (one in ten thousand). 
Therefore, large sample sizes are required in order to gain a reasonable 
degree of confidence in the sampling. 

The conclusion, based on these factors, is that small sample size de- 
structive testing alone cannot assure quality levels of less than one 
defect in ten thousand. The sample size required to assure this degree 
of connection reliability, with a reasonable degree of confidence, would 
be prohibitive from a cost standpoint. 


VI. SUMMARY AND CONCLUSIONS 


1. It is characteristic of the low-voltage percussive welding process 
that the relatively weak welds in a population will be contained in a 
long tail of the distribution. The welds in this tail are distributed nor- 
mally, with large variance, and are bounded on the high side by the 
wire strength. The vast majority of welds are not contained in the tail 
and have essentially the distribution of the wire strength. This double 
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nature of the over-all distribution, unfortunately, severely limits the 
effectiveness of small sample size destructive testing. 

2. In view of the difficulties encountered in statistical sampling by 
destructive testing and because of the lack of any effective and practical 
method of testing nondestructively, it was sought to attain consistently 
reliable percussive welds by continuous control of the process, i.e., by 
monitoring certain important welding parameters as each weld is made. 

3. Preliminary tests indicated that the time duration of the arc was 
the most sensitive over-all indicator of variations in the process. The 
effectiveness of are duration is further increased if the measurement is 
related to the approach velocity of the wire. Since both these param- 
eters are relatively easy to measure, arc duration and wire velocity were 
chosen as the most effective and practical measurements for purposes of 
monitoring. 

4. It was found that there is a range of arc duration for which the 
process is optimized. Weak welds which occur at longer arc durations 
are attributed to excessive heat conduction losses which result in prema- 
ture solidification of the molten metal film over the arcing electrodes. An 
analytical derivation is made for a maximum allowable arc duration, 
which is defined as that time when heat conduction losses begin to ex- 
ceed power input to the weld. 

5. At excessively short arc durations, weak welds can be attributed to 
insufficient arc energy, which results in insufficient burnback on the 
wire to produce a full-area weld. Accordingly, the minimum allowable 
arc duration is shown to be the time required for the moving wire to 
travel a distance equal to its tip length. 

6. In addition to satisfying the proper limits of arc duration, other 
conditions should be satisfied in order to attain good welding results. 
These include a properly and consistently made wire tip, no excessive 
spreading of the arc, and sufficient force between electrodes at contact to 
affect good mating and fusion. 

7. A test program of nearly 6000 welds confirmed that the low-voltage 
percussive welding process is optimized when the arc duration is between 
an upper and a lower value. Furthermore, these limits, as indicated by 
the data, agree with the predicted values from the analysis and earlier 
testing. 

8. Based on the data of weld strength for welds satisfying the moni- 
toring inequalities, the low-voltage percussive welding process is capable 
of meeting present quality assurance objectives of less than one defective 
connection in ten thousand, where a defective connection is defined as 
one weaker than 60 per cent of the wire strength. 
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A Theorem on the Distribution of Weights 
in a Systematic Codef 


By JESSIE MACWILLIAMS 


(Manuscript received September 4, 1962) 


A systematic code of word length n is a subspace of the vector space of all 
possible rows of n symbols chosen from a finite field. The weight of a vector 
ts the number of tts nonzero coordinates; clearly any given code contains a 
certain finite number of vectors of each weight from zero to n. This set of 
integers is called the spectrum of the code, and very little is known about tt, 
although zt appears to be wmportant both mathematically and as a practical 
means of evaluating the error-detecting properties of the code. 

In this paper tt is shown that the spectrum of a systematic code deter- 
mines uniquely the spectrum of tts dual code (the orthogonal vector space). 
In fact the two sets of integers are related by a system of linear equations. 
Consequently there is a set of conditions which must be satisfied by the 
weights which actually occur in a systematic code. If there is enough other 
information about the code, it is possible to use this result to calculate its 
spectrum. 


In most systems of error correction by binary or multiple level codes 
the minimum distance between two code words is an important parame- 
ter. (The distance between two code words is the number of coordinate 
places in which they differ.) Much attention has been given to devising 
codes which have an assigned minimum distance. 

The weight of a code word is its distance from the origin. The distance 
between two code words is the weight of the vector obtained by sub- 
tracting one from the other, coordinate by coordinate. If the code words 
form a vector space, this vector is itself a member of the code. Such codes 
are called systematic codes. The set of integers specifying the weight of 
each code word is then exactly the same collection of numbers as the set 
of integers specifying the distance between each pair of code words. 


+ This paper formed part of a thesis presented to the Department of Mathe- 
matics, Harvard University, in partial fulfillment of the requirements for the 
degree ‘of Doctor of Philosophy. 
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Thus it is customary to talk about weight properties rather than distance 
properties of systematic codes. 

In many cases, practically all that is known explicitly about the 
distribution of weights in a code is that the weight hasa certain minimum 
value. Recent studies have shown that it would be useful (e.g., in the 
study of real life channels) to have more information. We would like to 
be able to answer questions of the following sort: 


z. Given a method (implemented or theoretical) for constructing a 
systematic code, how many elements of each weight will be obtained? 
(It is a safe assumption that nobody will want to write out the code 
vectors and count them.) 


aw. Given a set, U1, Ue, -°* , Us, Of positive integers, is it possible to 
construct a systematic code with elements of these weights only? 


In theory there exists a method of answering these questions.’”" Un- 
fortunately this method is quite difficult to apply. The purpose of this 
paper is to give a different method which is in some ways more useful. 

We show that the spectrum of a systematic code determines uniquely 
the spectrum of the dual code (the orthogonal vector space). In fact, 
the two sets of integers are related by a system of linear equations. Our 
main theorem shows how to obtain this system of equations. 

In Section I we give definitions and statements of the main theorem 
and of some corollaries. Section II contains proofs of these theorems. 
Section III describes how the results of Section I may be applied. 


I. DEFINITIONS, NOTATION AND A STATEMENT OF THE MAIN THEOREM 


Let F be a finite field of g elements; q is a prime power. Let F” denote 
the direct sum of n copies of F’. F” is the set of all possible row vectors 
of length n, in which each coordinate is an element of F. Addition of two 
vectors is defined coordinate by coordinate, under the rules prevailing 
in F, 

F” is a vector space of dimension n over F’. Choose a basis consisting 
of the n vectors 


= (1, 0, 0, +++ , 0) 


An element u of F” is then expressed uniquely as 
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nr 
u= >. Wi €:, u,eF, 
t=1 


We write u = (tt, Ue, °+* , Un). 

The weight of u is defined to be the number of wu; which actually ap- 
pear in this sum —1i.e., the number of nonzero coordinates in the vector 
u. 

An alphabet is any subspace of F”; a vector belonging to the alphabet 
@ is called a letter of @. It may happen that every letter of @ has zero as 
the jth coordinate — this case is not excluded. 

The scalar product of two vectors, 


n nr 
UuU= > ees v= >> ve, Uz,v0iéF, 
i=1 i=1 


is u#v = >,7-1 Uv; , Where the multiplication and addition are carried out 
in F. If F is the field of two elements 0, 1, for example, the scalar product 
of (1, 1, 0) with itself is 1-1 -+ 1-1 = 0. 

Two vectors u, v are orthogonal if their scalar product is zero. In the 
example above, (1, 1, 0) is orthogonal to itself. 

The orthogonal complement of an alphabet @ is the set of all vectors 
of F” which are orthogonal to every vector of @. It is clear that these 
vectors also form an alphabet, say @, which is called the dual alphabet 
of @. If k is the dimension of @, the dimension of ® ism = n — k. 

The main result of this paper is as follows (the proof is given in Section 
IT): 

Let @ be an alphabet of dimension k, and @ the dual alphabet of di- 
mension m. Let A;, B; denote the number of letters of weight 2 in @, ® 
respectively. Of course, An = Bo = 1. Set y = q — 1. Let 2 be an in- 
determinate. 

Theorem 1: The quantities defined above are related by the equation 


dA ++ ye)” *(1 —z)'= qi’ Dd, Bz’. 
Remarks: 


t. The formula above is symmetric in the sense that, setting (1 — z)/ 
(1 + yz) = 4, we obtain by straightforward algebra 


dX By. + v2)" “(1 — 8)* = q” a Aé’. 
wi. Theorem 1 is a statement about equivalent classes;*" it is still 


true if @, ® are replaced by equivalent alphabets. 
An alphabet @ is said to be decomposable,“ with respect to the basis 
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a, €, ‘+: , & Of PF", if it is the direct sum of two alphabets Gi, Gz, 
where @, contains n columns of zeros and @, occupies these columns only. 
For example, the alphabet 


0 0 0 0 
1 1 0 0 
0011 
1. 4. 11 


is decomposable, with 


000 0 000 0 


OS fA. 020 2=9 011. 


In general, @; is a k,-dimensional alphabet in F'"’, and @2 a kp-dimen- 
sional alphabet in F™, with ny + m. = n, ky + ke = kiki Sm, ke S ne. 

The dual alphabet of a decomposable alphabet is also decomposable; 
in fact ® = @ + B®, where @B; is the dual alphabet of @; in F“,7 = 
1, 2. (The example above is self-dual.) 

Corollary 1.1: If @ is decomposable, say @ = @, + @:, the equation 


2, As aS yz)" “(1 = a = q° oz Bez 


is reducible in the obvious sense; the factors are the equations pertain- 
ing to @;, Bin F™,¢ = 1, 2. 
For the example above we have 


(i+z2)°+ d-—2)P= 214 2)% 


Corollary 1.2: A necessary condition for the existence of an alphabet 


containing letters of weights w;,7 = 1, 2, --- , s, and no other, is that 
there exists a set of integers a;, 7 = 1, 2, --- , s, such that the expres- 
sion 


(Vs ye)" ay dX ele) le) 
I= 
when expanded in powers of z, takes the form 
f +d 2 Be’ 


where the 8; are positive integers. 
Unfortunately, this condition is not sufficient. For example, 


(i+2)°+714+2)°1 —2’+7(14+2)(1—2)°+ (1— 2)° 
= 21472 + 72° 4+ 2°), 
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but it is not possible to construct a binary alphabet containing 7 letters 
of weight 2 and no letters of weight 4. 

If Ay = Ag = --+ = Ag; = O, every vector of weight <7 in F” appears 
as a coset leader for @, and conversely. Another way of saying this is 
that, for all pairs of distinct letters a, a’, of @ and any 2 S j, the set of 
vectors at distance 7 from a is disjoint from the set of vectors at distance 
t from a’. In this case we can enumerate these vectors by weights as 
follows: 

Let f;,; denote the number of vectors of weight s in F” which are at 
distance z from some letter of @. Write 


(t+ 2)U = 2) = Dwi 
Corollary 1.3: If @ contains no letter of weight < 27 + 1, then 
DY fsa” = D Bw, +72)" = ay 
The proofs of corollaries 1.1 to 1.3 are given in Section II. 


II, PROOFS 


If @ is an alphabet of f'” and 6 the orthogonal complement of @, the 
weights of the letters of ® are, of course, uniquely determined by the 
letters of A. However, a much stronger statement can be made: the set 
of integers specifying the number of letters of each weight in @ is related 
by a system of linear equations to the set of integers similarly defined for 
@. This section will consist of proofs of this statement and of some of 
its consequences. 

Two proofs are given. The first is short and easy; the second is longer 
and more sophisticated. However, it incidentally produces a more gen- 
eral result and gives some insight into what is going on. 

We make the following conventions for notation: @ shall be a k-di- 
mensional alphabet in F”; ® shall be the orthogonal complement of @ 
of dimension m = n — k; y shall denote the quantity q — 1. A;, B; 
denote the number of letters of weight 2 in @, ® respectively. The bi- 


nomial coefficient (") is understood to be zero if s > r. 


Let 6, €,°--, €n be the usual basis of F”. Let s = (81, 82, °°, 8) 
be a set of v different indices, 1 < s; S n, and let t = (4, th, -+: , tr) 
be the complementary set of indices. Denote by F,’, F;" ” the spaces 
generated by e,, °°: , €s, and e:,, ---, €,_,. Clearly, F,’, F;” ” are 
orthogonal complements in F”. Let | H | stand for the number of vectors 
in a space H. 
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Lemma 2.0: 
eae Qur ecg eae ls 


Proof: The orthogonal complement of @ M F,”” is the smallest space 
containing ® and F,’. This is the lattice theoretic union of 6 and F,’, 
which we write ®@ U F,’. Then 


|@ UF,’ | : \anF,”” | = q” = gq”. 
The number of vectors in ® U F,’ is q’"q’/| @N F,’ |. 


Hence 
gv’ iank.”’|=q""|aNF,’ | 
or 
lan F"| =" |@NF,’ |. 

Denote by {(e,, ---, €,), a} a pair consisting of » basis vectors of 
F” and a vector a of @ which is orthogonal to each of ¢,, +++, é, 

Lemma 2.1: 

t. For a fixed set of indices s,, --- , s, the number of pairs 


{(é,,°°', &,),a} is |@NF,’ |. 
a. The total number of such pairs for all choices of v distinct basis 


: —71 
vectors is vi Ai ( , ) : 


Proof: 

i. F,”” consists of exactly those vectors of F” which are orthogonal 
to &,, °°: , &, 3 hence @f F,”” consists of exactly those vectors of @ 
which are orthogonal to ¢,,,°-: , &,. 

au. If a e @ is of weight z, then a is orthogonal to n — 7 of the vectors 
&, °*+,é€,- A set of » vectors may be chosen from these n — 7 in 


(" > 2 ways. Hence the total number of pairs 


(és.s -+*, 6s,)5 a} is Sa), 


Let >., indicate summation over all possible choices of v indices 
Si, °°*, & 3 similarly, >>, denotes summation over all the comple- 
mentary sets 41, ---, tn. Lemma 2.1 is equivalent to 


Slanre| ie eee 


WEIGHTS IN A SYSTEMATIC CODE 85 
Replace @ by @, vy by n — vy and s by t. The same argument then gives 


Den F | =e a 


n—v 


. n—1 _ pee (n-4 
> As ( : \=4 ee): 


Proof: For a fixed set s (which determines, of course, a fixed set t) 
we have by 2.0 


Lemma 2.2 


\aNF’”’|=¢’|@NF,’ |. 
Thus 
DAeNEe | Se" ener | 
t Ss 


which, by 2.1 is the same thing as 


¥4("> ')=e7 ha ("7 ?). 


The equation of 2.2 holds for vy = 0,1, --- ,m — 1. This is, in fact, 
one form of the promised set of linear equations between the quantities 
A;, B;. 

We now give the second proof. 

Let ¢ be a finite Abelian group. A character x of S is a homomorphism 
of G into the multiplicative group of complex numbers of absolute value 
1. The characters of G form a group §* which is isomorphic to G, there 
being in general no canonical isomorphism. t 

If @ is a subgroup of G, the characters such that x(a) = 1 for all a of 
@ form a subgroup ®&* of ¢*. ®* is precisely the character group of S 
mod @Q. 

Suppose now that G is the additive group of a finite field. S* is just a 
multiplicative copy of S, and the characters can be labeled by the ele- 
ments of G in a symmetric way; that is, if r, s, --- are elements of S 
we have 


xXr(S) = x(r) = x(r, 8). 


If G is the additive group of a prime field of order q, we take r, s etc. 


{ For prime fields, the proof can be given without mentioning the word charac- 
ter. The presentation here is an uneasy compromise with conscience — we wish 
to indicate possible extensions to nonprime fields without doing too much work. 
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to be the integers mod q, and set x(r, s) = ¢" where ¢ is a primitive gth 
root of unity. 
We have from the general theory of characters 


x(r,0) = x(0,s) = 1, 
. 
>» x(r,8) = —1 ifs ~ 0. 
r=1 


Let a, €,°*°, én be the fixed basis of F”, and u = (1, Ue, «++, Un) 
the coordinates of a vector of F” with respect to this basis. The char- 
acter group of I” is, of course, a multiplicative copy of F”. We label the 
characters by elements of F’” as follows: 


valv) = TE x(us, %6) = valu) = Wu, 0). 


Let @ be a subspace of F”. The characters such that y%(a) = 1 for all 
a of @ form a subgroup ®&* of the character group. &* is exactly the 
character group of F” mod @. The elements b which label these charac- 
ters form a subspace 6 of F”, isomorphic to F” mod @. In our notation, 
the equation y(a, b) = 1 holds for all a of @ and all b of @, and given 
either @ or ®, the other is uniquely{ determined by this condition. 
Lemma 2.3: Let @, ® be related as above. Then 


4. > v(v, a) = gifve®. 
ae@Q 

it. >» vv, a) = Oifve®. 
aeQ 


Proof: Part 7 is obvious, since by definition Y(v, a) = 1 if v e B. For 
it we observe that for a e€ @, ¥(v, a) = Wa(v) is a character of F” mod B. 
If 6 denotes a coset of F” mod @, Yoace Wald) = 0 for 6 ¥ B. Now 
V.(v) = wW.(d) for any v in 5; hence 


dv, a) = 2 ¥a(5) =0 if vee. 


Lemma 2.4: If F ts a prime field, then @, 8 are related as in 2.3 af and 
only uf they are orthogonal complements. 

Proof: y(a, b) = ¢=%°' where ¢ is a primitive qth root of unity. Hence 
v(a, b) = 1 implies that a is orthogonal to b. Since @ is isomorphic to 
F” mod B the dimensions of @, @ add up to n. Thus B is the orthogonal 
complement of @. 


+ That is, if F is a prime field. Otherwise we must fix the basis of F over its 
prime field before we claim uniqueness. 
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Let f(z, s) denote a function of the integers 7, s with values in a ring 
R. The values of f(z, s) may be added and multiplied, and these opera- 
tions obey the two distributive laws. f(z, v;) denotes the same function 
of 7 and the zth coordinate of v. 

Lemma 2.8: 


YG 0 = S160. 
Proof: If n = 1 the statement is 
Ds, 8) = DIL), 
which is obvious. Assume the truth of the lemma for FP”. 


Let F.”,0 < r S y, denote the set of vectors of F” which have last 
coordinate r. Clearly the F,” are a partition of F”. Then 


Sse) = EY [L1G v9 1n, 9 | 


ve Fr {=1 r=O0 ve F 


= s(n, Ds 2 ve v;) 


= p2 f(n, 7) 1 > HG, r) (by induction) 


n Y . 
= IT X FG). 
Let 2” be a set of (commuting) indeterminates, r = 0, 1, ---, y. To 
each vector v = (v1, v2, «++, Un) of F” associate a monomial [J,%1 2°”. 


The monomial associated with v describes how many times each field 
element appears as a component of v. Let R be the ring of polynomials 


in ras ial ---,2™ over the complex numbers. Let wu = (1, Ue, °°* , Un) 
be a fea ee of F”. 
Lemma 2.6: 
(%1) (v9) (vp) 
Tvs a)etP2? «5 2 = TT DY x (us, 2”. 
veE j=1 r=0 


Proof: Set f(j, vj) = x(u;, Um)z?, which is in R. Then 
vu, verre” ... 2 = TT fi, 05). 
j=l 


By 2.5, 
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Eves opeP2? 2 = TD 9G, 2) 


j=1 r=0 


= TLD xtu;,r)2”. 


j=1 r=0 


Lemma 2.7: Let @, @ be orthogonal complements in F”, as usual. Then 


pS iat 3 x(u; rye” = q p27 gry (v2) , tek gen) 


ueqQ jel r= 


Proof: We evaluate the quantity 
F(u, v) = 2s. 2 x(u, vez (v2)... pn) 


in two ways, which give the two sides of the equation. 
By 2.6 


F(u,v) = do Il 3 x(u;, re. 


ueQ j=l r=0 


Also 
Fu, 0) = So eet? «.- 2 YY vu, v). 
ve Fn UE 


By 2.4 and 2.3 


DD Wu, v) = 


wuweER 


| ifveB 
0 ifves. 


Hence 


F(u, v) = g’ ¥ geVy2) ... gn) 


».EB 


Theorem 2.8: Let @ be a k-dimensional alphabet of F”, and @ the ortho- 
gonal complement of dimension m = n — k. Let A; , B; denote the number 
of letters of weight 1 in @, @. Then 


> A(l+ya)"*1—-2)i'=¢ p> Bye’. 


2 irs 0 
Proof: In 2.7 set 2 = 
1 afr =.0. 


If uz = 0, x(uj, 7) is 1 and Do% x(u;, r)z” becomes (1 + yz). 
If u; #0 >0%1x(u;,7) is —1, and reo x(u;z, 12 becomes (1 — 2). 
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Let | u| denote the number of nonzero u;. 
Then [[1 > 70 x(u;, rz goes into (1 + yz)" !“(1 — 2)", 


| wu | is of course the weight of wu = (t1, Ue, °-+, Un), SO that the left- 
hand side of 2.7 becomes 


Dy Ad +12)" = 2) 
The right-hand side of 2.7 is clearly 
q > Be, 
which proves the theorem. 
Innumerable sets of linear equations between the quantities A;, B; 


may be obtained from theorem 2.8. The following two are sometimes 
useful. 


Lemma 2.9: For vy = 0,1, ++: ,n, 
: > n-it\ wr n—t 
t. = As ( F \= Se). 


(These are the equations of 2.2.) 


oa 4 k—v i wi f{/n— 4 

t1. o 4:(*) = g do (-1) By Co 

t. is obtained by setting (1 + yz)/(1 — 2) =1+y, 

aw. by setting (1 — z)/(1 + yz) = 1+ y. The algebraic details are 
easy to verify. 

This process is reversible, i.e., (7) or (47) imply 2.8. Before exploring 
the consequences of theorem 2.8, we give a different specialization of 2.7. 

Theorem 2.10: Let B, be the number of letters of ® which contain s 
coordinates equal to 1. Let Ao; be the number of letters u in @ of weight s for 
which >>." ui = 0. Let Ay, be the number of letters win @ of weight s for 
which >> ;21 us ¥ 0. (Clearly As = Acs + As). 
Then 


n 


a Bers = (Au — As/y)(2 — 1)*(2 +)". 


s=0 
Proof: In 2.7 set 
Zz ifr = 1 


I ie I. 
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Then 


n 

k k 1 

g Dd nee eek glen) becomes fi » B,‘ gs 
ve s= 


Y Y 
Dd x(uz, rz becomes x(ui, Lz + D> x(ui,r) — x(ui, 1) 
r=0 r=0 


= x(u,1)(2 =D) + Yex(ass7) 


x (U3 , 1)(z =— 1) if u; ¥ O 
ety if u; = 0 


n 


Y 
Il Sx (ui,rz becomes (2 + vy)” (2 — 1)™! If xu, 1), 


t=1 r=0 
Now if uw is a letter of A, so are also the letters 2u, --- , yu, and these 
have the same weight as u. We sum first over these letters 


3 IE x (sus) =D x(s2u, 1) 


Y if Zu; = 0 
The sum of 
(2+ 7)" i" — 1)" I] x(us, 1) 
over all letters in @ of the same weight as wu is thus 


(Agia — (Atay) (2 +)" "ze — 2)". 
Hence the left-hand side of 2.7 becomes 


3 (dee = (An /))(@ +9)" = 1) 


We return now to the consequences of theorem 2.8. As remarked in 
the proof of 2.10, if @ contains a letter u, it contains also the letters 
2u, ---, yu; that is, the number of letters of weight 2 in @ is divisible 
by y forz > 0. We have then 

Lemma 2.11: A necessary condition for the existence of an alphabet con- 
taining letters of weights w;,7 = 1, 2, ---, s, and no other, ts the existence 
of a set of integers a;,t = 1, 2, --- , s, such that the expresston 
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(1 + yz)" + y¥ 2 ai( ye) teh ie 


when expanded in powers of 2, takes the form 
q +d 2d, Be’, 


where the B; are positive integers. 

It has been pointed out before that this condition is not sufficient. 

Suppose now that @ = Q@, + @: is a decomposable alphabet. @; is a 
k;-dimensional alphabet in F”’, 7 = 1, 2, with orthogonal alphabet B;. 
ki + ke = k, and m + m = n. Let A;”, A; and B;”, B;® be the 
number of letters of weight 7 in Qi, @:, 1, ®. 

Lemma 2.12: 


> Adl + y2)" “(1 — 2)° 
= P As? (1 + yz)" "(1 — -)'| Pb AOL + yz)? "(1 - 2)'| 


i=0 i=0 t=0 
Proof: The number of letters of weight s in @; + @e is 
> B® Bp®, 


otp=s 
which is the coefficient of 2° in > 745 Biz’ So725 B;z". Similarly, 


>? A A 


o+p=s 


is the coefficient of (1 + yz)” “(1 — z)* in 
[Saat ey] |S aa +a - 2]. 
7=0 7=0 


We define the coset leader of a coset of @ in F” to be an element of 
least weight in the coset. The weight of a coset is defined to be the 
weight of its coset leader. 

If A; = Oforz = 1, 2, --- , 2e every vector of weight < e¢ in F” ap- 
pears as a coset leader for @ and conversely. Another way of saying this 
is: for all pairs of distinct letters a, a’ of @, the set of vectors at distance 
t S e from a is disjoint from the set of vectors at distance 7 from a’. 

Let «°”, oo, --- , ¢,“” be the cosets of A of weight i; we assume that 


92 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


y= r(?), 1e., that all vectors of weight 7 appear as coset leaders. Let 

fs,; be the number of vectors of weight s contained in the set-theoretic 

union U ¢;‘”. The polynomial >~”_» f.,;x° is called the enumerator (by 
j=1 


S 
weight) of this set of vectors. We propose to show that theorem 2.8 gives 
a convenient expression for this enumerator. We need the following 
preliminary lemma. 

Lemma 2.18: Let u be a fixed vector of weight 1. Let ds: be the number of 
vectors of weight s which are at distance t from u. Then 


yy, Du doy! = (1+ yay)" fa ty + (y — Iayl’. 


Proof: Suppose first that u = (4, We, -*-, Ui) is a vector of weight 
tin F*, We show that under these circumstances 


Dd da = lw ty + (y= VWayl’. 


This is obvious for 7 = 1; we suppose it true for? — 1. Letv = (u,,%, 


-++ , 0;-1) be a vector of weight s distant ¢ from (uw, we, °+* , Usa) in 
F*", From v we obtain: 
t. One vector (v1, %2, °°: , ¥i-1, 0), weight s, distant ¢ + 1 from wu. 
a. One vector (v1, v2, °°: , Vin, Us) Weight s + 1 distant ¢ from wu. 


ait. y — 1 vectors (11, v2, °° , Vi, Vi) Vs X 0, vi ¥ Us which have 
weight s + 1, and are distant t + 1 from wu. 

Hence the enumerator for 7 is obtained by multiplying that for 7 — 1 
by [x + y + (Wy — 1)zy], and the lemma is proved for n = 7. 

We now apply induction to n — 7. Let wu be a vector of weight 7 in 
F", and u’ a vector of weight ¢ in F”” obtained from u by omitting one 
zero coordinate. Let v’ be a vector of F”* which has weight s and is 
distant ¢ from u’. From v’ we obtain in F” 

z. One vector of weight s, distant ¢ from u, by adding a zero coordinate 
to v’. 

a. y vectors of weight s + 1, distant ¢ + 1 from u, by adding a non- 
zero coordinate to v’. 

This corresponds to multiplication by (1 + yxy). Hence the lemma is 
proved. 

Lemma 2.14: Suppose that A; = 0,7 = 1, 2, +--+ , 2e, and take t S e. 
Then the enumerator, >.7-0 f.,1v°, of vectors in cosets of weight t of F” 
mod @ is the coefficient of y’ in 


a Axl + yay)” “(a2 ty + (y — Lay)’ 
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Proof: The cosets of weight ¢ in F” mod @ are disjoint, and contain all 
vectors of /” which are at distance ¢ from some letter of @. 

Set (1 + yz)" “(1 — 2)’ = Doty WG, n, r)2". Let ®, B; have their 
usual meaning. Assume the conditions of 2.14. 

Lemma 2.16: 


q” d, hat = > Bi WV(i, n, t)(1 + yx)" “(1 — 2)’ 





Proof: Set 
patty t Wy — ley 
1 + yay , 
then 
(1 + yx)(1 + vy) (l—2)(1 — y) 
1 as ee i o Yi 
1s (1 + yay) , : 1+ yay 


Make this substitution in the equation 


Dd Bal + ye)" “(1 — 2) = a" D Ae! 


we obtain 
2 Bal + ye)" = 2) tay = 9) 

= a" XA + vey) "e + y + (y= Day)’ 
Equating coefficients of y" gives us 


2, BUC, m,t)(L + ya)" — 2)" =” De fat’ 


III. APPLICATIONS 


The easiest application of theorem 1 is to a generalized Hamming 
alphabet, that is, a close-packed 1-error correcting alphabet over a field 
of q elements. Such an alphabet exists for n = (q” — 1)/y, all m > 1. 

The dual alphabet is of dimension m, and contains (q” — 1) letters 
of weight g”’. The spectrum of a generalized Hamming alphabet is 
thus given by the expansion of 


(Pye eg” 1 eye Ch Se 2) 


where n = (q” — 1)/y,u = q” 





7 j A similar formula (¢ = 2) is found by Lloyd® for close-packed codes which 
are not assumed to be group codes. 
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TaBLE I — DIstRiBUTION OF WEIGHTS IN THE Two GOLAY CODES 

The first table is for the 3-error-correcting (23, 12) alphabet over Z-2 , the second 
for the 2-error-correcting (11, 6) alphabet over Z; . In both cases, 7 stands for 
weight, B; for the number of letters of weight 7 in the dual alphabet, A; for the 
number of letters of this weight in the Golay alphabet. 


t B; A; 
0 1 1 
7 0 23 X Il 
8 23 X 22 23 X 22 
11 0 23 X 56 
12 23 X 56 23 X 56 
15 0 23 X 22 
16 23 X ll 23 X iil 
23 0 1 
t B; A; 
0 1 1 
5 0 2X 66 
6 2X 66 2X 66 
8 0 2X 165 
9 2X 55 2X 55 
11 0 2X 12 


Theorem 1 may, in fact, be used to calculate the number of letters of 
each weight in any close-packed code. The results for the two Golay°® 
codes are given in Table I. 

Anything which is known about the structure of an alphabet or its 
dual may be used with theorem 1 to limit the number of possible weight 
distributions. Such items of information are a very diversified character, 
and no general method has been developed. However, the results ob- 
tained by hand calculation indicated that it is probably worthwhile to 
make a systematic computer study of the known classes of alphabets. 
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Surface Effects of Radiation 
on ‘Transistors* 


By D. 5S. PECK, R. R. BLAIR, W. L. BROWN and F. M. SMITS 
(Manuscript received October 10, 1962) 


Observation of surface effects of tonizing radiation on several types of 
transistors indicates that in reverse-biased devices these effects occur at much 
lower radiation dosage than in unbiased devices or bulk semiconductor 
material. Further, the total radiation dose rather than dose rate seems often 
to be the more important factor in the effect. The type of particle used in 
trradiation 1s unimportant; the significant factor 1s the tonization tt pro- 
duces. The effects seem to arise from tonization of gases within a transistor 
encapsulation and interaction between the ionized gas and residual semi- 
conductor surface contaminants. This results in inversion layers at the 
device surface and thus in alteration of junction characteristics. The changes 
in device properties are not permanent, but the recovery after removal from 
radiation is complex and proceeds with characteristic times between seconds 
and days. 

Different types of devices may respond quite differently to exposure, and 
the response 1s different even between different batches and individuals, 
indicating a dependence upon device processing. 


I. INTRODUCTION 


A wide variety of effects of high-energy radiation on semiconductor 
materials and devices has been recognized and studied for a number of 
years. The major emphasis in this field has been on effects that involve 
the bulk properties of semiconductors. A great deal of progress has been 
made in understanding the processes that control bulk radiation phe- 
nomena! and in understanding the implications of these phenomena for 
semiconductor devices.? Radiation effects on semiconductor surfaces also 
have been observed,’ and this paper is concerned with some special 
aspects of surface phenomena that have recently come to light. In con- 
trast to bulk effects, the surface radiation effects are very poorly under- 


* Partly supported by contract with Electron Technology Laboratories, Aero- 
nautical Systems Division of the United States Air Force Systems Command. 
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stood and in general even poorly characterized. This paper will attempt 
to shed a little light on these complexities and indicate a type of meas- 
urement program that has been found appropriate for dealing with 
devices intended for use in a radiation environment such as that of 
the Van Allen belts in space. It will fall short of providing a satisfactory 
understanding of the processes involved. 

The paper is arranged as follows: in Section II, for perspective, a 
brief discussion of the two broad classes of bulk radiation effects is in- 
cluded. Section III describes the early observations that provoked the 
present work. Section IV is a proposed model of the basic process. 
Section V describes results of a number of experiments carried out to 
test the mechanisms of the process. Section VI characterizes effects 
with significant numbers of devices. Section VII describes the process of 
testing and selection undertaken for Telstar devices. Section VIII is a 
summary of the paper. 


II, BULK RADIATION EFFECTS 


Bulk radiation effects can be placed in two broad classes that arise 
from (a) hole-electron pairs produced in the crystal by ionizing radiation 
and (b) defects in the semiconductor lattice produced by high-energy 
particles. These phenomena often occur together but they result from 
quite different interactions of radiation with the solid and they have 
very different consequences in semiconductors. 


2.1 Hole-Electron Pair Formation 


Fig. 1(a) illustrates the first case. Any charged particle passing through 
a solid produces ionization through collisions with the bound electrons. 
These collisions excite electrons to the conduction band and leave holes 
in the valence band, producing electron-hole pairs in exact analogy with 
the production of pairs by light. Neutrons and gamma rays also cause 
ionization effects through intermediate reactions that produce charged 
particles. As far as effects which depend on ionization are concerned, the 
particular particle involved is incidental; all that matters is how much 
energy is lost in the solid. The number of hole-electron pairs produced is 
proportional to this energy loss. 

The generated pairs tend to recombine with a time constant that is 
the conventional lifetime. Hence, all effects in this class are transient 
and persist only for the order of a lifetime after the excitation is removed. 
The pairs produced alter the conductivity of semiconducting materials. 
They also contribute currents in p-n junction diodes and transistors. 
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Under pulsed ionizing radiation these effects can alter conductivities 
and currents by many orders of magnitude. On the other hand these 
effects can be very small in response to a single energetic particle and 
special p-n junction diodes may be required even to detect them.’ A 
case of intermediate magnitude has been considered by Rosenzweig5 
who has used silicon solar cells to measure the intensity of moderate 
radiation fields. 


2.2 Lattice Damage 


Fig. 1(b) illustrates the other type of bulk radiation effect in semicon- 
ductors that arises from collisions of energetic particles with the nuclei 
of the lattice. If such a collision transfers sufficient energy to the struck 
atom, it is capable of moving the atom from its normal lattice site to 
some interstitial position in the crystal. These events are rare by com- 
parison with the ionization events of Fig. 1(a), but they create perma- 
nent or at least semipermanent defects in the structure of the lattice. 
The most important consequence of these defects is reduction in the 
carrier lifetime of the material. Increases in diode reverse current and 
decreases in transistor current gain are produced as functions of the 
time integral of the flux of particles. In contrast to the pair production 
of the paragraph above, this radiation damage is extremely dependent 
on the particular particle involved. Energetic protons, for example, are 
much more effective in producing damage than energetic electrons. 
This type of radiation effect is of major importance to the long-term 
power conversion efficiency of solar cells in space. Detailed considera- 
tion of this problem can be found elsewhere.® 


Ill. THE SURFACE PROBLEM 


Surface effects on semiconductor devices have an illustrious history of 
subtlety and perversity and it is no surprise to learn that radiation is an 
environmental factor that must be considered. Several years ago, before 
the present sophistication in surface processing, experiments were carried 
out in an attempt to characterize radiation surface effects? (changes in 
junction current, breakdown voltage, current gain, etc.). No systematic 
picture evolved, although surface cleanliness seemed certain to play 
some role. 

Transistors that have evolved from refinements in junction formation 
and surface treatment techniques have, in the last few years, been 
examined in a variety of high-energy neutron and gamma-ray environ- 
ments. Particular attention has been given to bulk radiation damage 
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Fig. 1. — (a) Production of hole-electron pairs by collision of a charged parti- 


cle with electrons of a semiconductor. (b) Production of lattice defects by colli- 
sion of particles with the nuclei of a semiconductor. 


effects on current gain since these effects turned out to be serious (par- 
ticularly for silicon) at the flux levels of current interest in the vicinity 
of nuclear reactors. Very substantial improvements in radiation toler- 
ance have been found; however, no recent comprehensive work on the 
surface effects has been reported. 

The relative stability of the characteristics of a particular type of 
diffused silicon transistor under radiation is illustrated in Fig. 2. The 
collector reverse current is displayed because it is a particularly sensitive 
indication of surface stability at the very low currents that are conven- 
tional in present silicon devices. The gamma-ray radiation* used in this 


* Radiation dose is frequently measured in ‘‘rads,’’ which are a measure of the 
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Fig. 2. — The response of I¢zo of a diffused silicon transistor to either radia- 
tion or bias alone or to both together. Radiation dose rate 8.5 X 105 rads/hr. 


case is representative of ionizing radiation in general. The figure shows 
that after one-hour tests with neither bias nor radiation or with either 
alone, the current is essentially unchanged. (Actually much longer tests 
under bias are conventionally made in checking device reliability with 
the same result as shown here.) However, when bias and gamma radia- 
tion are simultaneously applied two effects appear: there is a sudden 
rise in current and then an upward drift of current over many minutes. 
The sudden rise of about 0.02 ya is largely due to ionization of the gas 
in the device can. The current drifts upward by another decade in the 
45-minute exposure following its rapid rise. When the gamma radiation 
is removed the current does not immediately drop back to its pre-irradia- 
tion value, but gradually declines and even after an hour is more than a 
decade too high. A drop of about 0.02 va must still occur at gamma-ray 
turnoff, but this is so small compared with the level to which the current 
has drifted, that the drop is not visible.t 

The drift up and the slow decay of the collector currents are surface 
effects produced by the radiation only when the collector junction is 
reverse-biased. This aspect of radiation sensitivity had not previously 
energy absorbed per unit mass of material. One rad is equivalent to absorption of 
100 erg/gm. Such a dose would produce approximately 2 X 109 ion-electron pairs/ 
cm? in atmospheric air or about 4 X 10!8 electron-hole pairs/cm? in silicon. 

+ Ion currents to the electrical leads either in the device encapsulation or in 
the gamma radiation chamber pose a serious problem to measurements on de- 
vices at very low leakage currents in high radiation fields (105-10® rad/hr). It 
has often proved necessary to remove the device from the radiation environment 


momentarily for measurement. Because of recovery effects care must be taken to 
obtain measurements quickly and at uniform times after removal. 
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been reported. It represents a factor that must be considered for relia- 
bility of devices that must operate in any radiation environment includ- 
ing the high-energy electron and proton belts found in outer space. The 
conclusion that these effects arise at the semiconductor surface can be 
reached in a number of ways. If they were bulk effects the marked in- 
fluence of applied bias is quite unreasonable. Radiation defects created 
in a solid have been influenced by extremely high electric fields, but only 
to the extent of causing them to migrate very short distances in long 
times.’ Furthermore, the decay of the effect occurs in a time that is 
very short for defect annealing in silicon. But more convincing is the 
sensitivity of the effect to the surface environment of the device inside 
its encapsulating can. 

Fig. 3 shows four typical npn diffused silicon transistors of two types, 
each with two kinds of ambient atmosphere. Increases in collector 
current are observed in all four, but at quite different integrated gamma- 
ray doses. In both device types early current increases are associated 
with gas filling. In type B the evacuated device shows no measurable 
change until the integrated dose is in a range expected to cause sub- 
stantial decreases in bulk lifetime. In this case, the influence of bias, 
although not shown in the figure, is practically nonexistent. 

Not only is there a variation in the response depending on the gas 
filling of transistor cans, but there is also great variability among devices 
with a single type of filling. All evacuated units are slow to respond, 
but in the type B transistor only about half of the gas-filled units re- 
spond quickly. The other half are almost as stable as the evacuated 


TYPE A 
DRY AIR + GETTER-~. 








Ico AT 5VOLTS (4A) 











10 10° 10° 10° 10 10 10° 
INTEGRATED ¥ DOSE (R) 


Fig. 3. — The radiation degradation of Iczo of two types of diffused silicon 
transistors, evacuated or with gas filling. Radiation dose rate 8.5 X 105 rads/hr. 
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units. This points again to a surface effect and to a broad spread in 
response arising from rather subtle differences in the surface chemistry. 
This lack of reproducibility complicates the study of the process in- 
volved and necessitates the use of statistical experiments, some of which 
will be reported in later sections. 

Before continuing to discuss the experimental observations we will 
introduce a simple model which describes some, but not all, of the effects, 
and provides a framework for the later discussion. 


IV. A MODEL OF THE PROCESS 


Radiation, gas encapsulation, and device bias that seem to be essential 
to the effects shown in Section III can be combined in a simple model of 
the process. Fig. 4 illustrates the ingredients: the fringing field of a 
reverse-biased collector junction on an npn transistor and ions and 
electrons produced by gamma radiation in the gas of the encapsulation. 
The fringing field separates the electron-ion pairs, depositing electrons 
on the collector side of the junction and positive ions on the surface of 
the base. On both regions these charges tend to produce inversion layers 
at the surface, the effects being analogous for pnp and npn devices. For 
simplicity only the inversion layer on the p-type base of a npn transistor 
will be considered. A magnified view of the edge of the device might be 
as shown in Fig. 5. The positive ions induce an electron-rich inversion 
layer or ‘‘channel’’ on the base and in effect extend the collector region 





Fig. 4. — A model of a reverse-biased transistor in a gas atmosphere ionized by 
radiation. 
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e 





Fig. 5. — Formation of a channel on the base of a transistor by positive ion 
collection on the surface. 


out over the base. The channel represents a grossly different surface than 
existed before. Since the junction between the electron-rich channel and 
the base material constitutes an extension of the collector-base junction, 
it contributes to the collector saturation current. In part this is simply 
because of the extra junction area. More importantly, because the 
channel junction is very close to the surface, the surface generation 
process can yield much more current per unit area than for the junction 
in the bulk. Furthermore, if the channel extends to the emitter it can 
add additional current to the collector by serving as a conducting path 
between the two. Channel effects on transistors and diodes have been 
studied previously in considerable detail’ entirely unrelated to the 
presence of radiation. The basic channel characteristics of emitter-to- 
collector conductance, high emitter floating potential, and channel 
pinchoff have all been observed in connection with the present surface 
radiation effects on diffused silicon transistors. 

We have so far considered only formation of a channel by charge 
collection in the collector fringing field, but a second possibility is 
shown in Fig. 6. In all the npn diffused silicon transistors used in these 
experiments the collector is electrically tied to the encapsulating can of 
the device. Under collector reverse bias a field then exists throughout 
the whole can of such a sign as to drift positive ions toward the surface 
of the device base. This feature can be expected to increase the ion- 
collection efficiency. 
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Fig. 6. — The enhancement of ion collection at a transistor surface due to 
electron fields between the device and its encapsulating case. 


One may ask if there are enough ions produced in the gas of the can, 
by the radiation doses that have been used, to provide large channel 
effects. At an integrated dose of 104 rads (at which these effects may be 
substantial) and at atmospheric gas pressure (the normal device filling) 
a total of about 2 < 10% ions per cc will have been produced. With a 
device enclosure of about 3 X 107° cc and a base layer area of about 
10-* cm?, if all the ions were collected on the base, their concentration 
would be about 10/cm?. On typical base material only about 10" surface 
charges/cm? will be required to produce a channel. These several orders 
of magnitude margin are probably quite important because the efficiency 
of the surface charging process seems likely to be quite low. Furthermore, 
the process is far from a perfect charge integrator as we shall have 
occasion to observe in connection with reciprocity experiments in Sec- 
tion 5.2. 

The lack of reproducibility noted in Section III in connection with 
devices of the same type and same gas filling seems to necessitate an 
elaboration on the model. If ions of the gas are sufficient to produce a 
channel, why are there some devices with gas that are as stable as de- 
vices without? The gas ions themselves must not be the tenacious 
charge on the surface that forms the channel. The gas ions probably 
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exchange their charge with residual contaminants on the device surface. 
Ionization of the surface contaminants directly is apparently too rare to 
be observed, since vacuum-encapsulated devices show uniformly high 
surface stability. This is roughly reasonable since the probability of 
ionization of any single atom is estimated to be only about 10-5 at 104 
rads. Even with a monolayer of residual surface contamination, the 
surface ion concentration would then be only about 10°/cm?2. Of course, 
if the gas ions are to do the job, they must be reasonably effective in 
finding and exchanging charge with the residual impurities. There seems 
to be margin for inefficiencies in these very rough numbers. 

This model predicts a number of effects that can be tested: 

1. The effect depends only on ionization, not on incident particle type. 

2. The simplest form of the model suggests that the effect is cumula- 
tive and depends on total dose, not on dose rate. 

3. The effect should be more pronounced at higher collector bias. 

4. The electric field between the can and the semiconductor may in- 
fluence ion collection. 

5. The decay of the effect should be faster if the device is not under 
continuous bias and should be accelerated in the presence of radia- 
tion without bias. 

The results of the experimental tests of these predictions are contained 
in the following section. 


V. TESTS OF THE MODEL 


5.1 Lonization 


If ionization in the gas of the device encapsulation is essential, then 
just as in the case of the ionization effect in a semiconductor discussed 
in Section 2.1, the type of energetic particle should not matter. The ex- 
periments in Section III were carried out with Co® gamma rays which 
ionize through photo or Compton electrons that they produce. We have 
tested the model by comparing the results of 18 Mev proton irradiation 
with the Co® gamma rays. The individual device response scatters so 
widely that the behavior of a number of similar devices is examined in 
each case. Fig. 7 shows the collector reverse current versus radiation 
dose, for gamma rays in dashed lines and protons in solid lines. The dose 
is calculated simply from the amount of energy deposited in a gas (the 
gas of the encapsulation) by gamma rays and protons. The radiation 
intensity (dose per unit time) was approximately 10° rads/hr for the 
protons as well as for the gamma rays. 

Within the spread of response observed, there is no significant differ- 
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Fig. 7. — Comparison of I¢go changes under radiation by gamma rays (dashed 
lines) and protons (solid lines). 


ence between protons and gamma rays. If bulk damage in silicon were 
involved, one could expect a factor of 10 to 100 greater effect for the 
protons than for the gamma rays. (The bulk damage effect per particle 
comparing 18 Mev protons and 1.25 Mev gamma rays would be a factor 
of 10* to 105, but the scale in Fig. 7 is not particles but ionization, or 
energy loss, and the protons lose energy at a much higher rate than 
gamma rays.) We conclude that ionization is essential to the mechanism 
of this surface radiation effect. 


5.2 Reciprocity 


The surface ionization effects may depend only on the total radiation 
dose or they may also depend on the dose rate. In the first case there is 
reciprocity between dose rate and time. Information on this point is 
relevant to understanding the process and is of immediate practical. im- 
portance as well. The experiments are most easily performed at high 
dose rates, but the device reliability is also of concern in modest radia- 

tion fields. 
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In Fig. 2 recovery of the collector current following radiation is illus- 
trated. The fact that recovery occurs at all when a device is removed 
from radiation but kept on bias proves that the observed result does 
not depend solely on the total dose. The ionized state of the surface at 
any time is not simply related to the total number of ions that have 
reached that surface. 

Fig. 8 shows the change of collector current vs dose at two different 
high dose rate levels. The experiment was started at 8 X 10° rads/hr. 
After approximately half a minute, when a dose of 6 X 10° rads had been 
delivered, the transistor was placed inside a lead shield that attenuated 
the gamma radiation by a factor of six. When a dose totaling 1.8 < 104 
rads had been given the device, the attenuator was removed and the 
higher level radiation continued. After each attenuator change, the col- 
lector current continues to rise in a smooth extension of the earlier por- 
tions of the curve. If the dose rate were important, the middle segment 
of the curve would be expected to have a different slope than the two 
ends. During the attenuator changes the device was out of the radiation 
for approximately a minute and the current had started to recover, as 
indicated by the first point taken at the lower dose rate. A comparable 
drop would have appeared at the second attenuator change, but the 
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Fig. 8. — Reciprocity of dose rate and time at two high radiation dose rate 
levels. 
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earliest measurement after reinsertion in the radiation field was not 
obtained soon enough for it to show. 

The rapid reestablishment of a previous high response to radiation, 
after some recovery from it, represents a memory in the process. This 
kind of memory has been seen repeatedly. Devices that were irradiated 
and have apparently completely recovered their original collector charac- 
teristics by standing out of radiation and even off bias, will still tend to 
reestablish their former response on a second radiation exposure. Even 
after several weeks, a device seems to retain a sensitivity to subsequent 
radiation as a result of an earlier exposure. There is some evidence that 
the memory can be removed by baking at approximately 100° C for a 
few hours. 

The ionization produced by the radiation apparently has two func- 
tions: first, to produce some chemical species that are capable of ioniza- 
tion and second, to keep these species in an ionized state. The number 
of these centers would reflect the total radiation dose, but some mini- 
mum level of radiation would be required to keep them active. 

We can draw the conclusion from Fig. 8 that, in the high-intensity 
region, dose is the important variable. Fig. 9 illustrates quite a different 


10-6 


EMITTER FLOATING 


1077 


10-8 


I¢Bo !N AMPERES 


1079 


10 RADS/HR 





19710 
102 10° 10% 10° 
RADIATION DOSE IN RADS 


Fig. 9. — Lack of reciprocity between high and low dose rates. 
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range. Here a device started its radiation history at only 10 rads/hr and 
established a pattern of current increase that is already determined at a 
dose less than 10‘ rads. Putting the device into the attenuated high- 
intensity source does not produce a simple continuation of the earlier 
curve, but rather produces a curve two and one-half orders of magnitude 
higher in current. This illustration is extreme and not many devices 
show this large a discontinuity, but essentially none give results that 
could be interpreted as simple reciprocity. Apparently over a dose-rate 
range this great, using the concept of the preceding paragraph, the 
ionization-sensitive chemical entities are not fully ionized in the low 
radiation field. 


5.3 Surface I’ffect vs Collector Bias 


Since the surface effects we have been considering are absent except 
under the simultaneous application of radiation and collector reverse 
bias, we expect to find that changing the bias will alter the effect. This 
can occur because of increases in efficiency of charge collection, but also 
because the increased junction field tends to bind the ions more tightly 
to the surface or distribute them to form a more extensive inversion 
region. 

Fig. 10 shows the effect of a sudden change in bias from 5 to 15 volts. 
One gains the impression that the current suddenly adjusts to a level 
and to a rate of change that are what they would have been if the entire 
dose had been given the device at the higher bias. If this impression is 
valid it would indicate that the effect at higher bias had not been re- 
tarded by the initial dose at low bias, and hence that ion collection was 
no different. In a plot like Fig. 10 it must be realized that the significance 
of the first dose rapidly diminishes in comparison with the total as one 
moves out along the second branch of the curve. It does seem possible 
to infer that a given number of ions on the surface gives a larger current 
contribution at higher bias. Rearrangement in the new field and effective 
extension of the length of a surface channel can occur. 

Vig. 11 shows the response of transistors under radiation at different 
biases. Up to a dose of about 10* rads the current increase depends 
strongly on bias. Between 10‘ and 10° rads the slopes of the curves on 
this log-log plot are quite similar. Although only 4 units are illustrated 
in Fig. 11, the total group in the experiment was 80 and this pattern of 
behavior was quite consistent. Even beyond 10° rads where the slopes 
increase there is considerable uniformity between individuals and across 
the voltage range. This point is illustrated in a different way in Fig. 12 
where the median behavior of a group of 10 to 12 devices at each bias is 
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Fig. 10. — Influence of change in collector bias on degradation of Iczo . Dose 
rate 3. 5 xX 10° rads/hr. 


plotted at 6 X 10° rads and 6 X 10‘ rads. The pattern of bias dependence 
is already established at the lower dose and is smoothly maintained 
after the further order-of-magnitude change in dose. The lump at 10 
volts appears to be significant. If it is real, its explanation will require a 
clever refinement in the model. 

We do indeed observe the anticipated increase in device degradation 
at higher bias. The consistency of the shape of the voltage dependence 
at doses greater than 10‘ rads strongly suggests that the dependence 
does not arise only from a difference in the numbers of available surface 
ions, but also from variation in the arrangement of ions and the influence 
of this arrangement on the current characteristics of the channel. The 
collection of charge does not seem to be the dominant variable over the 
voltage range examined. 

Another indication that junction bias affects the arrangement of 
charge on a device surface is shown in Fig. 13. The characteristics ob- 
served after a total dose of about 10’ rads and at a high dose rate may 
show considerable structure, depending on the rate at which the charac- 
teristic is swept out. I'urthermore, the differences between the charac- 
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Fig. 11. — Dependence of degradation in I¢go on collector bias. 


teristics with slow and fast sweep show that rearrangements are not 
instantaneous. Note that here also a peak in Igo occurs in the vicinity 
of 10 volts. 


5.4 Influence of Can Potential 


It has been suggested that the electric field existing between the can 
and the transistor base by virtue of the bias on the collector junction 
may alter the collection of ions at the semiconductor surface. This 
possibility has been examined using diffused silicon diodes, encapsulated 
in the same gaseous atmosphere as the diffused silicon transistors. De- 
vices with the p-region common to the can and others with the n-region 
common were studied. There was no substantial difference in sensitivity 
to ionization in these cases. On the other hand it was found that by 
alteration of the device processing, devices of either polarity were in- 
sensitive. It is clear that in this case the chemical surface condition is 
of more vital importance than the can-to-diode bias. 
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Tig. 12. — Median behavior of Zczo with collector bias at two integrated doses. 


No silicon transistor with reversible polarity to the can has been 
available, but a germanium transistor with all leads insulated from the 
can has been irradiated with reversal of the sign of can bias. The results 
are shown in Fig. 14. The device was irradiated at 8 X 10° rads/hr for 
one minute with the can negative with respect to the transistor base 
(the collector of which was continuously reverse biased) and then for one 
minute with the can positive and so on for longer times as the irradiation 
proceeded. The leakage current is very clearly much more sensitive to 
positive potential on the can than to negative, a result consistent with 
the original picture of the role of positive ions in a channel on the transis- 
tor base. This experiment has been repeated several times with some 
lack of reproducibility. In some cases can bias makes much less differ- 
ence than as shown in Fig. 14 and in two cases where the whole radiation 
effect was smaller, the dependence on can bias polarity was reversed. 

There seems to be no question that the can-to-device potential can be 
important, but clear-cut evidence that positive ions rather than electrons 
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Fig. 18. — Structure and drift in the reverse characteristic of a heavily irradi- 
ated transistor under radiation. 


(or negative ions) are always the important particle in transistor surface 
ionization effects has not been demonstrated. 


5.5 Recovery 


Recovery of the surface effects after an exposure to radiation has 
already been mentioned in connection with Figs. 2 and 8. The recovery 
is not exponential. In many cases, recovery under bias is well repre- 
sented by a straight line on a log Iso vs log ¢ plot, although the points 
at less than 0.5 min tend to fall below such a curve. The recovery repre- 
sents a loss of charge at the surface, perhaps by neutralization from the 
interior in the absence of a radiation flux, but the memory effects suggest 
the active species do not actually leave the surface. 

The comparison between recovery with and without bias is shown in 
Fig. 15. A transistor given a dose of 4.6 10° rads at 15 volts is shown, 
first recovering at this same bias and then recovering with the collector 
circuit open except momentarily for measurement. The decay is sub- 
stantially accelerated in the zero-bias condition. This is consistent with 
the idea that the ions are bound in place by the field but in its absence 
are free to diffuse on the surface and will no longer be concentrated to 
produce an inversion layer It has also been observed that with reestab- 
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Fig. 14. — The influence of can-to-semiconductor potential on degradation in 
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lishment of the bias the current rises before continuing its decay, sug- 
gesting that the ions are subject to recapture by the field. There is a 
loss in this process however, that may arise from the enhanced neutrali- 
zation of the ions by electrons from the inversion layer when they are 
more numerous in the absence of reverse bias. 

Fig. 16 illustrates the effect of radiation on recovery. The upper curve 
was taken first and is a normal recovery with bias. The lower curve 
shows recovery from the same starting value in the 8 X 105 rads/hr 
gamma. field without bias. The device is momentarily removed from 
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Fig. 15. — Decay of the radiation effect with and without bias. 
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Fig. 16. — Enhancement of the rate of recovery of I¢zo with radiation but 
without bias. 


radiation for measurement at 15 volts reverse bias. The recovery is 
substantially enhanced by the ionization in the absence of bias. This 
seems to be explained through neutralization of the charges on the semi- 
conductor surface by the newly formed ions. Under bias the gas ions are 
directed by the junction field to strike the surface in places where they 
add to an inversion layer. In the absence of a field from applied bias, 
these ions tend to go wherever they can reduce fields produced by surface 
charges. 


VI. CHARACTERIZATION OF EFFECTS WITH SIGNIFICANT NUMBERS OF DE- 
VICES 


With the above background concerning the surface effects of ionizing 
radiation, and with recognition of the limited knowledge of the range of 
exact surface conditions existing in semiconductor devices, further tests 
were planned to establish the extent of the ionization effect in larger 
samples of devices. Jor the purpose of such a test program, a transistor 
type was selected which showed particular sensitivity to the combined 
effects of electrical bias and ionizing radiation. These transistors nor- 
mally have quite low collector reverse currents, Iczo (in the order of 
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10-*° amperes). They are therefore good subjects for study at relatively 
low gamma dose levels since small increases in J¢go can readily be rec- 
ognized. A large number of devices were available with sufficient power 
aging to indicate stability of characteristics, so that changes could 
clearly be attributed to the radiation exposure. 

Although the measurements of gain of these transistors under gamma 
radiation indicate some change with dose, the changes are relatively 
smaller than those in the junction reverse current and are hence less 
subject to recognition of a specific pattern of change. The data presented 
here, therefore, relate specifically to the changes in Iczo . 


6.1 Typical Pattern of Degradation with Dose 


In order to determine the proper conditions for large scale evaluation, 
and to have some estimate of changes to be expected, it is desirable to 
establish the pattern of change in characteristics. From the data pre- 
sented in Section V on the various factors affecting degradation, there 
develops a typical pattern of degradation, at least for the type of diffused 
silicon transistor used, as shown in Fig. 17. In this log-log plot of 
Iczo vs integrated dose, four recognizable regions are indicated, an 
initial region of stability followed by three regions of change. The bound- 
aries of the regions as indicated in Tig. 17 depict only representative 
variations and do not indicate actual limits of the device type. 
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Fig. 17. — The typical pattern of Iczo degradation under radiation for a par- 
ticular type of diffused silicon transistor. 
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As discussed in Section 5.3, the slopes of the degrading I¢xo’s become 
quite uniform among all units of this type beyond a total dose of about 
10* rad, or beyond the variable Region 2. In Region 3 all units show 
slopes of approximately 3, and at about 10° rads there is a fairly distinct 
increase in slope which distinguishes Region 4. Although these two 
regions are recognizable, they are of lesser importance in judging com- 
parative usefulness of the devices in radiation than are the region of 
stability and the slope in Region 2. 

The region of stability is indicated even in the early tests (see Fig. 3) 
and continues to show up, even at the lowest initial current levels, 
when the dose rate is low enough to allow measurements at sufficiently 
low total dose. The extent of Region 1, or the dose value at which an 
individual unit starts to degrade, is quite variable among units of a type, 
as well as between types, and is a significant point for further considera- 
tion of devices. 

The transition into Region 2, where the units show a typical linear 
increase in log Iczo with log of dose, is quite distinct in all observed 
cases. For this reason, data which show only the degrading slope in 
Region 2 (perhaps because a high dose rate causes the first measurement 
to be beyond Region 1) can be extrapolated with some confidence back 
to the preradiation value to achieve an estimate of the transition dose. 
Major features of Region 2 are extreme variations between individual 
devices of a type and considerable dependence on collector voltage. This 
dependence is shown in Fig. 11 by the variations in I'¢zo values achieved 
at about 10‘ rad, apparently about the end of Region 2 for this type. 

In general, there appears to be a consistency between Regions 1 and 
2 in that those units which show the greatest change in Region 2 also 
tend to indicate the shortest period of stability in Region 1; and con- 
versely, those showing smaller changes in Region 2 appear to have the 
longer period of initial stability. Within the range of actual initial values 
(from 8 X 10-" to 2 X 10-* amperes) however, there appears to be no 
correlation between the initial value and the subsequent severity of 
degradation. 


6.2 Distributions of Dose at Initiation of Degradation 


With the recognition of the typical response of this transistor type to 
radiation, a question can first be asked regarding the variability in the 
extent of the region of stability, Region 1. 

Fig. 18 shows several plots of the distribution, on a normal probability 
scale, of the integrated dose at which degradation is initiated (the end of 
Region 1). The three lowest curves show the results obtained from ex- 
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Fig. 18. — Distribution of the radiation dose at initiation of degradation. 


posure to gamma radiation at 8.5 X 105 rads/hr., 10 rads/hr., and 1 
rad/hr., all with a collector bias of 15 volts and with no connection to 
the emitters. There appears to be only about a factor of 2 difference 
between the 8.5 X 10° rad/hr. and the 10 rad/hr. results. This evidence 
of reciprocity between the 8.5 X 10° rad/hr. dose and the 10 rad/hr. dose 
at the end of Region 1 is in contrast to the lack of reciprocity indicated 
in Fig. 9 which compares reverse currents in Region 2. On the other 
hand, the difference between the 1 rad/hr. and 10 rad/hr. distributions 
seems to indicate a breakdown of reciprocity, at this level, in the dose 
required to initiate degradation. 
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The two upper curves of Fig. 18 are obtained at 5 rad/hr. and with a 
forward bias on the emitters during radiation in addition to the reverse 
collector bias. One of these curves is on the same transistor type 
(F-54278) as that of the lower distributions. The other distribution at 
5 rad/hr. is on another transistor type (F-54647) which is basically the 
same but with a somewhat higher distribution of initial current gain. 
It would appear that the application of forward emitter bias causes an 
increase of roughly two orders of magnitude in the dose required for the 
onset of Iczo degradation. This could result from the neutralization of 
positive charges on the surface of the base region by electrons injected 
into the base by the forward biased emitter. 


6.3 Distributions of Icro Increase 


Fig. 19 shows the distribution of the increase in collector reverse cur- 
rent resulting after 1.4 X 10‘ rads, which dose is past the completion of 
the variable Region 2. Here the Iczo , on a log scale, is plotted on the 
normal probability scale. These data are on the same types of units 
represented in Fig. 18 and are obtained from radiation with the emitters 
forward biased. As in Fig. 18, the higher-gain type evidences a lesser 
degradation, having a longer region of stability as well as a smaller 
increase in Iczo during Region 2. The excellent match to a log-normal 
distribution (except at very low currents where the data reflects the 
inaccuracy of measuring the difference between two nearly equal num- 
bers) lends confidence to the interpretation that the break at the upper 
end of the distribution is caused by the transition into Region 3 of those 
units which have changed the most in the variable Region 2. 

The broad distributions of Figs. 18 and 19 point up the necessity 
for relatively large-scale experiments in order to make valid comparisons 
between different test conditions or device types. 


6.4 Tests of Other Semiconductor Devices 


Similar kinds of distributions and general responses to the various 
conditions affecting degradation under radiation have been found in 
other codes of diffused silicon transistors and silicon diodes which have 
been encapsulated with a gas filling. 

Fig. 20 shows, for example, the distribution of reverse current increase 
of a 4-watt diffused silicon diode after 3.7 X 104 rads at 7.4 X 105 
rads/hr. The broken shape of the distribution may simply indicate that 
the devices did not come from a single product run. The lack of the 
distinct break in the distribution at high currents suggests that the 
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Fig. 19. — Distribution of current degradation after a dose of 1.4 X 104 rads 
at 8.5 & 10° rads/hr. 


transition to Region 3 may not exist, or exists at a different dose. It 
would appear that each type must be examined in detail for a good 
understanding, even empirical, of its performance under radiation. 

One interesting variation to the kind of response indicated above is 
that of a type of silicon alloy transistor which contains a silicone grease. 
In this case a relatively minor degradation is observed during short 
periods at a high gamma dose rate, but more severe degradation occurs 
during the period subsequent to completion of the exposure. This is in- 
dicated in Fig. 21, a plot of collector reverse current vs time after re- 
peated exposures on an individual device. The device was exposed for 
several consecutive periods at 1.4 X 10° rad/hr., with current measure- 
ments subsequent to each exposure. The procedure was then continued 
at 8.5 X 10° rad/hr. Here each curve represents the Jcz0 measurements 
subsequent to radiation for the total time indicated on the curve. Icgo 
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Fig. 20. — Distribution of degradation of current in a type of diffused silicon 
diode at 3.7 X 104 rads. 


increases subsequent to the exposure, but may actually be decreased 
during the next exposure. After sufficient exposure, the devices evidence 
a saturation of the subsequent [cso , and also a saturation of the effect 
with continued exposure. 

This saturation effect is confirmed at the 5 rad/hr. dose rate in Fig. 22, 
showing the Iczo curves vs integrated dose for two typical transistors 
of this type. Both show an eventual saturation of the curve and one 
shows an ultimate reduction which appears also to be typical of this 
type of transistor. 

Several types of diffused germanium transistors have also been ex- 
amined, with the general result that increases in Icxo are relatively 
small (of the order of a factor of 3 or 4), and relatively consistent, until 
the dose reaches about 10° rads, where more drastic increases may be 
expected. The dependence on V cg was relatively insignificant. The effect 
of processing is evident in these types, however, with one group showing 
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Fig. 21. — The pattern of collector current response in a grease-covered silicon 
transistor. 


marked changes at about 10! rads, and another group having essentially 
no changes at 10’ rads. 


VII. TESTING AND SELECTION FOR TELSTAR DEVICES 


The process of selecting the semiconductor components for use in the 
Telstar experimental communication satellite consisted of: 


1. 


2. 
3. 


4, 


5) 


Qualification of design for reliability and performance as required 
in each application, 

Fabrication, 

Screening and preaging to assure satisfactory operation in the 
system environment, 

Life testing, 

Selection of the most stable devices. 


At the time of recognition of surface effects from ionizing radiation 
(about October 1961), much of the life testing was already in process. 
Since the Telstar satellite was to orbit through the Van Allen radiation 
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Fig. 22. — The response of a grease-covered silicon transistor under prolonged 
radiation at 5 rads/hr. 


belt, the addition of radiation as an environmental factor was essential. 
Steps were therefore taken to determine the qualification of all types in 
the program and to study the screening and selection techniques that 
would be useful. 

Fig. 23 shows a plot of the estimated ionizing radiation intensity for 
components inside a satellite, as a function of the shielding of the com- 
ponents from the external environment. These curves refer to the ioni- 
zation produced by protons only and were developed from existing esti- 
mates of the Van Allen belt particle fluxes. The upper and lower limit 
curves reflect the uncertainties in this estimate. For shielding thicknesses 
of less than 0.1 inch the electron contribution to ionization should be 
considered as well, but for Telstar, typical components are shielded by 
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Fig. 23. — The maximum radiation intensity anticipated in space due to high- 
energy protons, as a function of the thickness of aluminum shielding. 


the equivalent of 0.5 to 1.0 inch of aluminum and only high-energy 
protons are significant. The curves of Fig. 23 indicate that at this 
typical shielding, the radiation intensity would be a maximum of 10-20 
rads/hr in the heart of the Van Allen belt. This results in a maximum 
average of 3-5 rads/hr over an entire orbit, considering Telstar’s approxi- 
mate 25 per cent effective exposure to Van Allen radiation. 

With this estimate of the maximum radiation environment, all device 
types were given a gamma exposure at 8.5 X 105 rads/hr for one minute 
(the equivalent of at least three months in orbit), followed by at least 
one week at 3 rads/hr. Any device type showing evidence of change in 
either condition was replaced or subjected to individual selection or 
screening. Device types showing no change were considered satisfactory. 
All of these devices were subject to a 15-v reverse bias on the collector 
and an emitter current corresponding to the application. 

The diffused silicon transistors I'-54273 and F-54647 were used for 
experimental studies of screening and selection procedures. It is noted 
that the one-minute dose of 1.4 X 10‘ rads (at 8.5 X 10° rads/hr) is not 


124 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


much beyond the variable Region 2 of Fig. 17 and should be very effec- 
tive in providing a comparison between individuals. The study of the 
effectiveness of such a dose as a screening procedure was carried out 
through the following program: 

1. Preradiation of a number of these transistors to screening expo- 

sures between 10 seconds and 6 minutes at 8.5 X 10° rad/hr and 

2. Subsequent exposure at 5 rad/hr. 

This permits an evaluation of the effect of the low dose rate after a 
screening dose. 

Fig. 24 shows the collector reverse current measurements of two units 
from this program, these being generally typical of the results of all of 
the devices. In this case, the measuring equipment was limited in sen- 
sitivity to about 10-® amperes and the initial values prior to radiation 
are shown at this value although they may have been somewhat lower. 
Unit No. 762 is representative of those units which suffered a relatively 
minor increase in reverse current during the preradiation dose. During 
the subsequent radiation at 5 rad/hr, this unit returned quite rapidly to 
its original value, remaining there until the dose became somewhat 




















Togo IN AMPERES 





























AFTER 7x103R DOSE AT 5R PER HR IN RADS 
AT 8.5x105 R/HR 


Fig. 24. — The typical pattern of change in I¢go in low-level radiation after a 
screening dose. 
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greater than 103 rads, at which time it began to degrade in a fashion 
very similar to that expected from earlier tests. Unit No. 722 suffered 
a much greater change during the preradiation and took an appreciably 
longer time recovering toward its initial value. Before it fully recovered 
it reached the point of onset of final degradation and began to change 
quite rapidly. (A third type of response was seen in an occasional unit 
which degraded so severely during the preradiation that the subsequent 
low dose resulted in further increase in current and no recovery was 
evident at all.) 

The influence of the screening dose on the distribution of onset of the 
final [cso degradation (at 5 rad/hr) is shown in Fig. 25 for four similar 
groups of 12 units each. One of these groups was irradiated at 5 rad/hr 
with no initial radiation at the high dose level. The other three groups 
were given an initial radiation at 8.5 X 10° rad/hr for different lengths 
of time to achieve the initial radiation dose indicated in the figure. It 
is seen that all of these distributions are essentially the same, indicating 
that the 5 rad/hr dose will cause a modification of the surface condition 
established by the high initial dose, causing the units to look eventually 
as if they had not received the initial dose. This is evidence of a contra- 
diction to the principle of reciprocity of dose rate and time, in that the 


ZERO 
2.3x 103 
1.4x 104 
8.4x 104 





pe ef LT 
pte 
as Oe 
eas ee ae 


DOSE IN RADS AT INITATION OF Icgq DEGRADATION 
3 
w 


wn 


10 20 30 40 50 60 70 80 90 95 
PER CENT EQUAL TO OR LESS THAN ORDINATE 


Fig. 25. — Distribution in dose for initiation of degradation with and without 
pre-irradiation. 
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additional dose added at the 5 rad/hr rate does not normally cause a 
continuation in the degradation produced in the initial dose. 

Another point of interest in Fig. 24 is the comparison of the Iczo 
value resulting from the preradiation dose with that subsequently occur- 
ring after an equal dose at 5 rad/hr. Inspection of the two curves on 
Fig. 24 reveals that the currents at 7 X 10° rad are approaching those 
resulting from the preradiation dose. I‘ig. 26 shows a plot of the results 
of all the units so treated. The dashed line indicates the one-to-one 
correlation between the two current values, those points to the right of 
the line representing the units which during the 5 rad/hr exposure did 
not develop reverse currents as high as those obtained in the initial 
dose. This line was not extended below about 4 X 10-® amperes because 
this was approaching the limit of sensitivity of equipment which was 
then in use. It is noted that all but two of the significant readings fall 
on the side of the one-to-one correlation line corresponding to larger 
changes after the preradiation dose than after the subsequent low dose 
rate exposure. 

In some cases the lack of correlation is quite significant, the currents 
after the low dose rate exposure remaining below 10-° amperes although 
quite large changes were observed after the screening dose. In these 
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Fig. 26. — Correlation of change in I¢go after a screening dose at 8 X 105 rads/ 
hr with subsequent changes at 5 rads/hr. 
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units, however, the changes in current gain after the low dose rate ex- 
posure were found to be appreciably greater than in those units which 
were more stable in the initial screening dose. Consequently, it is found 
that selecting of those units with Icgo less than 10-8 amperes after the 
screening dose would have eliminated 9 of the 10 units which subse- 
quently degraded to either Iczo values greater than 10-’ amperes or 
gains less than 50 per cent of the initial gain. Of the 12 units good in 
these respects after the low dose rate exposure, only two would have 
been eliminated by selection following the screening dose. It is thus 
shown that selection on the basis of Icso after a screening dose is effec- 
tive for this device type. 

An indication of the effect of this screening on the distribution in 
degradation of [czo under low-level radiation exposure is given in I'ig. 
27. Both distributions are of [cso after exposure, first to 30 seconds at 
8.5 X 105 rad/hr and then to 10‘ rad at 5 rad/hr. One distribution is of 
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Fig. 27. — Distribution of radiation response in screened and unscreened devices. 
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the entire sample of devices and the other is of only those which were 
selected for Ico less than 10~-* amperes after the initial dose. It is seen 
that an appreciable improvement is achieved by the screening procedure. 


VIII. SUMMARY 


It has been found that changes can occur in semiconductor charac- 
teristics, because of surface effects of ionizing radiation in a device under 
electrical bias, at much lower doses than those required to produce 
surface effects in the absence of bias or to produce changes in the bulk of 
the semiconductor. These effects apparently arise from ionization in the 
gas of a device encapsulation and collection of ions on the device surface. 
The gas ions probably serve to produce and activate chemical species on 
the semiconductor which induce surface inversion layers that alter the 
junction characteristics. The chemical condition of the surface previous 
to irradiation is apparently involved in the process in a sensitive way. 
The effects observed depend on junction bias, envelope potential, and 
in many respects on total radiation dose rather than dose rate. The 
major features of these observations and of observations of the char- 
acteristic recovery of the surface effects after radiation under bias can 
be described by the ion-induced inversion layer model of the process. 

In the time since the first direct observation of this effect, a compre- 
hensive study of all types of semiconductor devices has not been at- 
tempted. Observations of several types, however, indicate that different 
types may respond quite differently to radiation, that the response may 
be quite dependent upon processing (and therefore upon production 
periods or batches), and that within a type the response may be quite 
variable between individual units. 

One type of transistor was used for much of the experimental work of 
defining the radiation effect, and the cso response was studied because 
the large changes facilitated comparative measurements. It is recognized 
that current gain and other surface-dependent characteristics can also 
be affected by radiation, and they should be observed in any evaluation 
of types, if critical in specific applications. It is hoped that the studies 
presented here will be a guide in formulating such evaluations. 

It has also been found that at least some devices can be screened for 
sensitivity to radiation by means of a short-time, high-level dose, with 
correlation to subsequent low-level exposure results. Here, too, it is 
recognized that variations exist among types, and the usefulness of such 
screening operations should be evaluated for each type of interest. 

It is hoped that further work will serve to provide more definitive 
results relating the radiation effects to a physical model and also to the 
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surface conditions or processes contributing to the effect. These prelimi- 
nary results can at least serve to expose the problem and to suggest the 
lines of further study and action. 
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Development of Solderless Wire Connector 
for Splicing Multipair Cable 


By H. J. GRAFF, J. M. PEACOCK, and J. J. ZALMANS 
(Manuscript received September 17, 1962) 


This paper describes the development and laboratory testing of a new 
high-reliability solderless connector for splicing cable conductors. It in- 
cludes discussion of the physical parameters which influence the perform- 
ance of electrical contacts in general. The experimental and analytical 
techniques which evolved as part of this project permit important reductions 
in the amount of experimental data required to make reliability and aging 
predictions, and should be useful in other problems dealing with the ap- 
pratisal of electrical contacts. 


I. INTRODUCTION 


The development of new connectors and new techniques for making 
wire connections has been actively pursued at Bell Telephone Labora- 
tories for some time. Current estimates indicate that well over two 
billion wire connections are made each year in the Bell System. These 
connections obviously represent a large investment, and the problems 
of cost and reliability are becoming increasingly difficult as modern 
circuits grow in complexity. Performance testing of a connector to 
establish its reliability is a challenging task. 


1.1 B Wire Connector 


The purpose of this paper is to describe the development and per- 
formance appraisal of the B Wire Connector recently released for splic- 
ing outside plant multipair cable conductors. As illustrated in Fig. 1, 
this connector consists of a thin springy liner with sharp tangs on the 
internal surface, and a brass outer shell encased in a plastic jacket. In 
use, the conductors to be joined are inserted without removing insula- 
tion, and the connector is pressed with a pneumatic tool. The sharp 
tangs penetrate the insulation and establish contact with the conductors 
as illustrated in Fig. 2. The individual tangs act as springs, storing en- 
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Fig. 1 — B Wire Connector. 


ergy to maintain intimate contact with the conductors over the long 
service life of the device. Electrically, the contact is equivalent to a 
soldered joint, and therefore is suitable for use in any type circuit where 
multipair cable may be used. Demand for the B Wire Connector is 
expected to reach about 250 million per year with resulting annual 
savings of several million dollars. 


1.2 Reliability and Aging Requirements 


The application for which the B Wire Connector is intended demands 
the utmost in reliability. The connector, of course, is a series element; 
a typical cable circuit has a splice every 500 feet or so, often as many as 
20 splices in a line. The connections are made under field conditions, 
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Fig. 2— B Wire Connector — cross section after pressing on wires. 
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and are expected to match the service life of the cable, frequently well 
over 40 years. The splices are enclosed in splice cases or lead sleeves, 
many of which are hermetically sealed. 

When reduced to quantitative values, the reliability requirements 
are very high indeed. Consider a one-mile run of 2424-pair cable, with 
ten splices. 48,480 connections are needed, and if 0.1 per cent or 48 of 
these are bad initially, or become bad in time, 48 or 2 per cent of the 
circuits would be lost, neglecting the possibility that two bad connections 
could occur in the same circuit. Similarly, at a failure rate of 0.01 per 
cent or 1 in 10,000, five circuits would not be available. To assure satis- 
factory field acceptance, and to allow margin for other sources of failure, 
target values for the inherent reliability of the connector were estab- 
lished at 1 in 100,000 after 40 years of aging in the outside plant environ- 
ment. 


1.3 Performance Appraisal Methods 


The appraisal of a device for reliability and performance requirements 
such as these is always a very difficult and challenging task. As part of 
this project, a number of advances have been made in the techniques 
usually used for contact appraisals so as to reduce the amount of experi- 
mental data required and to improve the precision of the aging predic- 
tions. Some of these techniques are applicable to development work on 
other types of connectors, and will be discussed in this paper. First it is 
necessary to describe the mechanisms which influence an electrical con- 
tact. 


II. CONTACT MECHANISMS 


Wire connectors in use in the Bell System today take many forms, 
ranging from the time-proven binding post, through the solderless 
wrapped devices, to the newer insulation-piercing, quick-connect ter- 
minals. Although radically different in many of their features, these de- 
vices have much in common — the basic mechanisms that influence the 
behavior of the contact are the same. In all of them the circuit is estab- 
lished by clamping together two metallic elements, the wire and the 
terminal. Circuit continuity is established through the areas of metal- 
to-metal contact. This contact can be considered to have two compo- 
nents of ‘‘contact resistance”: (1) constriction resistance due to the con- 
vergence of current paths within the parent metal to the relatively 
small points of contact, and (2) interfacial resistance due to the possible 
presence of resistive surface films. 
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A number of variables influence the actual value of the contact re- 
sistance in a given connection. Some of these can be controlled quite 
well through proper selection of materials and platings, geometry, etc. 
Others tend to be more random in nature: for example, the fine-grain 
details of surface roughness or oxide coverage which would influence 
the actual areas in contact at a given loading. Due to the random nature 
of some of these mechanisms, a distribution of contact resistance may 
be expected. The spread or range of resistance is a measure of the con- 
trol or reproducibility of a connector, and thus tends to indicate its 
reliability. To hold the degree of variability to a minimum, it is neces- 
sary to understand the various environmental and aging mechanisms 
applying to a particular design. 

Having established an initial contact of a given quality, one must be 
alert to possible degradation of the contact with time. This is likely to 
happen as a result of corrosion, stress relaxation, ‘“‘breathing”’ of oxidiz- 
ing atmosphere due to contraction and expansion of connecting parts 
consisting of different materials, or mechanical disturbance. A good 
degree of mechanical rigidity, with adequate contact force and energy 
storage to maintain these forces over the years, are necessary design 
features to achieve long service objectives. With inadequate provision 
for these requirements, high resistance and noisy connections are in- 
evitable, and under dry circuit conditions (low-voltage, low-current 
signals) barrier films can easily cause complete circuit failure. The re- 
quirements for sufficient areas of intimate “gas tight”? contact and for 
energy storage have been well known for some time and were studied in 
some detail during the development of the solderless wrapped connec- 
tion.} 


III. MEASURING TECHNIQUES 


Measurement of the resistance distribution and aging trends sug- 
gested above poses quite a challenge. For a “good” connector design, 
i.e., one with low contact resistance, the variation of contact resistance 
from sample to sample may be in the order of a few microhms. Resist- 
ance measuring techniques commonly used in the past for the appraisal 
of contact devices have been capable of resolving joint resistances only 
to the low milliohm range. With these techniques it is possible to single 
out “high-resistance”’ joints but not to study the statistical distribution 
of “good” joints. The sample sizes needed to establish reliability rates 
and to draw conclusions about designs and design variables with these 
relatively insensitive measuring techniques are prohibitively large. On 
the other hand, with experimental data of sufficient precision, the be- 
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havior of the many variables can be studied in much greater detail, and 
analytical techniques can be brought to bear so as to reduce drastically 
the amount of experimental data required. For effective analysis by 
statistical techniques with reasonable sample sizes, the data must be 
collected with instrumentation sensitive enough to permit the quantita- 
tive measurement of differences in resistance levels from joint to joint. 
Also, for analysis of the aging mechanisms in types of connectors where 
this is of concern, changes in contact resistance with time must be meas- 
urable at short time intervals (several days or a week). If the above re- 
quirements are met, statistical study of contact resistance distributions 
is possible — both initially and after any desired interval of normal or 
accelerated aging — and reliability rates can be estimated with small 
sample sizes. 

Techniques for measuring fractions of a microhm of resistance have 
been known for years. The problem in measuring contact resistance on 
actual devices to these levels of accuracy lies principally in separating 
the variables in such a way that the experimental errors do not swamp 
out the resistance to be measured. A typical device, for example a con- 
nector joining two wires, may be thought of as a two-terminal black 
box, on which only external measurements may be made. The various 
resistive elements within the box include the internal resistances of the 
lead wire and connector in addition to the contact resistances between 
wires and connector. The normal tolerances on the internal resistances, 
due to uncertainties in wire diameter, length, temperature and resistiv- 
ity, are generally larger than the true contact resistance of connectors 
of high enough quality to be of interest for long-term service on dry 
circuits. These uncertainties must be either eliminated or accurately 
compensated if the desired measurements are to be obtained. 

A very simple circuit has been developed to permit resistance meas- 
urements of the kind contemplated here. This circuit, as shown in Fig. 
3, is a modified Kelvin bridge in which the unknown resistance of the 
contact is compared to a “standard,”’ which is an accurately controlled 
length of one of the lead wires. 

The circuit has a number of features which reduce experimental 
errors: 

1. Variations in lead wire diameter and resistivity are compensated 
exactly by proportionate changes in resistance of the standard. 

2. Resistance measurements are insensitive to temperature changes, 
at least to the extent that the unknown has the same temperature co- 
efficient as the standard. 

3. The effects of thermal potentials can be balanced out in de bridge 
measurements by using normal and reverse polarity of test current. 
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Fig. 8 — Modified Kelvin bridge — setup for loop resistance measurement. 


4. The effect of variable resistances due to the connection of the sam- 
ple to the test set (frequently a source of appreciable error) may be made 
arbitrarily small by putting them in series with the high impedance 
arms of the bridge. 

The loop resistance measurement, which includes the internal re- 
sistance of the lead wires and connector in addition to the contact re- 
sistance, lig. 3, can be used effectively to determine variations in con- 
tact resistance between connections, and to measure changes in contact 
resistance with time. To keep the experimental errors to a minimum, 
these measurements require very close control of the reference points 
from which resistance measurements are made. With the measuring 
equipment used, mechanical tolerances are such that the average varia- 
tion in loop resistance due to these factors is about +4 microhms (for 
19-gauge conductors and 100 milliamperes test current). 

In some cases it is possible to make a more cumbersome measurement 
than the simple “Joop”? measurement shown in Iig. 3 so as to separate 
the various “black box’ resistance elements from the contact resistance 
itself. Modifications in the test set connections as shown in I'ig. 4 have 
been made for this purpose. The unknown bridge arm in Fig. 4(a) in- 
cludes the resistance of one of the lead wires in to the point of contact 
and the contact resistance between the two wires. In Fig. 4(b) the un- 
known contains merely the resistance of the lead wire to the “centroid” 
of the contacts. The difference between these two readings gives us the 
contact resistance by itself. These multiple measurements are quite 
accurate for simple point-contact configurations. or more complicated 
structures having distributed contacts, the observed resistance must be 
interpreted in terms of the distribution, somewhat analogous to the 
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Fig. 4 — Modified Kelvin bridge. Circuit (a) is used to measure the resistance 
of the top lead wire in to the point of contact, plus the resistance between the 
two wires. Circuit (b) measures merely the lead wire resistance. The contact 
resistance is calculated from these two measurements. 


input observations on the distributed parameters of a transmission line. 
In many cases the primary interest is in changes that occur with time 
rather than in the absolute value of contact resistance, and the more 
straightforward loop resistance measurement is sufficient for these pur- 
poses. 


IV. AUTOMATED MEASURING APPARATUS 


The need for automated measurements and streamlined data handling 
was recognized in the early development stages of the B Wire Connector. 
The manual bridge described in Section III was automated to read and 
record joint resistances of groups of up to twenty joints automatically. 
This equipment is shown in Fig. 5. A block diagram of the automated 
bridge giving major components is shown in Fig. 6. The basic circuit is 
the same as for the manual bridge. The self-balancing feature is obtained 
by using the bridge unbalance voltage as the input to the servo-system 
of a self-balancing potentiometer. The shaft position of the potentiom- 
eter, directly related to the resistance measured, is encoded and con- 
verted into digital form for direct recording on punched cards or for 
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Fig. 5 — Automated, self-balancing bridge for repetitive joint resistance meas- 
urements. 


visual readout. A timing unit controls the readout, current reversal 
and the scanner switch by which connections to individual joints 
mounted on a board are made. 

Fig. 7 shows a mounting board with 16 joints. The board also includes 
four looped conductors, which are used as controls to detect possible 
errors due to misalignment or shifting of the mounting board relative 
to knife-edge contacts on the contact board. The use of controls is arbi- 
trary, of course. In this instance, it will be noted that the length of the 
loop was random. This mounting arrangement provides good control of 
reference points and saves considerable time on repetitive measurements. 

As noted in Fig. 3, the loop resistance of the joint is compared to the 
resistance of a “standard” length of wire, in this case, one inch. It has 
proven convenient to read out the resistance in thousandths of an inch 
of conductor, measured from the entrance end of the connector, rather 
than in ohms. Tor the purpose of this paper, the resistance values 
plotted have been converted to milliohms. It should be noted that the 
wire resistance varies with the wire gauge used, and this must be taken 
into account in comparing one plot with another. The precision or re- 
peatability of the system is well within ten microhms, the actual value 
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Fig. 6 — Block diagram of automated bridge. 


depending on the gauge of the wire being tested and the test current 
used. 


V. EVALUATION OF B WIRE CONNECTOR 


The objectives in the design of the new connector were to improve 
the electrical performance and reliability of the previous standard 
“pigtail” cable splice, and to achieve economies through time savings 
by simplifying the splicing methods. The program which led to the final 
design of the B Wire Connector shown in Fig. 1 will be described in this 
section. 


5.1 Evolution of Basic Design 


After considering a number of wire joining methods, it appeared that 
a device which could be pressed onto the ends of unskinned conductors 
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Fig. 7 — B Wire Connectors on mounting board used in automated bridge. 


offered the most attractive solution. For this particular application, 
initial studies showed that such a device should consist of two metallic 
members: an inner member of relatively hard material with some form 
of penetrators to pierce the insulation, and a soft, ductile outer member 
to hold the penetrators in contact with the metallic conductors after the 
pressing force is removed, and thus minimize the effects of springback. 
With this principle in mind, two alternative designs were proposed. 
These two designs used the same soft outer brass sleeve, but differed in 
the design of the inner member. The liner in connector X consisted of a 
helix formed from hard drawn rectangular brass wire that had been 
twisted on its own axis. Connector Y had a liner of thin spring bronze 
pierced with numerous holes so as to form a multiplicity of tang-like 
penetrators on the inner surface. At this early stage of the design effort, 
a number of questions arose, the principal one being, ‘‘Which design 
approach offers the most promise?” 

To appraise the performance of connectors X and Y, loop resistance 
measurements were obtained on small groups of comparably made 
joints of both types. Measurements were repeated periodically, during 
which time the connections aged naturally at laboratory ambient. Mi- 
nute differences in joint resistance resulting from statistical variations in 
contact characteristics and the aging phenomena were readily discerni- 
ble with these measurements. Mean connector resistances and standard 
deviations were plotted vs log time and are shown in Fig. 8. 

A study of the aging characteristics, as illustrated by the mean re- 
sistance vs time, shows a significant difference between the two con- 
nectors: connector Y has a lower aging rate. The plot also gives some 
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Fig. 8 — Comparative resistance plots for X and Y connectors showing per- 
formance under natural aging conditions. 


insight as to the fundamental mechanisms affecting the performance of 
connector Y. The linear relationship of resistance and log time indicates 
the aging phenomenon is following some exponential law. Further 
analyses of these and other data suggest that the principal component 
of aging results from stress relaxation of the soft copper conductor at 
the various points of contact. Stress relaxation has been studied exten- 
sively, is predictable, and can be compensated by the provision of in- 
creased energy storage in the connector. 

The standard deviation plots are also of interest since they provide 
an indication of the reliability that may be expected from larger samples. 
Any significant increase in the deviation with aging is an indication of 
instability and poor reliability. In view of these considerations, it can 
be concluded from Fig. 8 that connector Y has both a potentially longer 
service life and a higher reliability rate than connector X. Fortunately, 
connector Y is also more manufacturable. 

It is of interest to note that the initial performance of the two con- 
nector designs is not radically different. Based on the initial results, 
either design might be expected to give good results. In less than one 
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month of natural aging, however, significant differences in the aging 
rates of the two designs became apparent, in spite of the fact that only 
a handful of samples were under study. As verified by subsequent data, 
a valid choice could have been made between the connectors at this 
time. It is apparent that the ability of the test method to resolve and 
detect small changes in resistance is of considerable assistance in pre- 
dicting the long-time performance of connector designs with a minimum 
number of test samples. 

If the resolving power of the test method were degraded to -+1 milli- 
ohms, as an example, a significant difference between the two designs 
could be detected eventually, after the lapse of a much longer aging 
time. Also, it is possible to improve the accuracy of the predictions of the 
standard deviation, in spite of the poor accuracy of the test method, by 
testing many more samples, and by applying the appropriate statistical 
principles. Neither alternative offers an attractive substitute for good 
resolving power in the test method. 

Still another alternative, in case the resolving power of the test method 
is poor, is to resort to accelerated aging. A good example of an acceler- 
ated test was used in developing the solderless wrapped connector.’ 
Here, the connections were exposed to 105°C for 150 days, with inter- 
mittent mechanical disturbance. One criterion for judging the perform- 
ance of the connector is the magnitude of resistance instability resulting 
from the mechanical disturbance. The measurements are repeated peri- 
odically as the stress relaxation proceeds. Due to the mechanical stabil- 
ity and rigidity of connectors X and Y, mechanical disturbance does 
not produce significant resistance fluctuations, certainly not any that 
could be detected by measurements of 1.0- or even 0.5-milliohm sensi- 
tivity. 


5.2 Development of Final Design 


To complete the design of connector Y, now known as the B Wire 
Connector, a number of experiments were conducted in which the 
pertinent design parameters were varied to determine their effects on the 
performance. An interesting example is a Latin Square experiment per- 
formed to study the geometric variables involved in the design of the 
perforated liner. In this experiment, three levels of pierced hole size and 
spacing were studied for each of three different material thicknesses. 
The results of this experiment are shown in Table I, which is a tabula- 
tion of the resistance change after 1400 hours of aging at laboratory 
ambient for each of the design combinations. An analysis of variance 
indicates that the 0.005 inch thickness is inferior from the standpoint 
of contact resistance to either the 0.006- or 0.007-inch thickness at the 
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Liner Thickness in Inches 
Hole Size in Inches: Row Sum: 
0.005 0.006 0.007 

0.025 0.093* 0.001 0.012 
0.047 0.004 0.001 
0.026 —0.005 0.001 
0.020 0.019 0.000 

0.186 0.019 0.014 0.219 
0.030 0.045 0.024 0.009 
0.024 0.019 0.018 
0.076 0.019 0.010 
0.080 0.002 0.003 

0.225 0.064 0.040 0.329 
0.035 0.042 0.003 0.007 
0.025 0.012 0.005 
0.019 0.009 0.007 
0.034 —0.001 —0.002 

0.120 0.023 0.017 0.160 

Column Sum: 0.531 0.106 0.071 0.708 


* Resistance increase in equivalent inches of 26-gauge conductor after 1400 
hours of aging. This conductor has a nominal resistance of 3.48 microhms/mil. 


1 per cent level of significance. The analysis further indicates that the 
variations occurring between pierced hole sizes are probably due to 
chance causes. Although differences in performance could be noted with 
the measuring precision available, this experiment, and others of a 
similar nature designed to study other variables such as material hard- 
ness, tang shape, etc., showed quite clearly that the connector perform- 
ance was not critically dependent on the minor details of the liner 
fabrication, and that reasonable manufacturing tolerances were not 
likely to upset the reliability estimates. 

In other experiments, accelerated aging tests were used to good ad- 
vantage to magnify the effects being studied. One such experiment was 
used to demonstrate the advantage of tin plating the liner. In this ex- 
periment, a group of connections was subjected to cyclically varying 
humidity. The pulp insulation used on cable conductors is quite sensi- 
tive to changes in humidity. Pulp is trapped under compression in the 
connector, and humidity changes cause an alternate shrinking and 
swelling of the pulp and consequently a mechanical working of the 
joint. This test has proven very effective in evaluating the stability of 
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the various contacts within the connector. Fig. 9 shows graphically the 
behavior of two small groups of joints under this accelerated test. One 
group has liners with tin plating; the other group has liners without 
plating. 

Analysis of these results led to a good understanding of the manner in 
which the individual contacts are made. I’ig. 2 is a cross-sectional view 
which shows the liner tangs contacting and digging into the conductor. 
Due to spring-back, immediately after pressing, small voids occur here 
and there between the tang and conductor. The mechanical working of 
the joint by humidity changes causes sufficient movement to dislodge 
the contacts. Tin plating, however, provides a soft, smearable surface 
which is evidently quite effective in filling these voids and preventing 
motion. In addition to better stability, the increased contact area due to 
tin plating reduces the contact resistance quite significantly, and at the 
same time reduces the statistical deviation among similar joints. Also, 
previous work has shown that tin-plated surfaces under sufficient con- 
tact pressures are subject to solid-state diffusion, which, of course, 
would improve the electrical contact as a function of time. 


5.3 Accelerated Testing 


The design and interpretation of accelerated tests, although neces- 
sary in many cases, are always both difficult and controversial. Under 
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Fig. 9 — Humidity cycling test results which show the superior performance of 
tin-plated contact surfaces. 
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all normal environments, a number of aging mechanisms are involved, 
with various degrees of interaction between them. An accelerated test 
should ideally accelerate each of the mechanisms equally and preserve 
the interaction factors; all the pertinent mechanisms must be included 
or accounted for, and any new mechanisms that may emerge in the 
future must be anticipated. The development of a test that will satisfy 
all these requirements is a Utopian goal, and asa practical matter, some 
degree of engineering risk is inevitable in the use of accelerated tests. 

During the development of the B Wire Connector, considerable effort 
was devoted to the identification of the pertinent aging mechanisms 
and to the design of accelerated tests to study their influence. One in- 
teresting approach has been to design a test that is obviously too severe; 
if the connector survives this test, then a good factor of safety is assured, 
and the quantitative calibration of the test becomes unnecessary. The 
humidity cycling test already mentioned is an example. In its proposed 
usage on pulp-insulated conductors, the B Wire Connector will always 
age in a hermetically sealed, low-humidity cable system; and by com- 
parison, the accelerated humidity cycling test is quite severe. The good 
performance of the tin-plated inserts gives us confidence that the con- 
nector has an adequate margin of safety against the effects of changing 
humidity without further calibration. 

In other cases, it has been necessary to attempt a more precise cal- 
ibration. For instance, we are concerned with the aging mechanism re- 
sulting from daily and seasonal temperature fluctuations. These changes 
cause a mechanical working of the contacts brought about by the differ- 
ences in coefficient of thermal expansion of the conductor-connector- 
insulation components. This type of aging could result in contact in- 
stability. In addition, the mechanical working of the contact parts may 
expose the contact areas to the entry of gas and to the formation of 
oxides or insulating films. 

In evaluating the performance of the connector on polyethylene- 
insulated conductors, the temperature change becomes an especially 
important variable. The insulation trapped under compression in the 
connector causes considerable mechanical working of the contacts due 
to its coefficient of thermal expansion, which is much larger than the 
coefficients of the metallic parts of the connector. To determine the 
effect this mechanism has on connector performance, two tempera- 
ture cycling tests have been devised. These tests have the following 
temperature ranges: (1) —320° to +160°F, and (2) —40° to +140°F. 

The first test was chosen in order to subject joints to a thermal shock. 
It consists of heating the joints to +160°F and chilling them rapidly by 
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immersion in liquid nitrogen. This test proved effective in causing an 
increased rate of aging. Fig. 10 shows the accelerated aging curve of 
polyethylene-insulated 19-gauge copper conductor joints as compared 
to pulp-insulated conductor joints. Correlation of this and natural aging 
data helps to estimate the number of cycles representing the expected 
service life. Relationship between number of temperature cycles of 
—320° to +160°F and hours of natural aging for pulp- and polyethylene- 
insulated cable conductors is shown in Fig. 11. These curves indicate 
that 1 year of natural aging is approximated by something in the order 
of 50 cycles, and that the equivalent natural aging time increases as the 
cube of the number of temperature cycles. Assuming this correlation 
exists for future time, extrapolation of these curves shows that 40 years 
natural aging is approximated by some 200 cycles. 

Since the temperature changes in natural environment are more 
gradual, the connectors will not be exposed to thermal shocks of the 
magnitude present in the first temperature cycling test. To represent 
actual temperature extremes in the field, the —40° to +140°F tempera- 
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Fig. 10 — Liquid nitrogen temperature cycling test results showing the accel- 
eration of aging for polyethylene-insulated conductors as compared to pulp. 
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Fig. 11 — Correlation of liquid nitrogen cycles to natural aging. 


ture cycle was chosen. The smaller temperature range has the advantage 
of simulating the natural environment more closely than the thermal 
shock test and, therefore, is being used exclusively for the final evalua- 
tion of plastic-insulated conductor joints. The experimental correlation 
obtained between natural aging and cycles of temperature cycling be- 
tween —40°F to +140°F is shown in Fig. 12. Early in the test fewer 
cycles are required, but the straight-line relationship observed between 
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Fig. 12 — Correlation of —40°F to +140°F temperature cycles to natural aging. 
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liquid nitrogen cycles and natural aging does not hold for the longer 
—40°F to +140°F temperature cycles. The explanation for this devia- 
tion is not known, but it is suggested that presence of nitrogen within 
the connector may exclude oxygen or other corrosive gases from the 
critical contact regions of the joint; also, the considerably longer cycle 
time for the —40°F to +140°F test may have a bearing on this devia- 
tion — about 5 cycles per day against 50 cycles for the liquid nitrogen 
to +160°F test. 

An important factor contributing to aging and not included directly 
in the above tests, is corrosion. Tests where joints are subjected to cor- 
rosive atmospheres at constant temperature do not simulate the ‘“‘breath- 
ing” that occurs between contact surfaces due to uneven expansion and 
contraction of contact members with changing temperature. Calibration 
of room temperature corrosion tests based on corrosion film thickness 
may therefore be misleading, because the penetration of the corrosion 
products between contact surfaces during the short exposure time is 
likely to be insignificant. 

To simulate the “breathing” of corrosive gases between contact sur- 
faces, temperature cycling in the presence of corrosive atmosphere has 
been tried with reasonable success. This test approaches more nearly a 
true accelerated aging test because, in addition to forcing a mechanical 
working of the contacts, it accelerates the corrosion process. Ideally, the 
same degree of acceleration should be used for the two mechanisms. 
Attempts to correlate the aging of corrosion-cycled B Wire Connectors 
and natural aging have not been successful to date. The heavy corrosive 
films that form in a short period of time interfere with the test-set probes 
and also significantly change the resistance of the conductor used as a 
“standard” in the bridge measurement. The test has been useful, how- 
ever, for comparison testing of different types of connectors. 

As indicated above, stress relaxation of the soft copper conductor at 
the various contact points in the B Wire Connector is a possible aging 
mechanism. Fundamental experiments have established that relaxation 
does occur in stressed copper structures and that it is possible to accel- 
erate this process with elevated temperature exposure. For example, 
previous research has shown that 3 hours at 175°C is equivalent to the 
stress relaxation during 40 years of natural time.! This test when used 
on the B Wire Connector indicated that the connector contacts were 
not affected. Although encouraging, we do not consider that this test 
necessarily approximates 40 years relaxation in the connector, since the 
elevated temperature may well have changed the true force pattern 
significantly. For instance, in polyethylene-insulated conductor joints 
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the plastic insulation softens and flows within the connector, changing 
the internal force distribution. 

The philosophy used in applying accelerated tests to the B Wire Con- 
nector, and attempts to correlate these with natural aging, are based 
on correlating equal increments of change in resistance rather than on 
correlation of a known mechanism, such as stress relaxation, to expected 
service life. The aging mechanisms are complex, but observations of the 
resistance behavior itself with time is believed to be more meaningful in 
this case than a fundamental analysis of the contact members them- 
selves. About a year is required to establish with reasonable accuracy 
the correlation between accelerated and natural aging relationships 
using equal increments of resistance change as the criterion. This cor- 
relation has been made possible by the development of the precise meas- 
uring techniques already described. In any of these accelerated tests, 
however, a certain degree of engineering risk is taken when extrapola- 
tion to long periods of time is made. 


5.4 Reliability Appraisal 


As discussed earlier, the reliability and life requirements for a con- 
nector to be used in cable splicing are very demanding indeed. Con- 
siderable effort has been expended to determine if the final design of the 
B Wire Connector meets these requirements. Statistical techniques have 
been used to predict the performance of large numbers of connections 
based on the performance of a relatively few connectors. The validity 
of these predictions depends on the accuracy with which one can define 
the nature of the distribution which describes the population of the 
connections and on how this distribution changes with time. In studying 
the resistance data for the B Wire Connector, it appears that the log 
normal distribution provides a fairly good fit. Analysis of contact re- 
sistance data from simple configurations, such as the classical crossed 
cylinder, supports this conclusion. In this respect, our conclusions are 
not directly comparable to those reached in the evaluation of the solder- 
less wrapped connections, where a Poisson distribution is assumed.? It 
should be noted that in the present work we are describing the absolute 
value of the resistances with a log normal distribution, whereas in the 
work on the solderless wrapped connection, the degree of resistance 
instability with mechanical disturbance was described with a Poisson 
distribution. 

The reliability of the connector after 40 years of aging in the outside 
plant environment is of prime importance. As a starting point we have 
established from resistance measurements on a few hundred connections 
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the initial reliability of the connectors. From the log-probability plot 
shown in Fig. 13, it would appear that less than one connection in 100,000 
would have an input “loop” resistance as great as two milliohms when 
first made. As discussed previously, the loop measurement, in addition 
to measuring contact resistance, also includes the resistance of the con- 
ductors that are within the connector. Due to the complexity of the 
contact configuration, it is difficult to determine precisely what portion 
of the loop resistance is wire resistance. However, from a study of the 
geometry, the maximum and minimum possible wire resistances within 
the connector can be calculated, and are 1.73 and 0 milliohms, respec- 
tively, for 24-gauge conductor joints. In practice the wire resistance 
will always be greater than zero, or conversely, the pure contact resist- 
ance will always be less than that plotted in Fig. 18. To this extent, the 
reliability prediction in terms of pure contact resistance would be better 
than that quoted above. 

Although it may be somewhat premature, we believe predictions of 
the future reliability of these connections can be made with reasonable 
confidence based on the early test results of the final design. Irom re- 
sistance measurements during a 7000-hour period on a small sample of 
connections, aging characteristics for the mean loop resistance and 
standard deviations have been established, as shown in Fig. 14. From 
the extrapolation of these plots, it is predicted that after 40 years the 
mean loop resistance and standard deviation will be 0.780 and 0.110 
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Fig. 13 — Initial resistance distribution of 100 B Wire Connector joints. 


SOLDERLESS WIRE CONNECTOR 151 








24-GAUGE POLYETHYLENE 
INSULATED CONDUCTORS 

















INPUT LOOP RESISTANCE IN MILLIOHMS 























HOURS AGING AT LABORATORY AMBIENT 


Fig. 14 — Performance of B Wire Connector joints under natural aging. 


milliohm, respectively. The distribution of resistance of these connec- 
tions initially and after 7000 hours of aging at laboratory ambient are 
shown in Fig. 15. From the predicted values of the mean resistance and 
standard deviation obtained from Fig. 14, a probability plot showing 
the distribution after 40 years has been constructed, as shown in Fig. 15. 
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Fig. 15 — Predicted changes in the distribution of B Wire Connector joints 
with aging. 
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T'rom the extrapolation of the log-probability plots, it is predicted that 
about 1 connection in 10,000 can be expected to have a loop resistance 
as great as 1.3 milliohms initially, and 1.5 and 1.7 milliohms after ap- 
proximately 1 year and 40 years, respectively. 

The tests referred to in this paper are only a few of the many that have 
been made to optimize the design and study the performance of the B 
Wire Connector. The other tests include studies made with connectors 
having only a few penetrators to downgrade the connector performance 
and make the resistance more sensitive to contact degradation, connec- 
tors made with reduced holding forces as a result of thin outer shells, 
and connectors without tinned liners to accelerate oxidative effects at the 
critical contact areas. This whole body of test results supports the reli- 
ability estimates given here. It is believed that all of the important aging 
mechanisms have been identified and sufficiently taken into account in 
the design. The device has been released for field use and so far the in- 
herent reliability is at least as good as expected. 


VI. CONCLUSION 


The engineering risks involved in appraising a device design for very 
high reliability and long service life can probably never be reduced to 
zero. The techniques described here are believed to have minimized 
these risks in the B Wire Connector, however, and at the same time they 
have permitted substantial reductions in the laboratory effort required 
to develop and prove-in this device. Early field experience is considered 
good confirmation of the basic statistical predictions of initial reliability 
that were made. Final confirmation of the aging predictions must, of 
course, await the passage of time. 

The experimental techniques described are based largely on the use 
of laboratory measurements of improved sensitivity which permit the 
collection of data suitable for more rigorous statistical and experimental 
analysis than otherwise possible. The data analysis itself has two distinct 
goals: (1) the study of the fundamental mechanisms which influence the 
physical behavior of the device, and (2) the formulation of reliability 
and aging estimates based on a minimum of experimental data through 
the use of various established statistical and other analytical techniques. 
In this broad sense, then, the general principles employed in the de- 
velopment of the connector are not new. The degree to which these basic 
principles have been applied in the past to the design of connectors has 
been limited, at least in part, by the experimental difficulties in making 
resistance measurements on actual devices to a sensitivity of much better 
than +0.5 milliohm. The measuring circuit and apparatus discussed 
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here are simple and permit an improvement in resolving power or re- 
finement in sensitivity by at least two orders of magnitude. The work 
on the B Wire Connector clearly shows that important information is 
lost without this refinement. Many of the specific methods presented 
here are directly applicable to a broad class of engineering problems 
having to do with electrical contacts. 
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The Fabry-Perot Electrooptic Modulator 


by E. I. GORDON and J. D. RIGDEN 
(Manuscript received October 10, 1962) 


The Fabry-Perot modulator, consisting of Fabry-Perot etalon plates 
separated by an electrooptic material such as KDP, is analyzed in detail. 
Time-dependent perturbation theory is used to describe the coupling of the 
axial modes by spatial and time varying perturbations in the dielectric 
constant. The perturbations are produced by the applied microwave modulat- 
ing field. It 1s shown that the correct choice of the spatial variation of the 
microwave modulating field is essential to achieve efficient modulation and 
the choice 1s equivalent to matching the phase velocities of the microwaves 
and the light. 

Power requirements, heating, and bandwidth are discussed and a com- 
parison 1s made to the traveling-wave modulator described by Kaminow. 

Calculations indicate that bandwidths of several hundred megacycles, 
centered at any microwave frequency, can be obtained with the expenditure 
of several watts of modulating power. 


I. INTRODUCTION 


Light modulators, operating at microwave frequencies, have been 
receiving considerable attention. Recently, Kaminow,! using the Pockel’s 
effect in KDP, produced usable amounts of modulation at an X-band 
frequency. The modulator was operated on a pulse basis because of the 
large power dissipated per unit volume in the KDP. One solution for the 
problem of heat dissipation as discussed by Kaminow! and others, 
requires careful matching of the phase velocity of the microwaves to the 
light velocity in the electrooptic material, and calls for the construc- 
tion of a long modulator. The resulting large coherence volume decreases 
the power dissipation to a value low enough to allow continuous opera- 
tion. 

The mechanical difficulties, as well as the transmission loss associated 
with long rods, among other problems, has made the realization of such 
low power modulators impractical. On the other hand, the advent of 
sources of monochromatic light energy, such as the optical maser, pro- 
vides alternatives not available with conventional light sources. 
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Fig. 1 — Schematic of Fabry-Perot modulator. 


We refer, in particular, to resonance or cavity techniques analogous 
to those employed at microwave frequencies. Here, the coherence 
volume of the material is made effectively larger by virtue of the multi- 
ple reflections within the resonant cavity. The particular technique 
which we wish to describe is the use of a Fabry-Perot etalon (FPE) in 
which is included an electrooptic material, as shown in Fig. 1. Although 
the bandwidth over which modulation can be accomplished is limited 
by the use of an optical cavity, a bandwidth of several hundred megacy- 
cles centered at any microwave frequency should be readily attainable. 

Before proceeding with an analysis of the modulator it is worthwhile 
to consider some simple points of view which indicate the correctness 
of this approach. The relative transmission J; through a lossless FPE 
can be written’ 


I, = {1 + 472(1 — I’) sin’ g]” (1) 
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in which I” is the power reflectivity of the reflectors and ¢ is the phase 
shift for a single pass. The phase shift, for normal incidence, can be 
written ¢ = 2afLu/c + 6, in which f is the optical frequency, L is the 
etalon spacing, c is the vacuum light velocity, u is the index of refraction 
and @ is the phase shift upon reflection.’ It should be possible to con- 
struct a modulator in which L, @ or yu is varied, thereby varying ¢ and 
modulating the transmitted intensity. We shall be concerned with 
variations in » through the electrooptic effect. In particular we will 
be interested only in phase shifts associated with volume effects rather 
than the surface effects associated with variations in reflectivity. The 
latter are many orders of magnitude less significant. 

In the absence of modulation, the transmission characteristic as a 
function of frequency is given by (1) and is shown in Fig. 2. The trans- 
mission peaks occur at frequencies fa = ac/2Ly for which g = az. 
The integer a is the number of half wavelengths in the cavity. The 
frequency difference between adjacent peaks is c/2LZy. For the purpose 
of discussion assume that the frequency of the incident light corresponds 
to one of the transmission peaks. Assume also that the index of refrac- 
tion has a varying component that takes the form of a cosinusoidal 
standing wave as shown in Fig. 3(a), with a frequency wm/2r = c/2Lp 
and wavelength 2D. The variation is equivalent to oppositely traveling 
waves with phase velocity equal to the light velocity, c/u. With a peak 
variation du COS wmt, the amplitude of each wave is 46u. Under these 
circumstances, the perturbed phase shift travels in synchronism with 
the light producing a cumulative phase shift, 


50 = afhalc ‘du COS Omt = 4am(du/L) COS Wt. 


TRANSMISSION 





fa-2 fa-; fy Fass 2 fa42 
FREQUENCY 


Fig. 2 — Transmission vs frequency for a Fabry-Perot etalon. 
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Since 6u/z is of order 10-° for practical cases, and az is of order 10°, 
it follows that 5g < 10‘ and (1) can be written 


Let + MP1 — 1)? (arbp/p)* cos” wnt)’. (2) 


Choosing the reflectivity T’ = 0.9 indicates that the coefficient of 
COS’ wmt can easily achieve the value unity, leading to large intensity 
modulation. Note that the modulation rate is twice the modulating 
frequency, which is consistent with the fact that the transmitted in- 
tensity is an even function of the induced phase shift. 

From (1) the separation of half-power frequencies or pass band of 
the FPE, shown as Af in Fig. 1, is given by 


Af = fa(1 — T’)/anT. (3) 


For a one centimeter long etalon, with I” = 0.9, and c/2Ly = f,/a = 10", 
Af ~ 300 mc/s. One wonders whether this represents a limiting rate 
at which the light can be modulated, since the pass band Af represents 
the maximum rate at which the stored energy in the cavity can be 
varied. When the perturbation in dielectric constant is uniform across 
the etalon, it will be shown that the maximum possible modulation 
rate is of order Af/2. However, when the phase velocities are matched 
as described above, there is no limiting modulation rate. 

We can obtain an intuitive appreciation of the lack of a limiting 
modulation rate by considering Fig. 3. The perturbed part of the dielec- 
tric constant ¢«, shown as de in Fig. 3(a), is depicted as a standing wave. 
This standing wave corresponds to the microwave modulating electric 
field which is assumed to have a field component normal to the reflectors. 
This component produces the variation 6«. The optical field, excited by 
the incident light, is shown as the standing wave, /, in Fig. 3(b). 
In addition to the displacement current, «d/,/dt, there will be a per- 
turbed displacement current, 05eH,/dt, shown in Fig. 3(c). It can be 
seen that the perturbed displacement current has the appropriate 
spatial distribution to excite the modes E,,;: shown in Fig. 3(d) and 
3(e). 

The microwave modulating frequency fm = wm/2m is chosen to have 
the value c/2Zy, so that the sum and difference frequencies appearing 
in the time dependence of the perturbed part of the displacement 
current correspond to the frequencies, fri:. AS a result, components 
of the perturbed part of the displacement current have the proper 
spatial variation and frequency to excite the sideband modes. The 
transmitted light will then have upper and lower sideband frequencies 
which result from the oscillating sideband modes. 
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Fig. 3 — Mode coupling resulting from perturbations in the dielectric constant. 


From the point of view of coupling of modes, there is no inherent 
limitation to the rate at which modulation can be accomplished, except 
that afforded by the electrooptic effect itself. However, the range over 
which the microwave modulating frequency can be varied must cor- 
respond to the bandwidth of the FPE modes, otherwise the sideband 
modes cannot be excited strongly. 

The conditions under which modulation is possible can be summarized 
by the equations 


Jai = ta + fim 
Bam = Ba = Bm 


in which 6, is the propagation constant of the nth mode, in this case 
fm and Bm are also the microwave frequency and propagation constant. 
Equation (4) is recognizable as the Tien w — £* relation for reactive 


(4) 
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frequency mixing. Equation (4) also implies that the microwave phase 
velocity equals the light phase velocity, assuming negligible dispersion. 

It is not clear yet that the variation in transmitted intensity given 
by (2), which follows from relatively simple considerations, is consistent 
with the coupled mode picture. The next sections will be devoted to an 
analysis of the modulator using time-dependent perturbation theory to 
determine the transmitted field amplitudes and phases. The variation 
in transmitted intensity will be computed and shown to be consistent 
with (2). Other modes of operation will be described which enhance the 
modulation efficiency and which can yield linear intensity modulation 
at the microwave frequency. Power requirements and heating will be 
considered for an X-band modulator for light at 6328 A.t Some of the 
optical requirements of the Fabry-Perot cavity will be discussed with 
respect to available maser beams. Finally a comparison will be made 
to the traveling-wave modulator described by Kaminow. 


II, THEORY 


The analytical description will be similar to that used in time-de- 
pendent perturbation theory. The optical fields within the IPE will 
be expanded in a set of eigenfunctions appropriate to perfectly reflect- 
ing boundaries. It will be seen that the expansion coefficients are time 
dependent and satisfy “harmonic-oscillator-like” linear differential 
equations which are coupled by the perturbed part of the dielectric 
constant. The mixing or coupling terms involve integrals over the cou- 
pled modes and the perturbed part of the dielectric constant. 

The microwave fields will not explicitly enter into the calculation; 
rather they will serve only to produce the perturbed part of the dielec- 
tric constant. Although the dielectric constant in electrooptic materials 
often exhibits tensor properties we will suppress this feature since the 
notation would be burdened without any real gain in generality. 

The basis of the calculation will be the normal mode formulation 
given by Slater.° In the absence of the perturbation in the dielectric 
constant the electric and magnetic fields within the FPE are divergence- 
free. Thus, we will use only solenoidal eigenfunctions which are defined 
in the volume V appropriate to the FPE by the equations 


k.Ea(r) = curl H,(r), k,H.(r) = curl E,(r) (5) 


in which k, is the eigenvalue of the ath mode. A consistent set of hound- 
ary conditions on E, and H, are taken to be 


+ As the modulation efficiency increases with decreasing wavelength, we have 
used the wavelength of the presently shortest wavelength maser.® 
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nXE,=0, n-H.=0 on S, 6) 
nxXH,=0, n-E,=0 on 8S, 


in which n is a unit vector normal to the outer surface of the volume V. 
The surface S, corresponds to the short circuit boundary while S, is 
the open circuit boundary; together they represent the total surface 
of the volume V. Under these conditions, Slater shows that the func- 
tions E, and H, represent an orthogonal normalized set satisfying the 
equations 


VE, tkhevE, =0 and VH, + k.H. = 0. 


The orthogonality relations take the form 
i E,: E,dV — / H,:HedV = Cab 
V V 
The starting point for the analysis is the Maxwell equations 


curl E = —y.dH/ot 


(7) 
curl H = de’E/dt + cE 


in which o is the volume conductivity associated with the transmission 
loss of the electrooptic material, wo is the vacuum permeability and ¢’ 
is the permittivity. As shown by Slater the various terms in (7) have 
expansions of the form 


E(r, t) = > e.(t)E,(r) 


a 


H(r, t) = 2 ha(t)Ha(r) 
curlE = )/ H, i 7 / (n X B)- Huds | (8) 


curl H = >) E, | ta + i (n X H)-E,dS | 
a So 


in which dS is an element of the surface. The mode amplitudes e, and 
ha are defined by 


ea(t) = | BCs, )-EdV 
(9) 
Oe i: H(r, t)-HydV 


in which dV is a volume element. Substituting the relations in (8) into 
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(7) and using the orthogonality and normality of the functions E, and 
H, yields equations for the coefficients of the form 


Ss 
and 


€0€a/0t + cea — kaha — / (n X H)-E.dS 
So 


(11) 
+0 >> (a | de | b)e/at = 0 
b 
in which we have written e’ = ¢ + de, « being the unperturbed part of 
the dielectric constant, and we have used the notation 
(a | 40) |b) = | By Bay(rav. (12) 


Iterating (10) and (11) yields 
cea /dt + ka €a + poodea/dt + ca" >> (a | de/e | b)en/at? 
b 


= Hod (J (n X H)-Buds) / a — be f (n X E)-H.dS (18) 


in which c’ = (oe)? is the unperturbed velocity of light in the dielec- 
tric. A similar redundant equation obtains for h,. The surface integrals 
contain two types of field terms: (i) external optical fields incident on 
the cavity as they appear on the inside boundary wall of the reflector, 
which serve to excite the FPE modes, and (ii) internal fields which 
account for energy lost by radiation, in particular the transmitted beam. 

In the following treatment of the modulator we will assume that the 
reflecting surfaces of the etalon can be characterized by a reflection 
coefficient for field amplitude I and that the surface is a “short,” i.e., 
a surface S,. The fields will be assumed to be parallel to the etalon 
which is of sufficient lateral extent that the axial modes of the FPE 
are essentially TEM waves with no diffraction loss. The dielectric 
modulation is assumed to be uniform across the cross section. Later we 
will consider nonuniform excitation of the dielectric in the finite cross 
section, and it will be seen that the lowest-order axial modes can be 
coupled to the off-axis higher-order modes. 

Under these assumptions, the electric and magnetic field on the 
inside face of the reflectors constituting the surface S, can be shown to 
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be related by 
nX E = c'mH(1 — )/(1+T) (14) 


neglecting the extremely small variation in I produced by the variation 
of the dielectric constant within the etalon. First we consider the loss 
terms in the cavity. For the right-hand side of (13) we may write, 
using (14) and (8) 


—ky {| (n X E)-H.dS 
“Ss 
(15) 
= — ac’ po( 1 — Te + 1 ie > ms | H,:-H.dS 
b Ss 


there being no surface S, over which n X H-E, is nonzero. Substituting 
(11) into (15) to eliminate the term in h, yields 


—Kea [ (n X E)-H.dS = —ke (1 —T)(1 +7)” 
x 2. ky * [de,/dt + o¢s/e + 9 » (b | de/e | d)ea/at] (16) 


[ H;-H.ds. 
Ss 


For the axial modes of the FPE, the appropriate normalized eigen- 
functions are given by 


H, = (2/AL)’ cos kaz 


, (17) 
E, = (2/AL)’ sin kz 


in which A is the cross-sectional area, L the length of the modulator, 
and ka = za/L. The integer a is the number of half wavelengths within 
the etalon. Remembering that there are two faces of the etalon com- 
posing the surface S, , we note that the integral 


: H,-HidS = 4/L (18) 
Ss 


when 0b differs from a by an even number, while the integral has the 
value zero when b differs from a by an odd number. The implication is 
that the surface loss terms couple only modes of the same longitudinal 
symmetry. The integral is zero when H, is one of the off-axis FPE 
modes because of the orthogonality. Defining 
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Qo = tma(1 + 1)’/(1 — Tr’) (19) 
(16) becomes 


seh (n X E)-H.dS 
(20) 
= —k,c De Qs [de,/dt + ce/e + 8 d (b | 5¢/e| d)ea/dt] 


in which the prime implies that only those values of b which differ from 
a by an even number are to be included in the sum. The third term in 
the square bracket of (20) is of order de¢/e compared to first. The second 
term is of order ¢/ew, ~ 10-* compared to the first. Since the first term 
is as small as any other term in (13) because of the factor Q,-! ~ 10-5, 
we can safely neglect the other terms. In the absence of internal losses 
the term Q, is known as the quality factor of the cavity. It defines the 
resolution of the IPE as an interferometer and defines the bandwidth 
by the relation Af ~ f./Q. which differs from the bandwidth defined in 
(3) by the factor 47/(1 + I)? & 1.f 

Considering the incident wave amplitude as transformed to the inside 
surface of the reflector it can be shown in similar fashionf{ that 


—h, | (n X E)-HidS = —4kuc’ "Qa ‘de/dt (21) 
Ss 


in which e is proportional to the magnetic field of the incident wave on 
the surface of the reflector. The factor 4 appears because the incident 
fields appear only on one face. Thus (13) becomes 


O°ea/dt + wala + waQad OCa/Ot 
+ wa >,’ Qs den/dt + 8° D> (a | 5e/e | b) en/at? (22) 
b b 


= —1y,Q, de(wt)/at 


in which w, = k,c’ is the resonant frequency of mode a and Qaa = wae/o 
is the dielectric quality factor or inverse loss tangent. An alternate 
expression can be written Qaa Y waf/c’ in which ¢ is the distance for 
which the energy decays to e. 

Next, we evaluate the quantities (a | (d¢«/e) |b). Since the spatial 
variation of de/e can take the form of sin mzz/L or cos mrz/L in which 
m is some integer, which corresponds to the appropriate microwave 


sl The factor 4r/(1 + I)? differs from unity by less than one per cent for I? > 
0.7 


t The mirror acts as a transformer with turns ratio (1 + 1T)3/(1 — T)3. 
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field variation, we evaluate two types of coupling terms: 








/_| sin wmz/L},\ 
\" | cos amez/L o 
' — (23) 
= (2/L) i (sin raz/L a da sin xbe/1) dz 
in which we have used (17). Performing the integration yields 
(a | (sin rmz/L) | b) 
=5|5 — cos (a —b + mq] , El — cos (b-—a+m)za] 
2 (a-—b+m)r (b—a+m)r (24) 
Uo = IG eel a Ee ee 
(a+b—m)r (a+b+m)r j 


Terms like (a + b + m)7 will be of order 10° and the last two terms 
can be neglected. When a — b + m = O or an even integer the remain- 
ing terms vanish. When a — b + m is an odd integer the terms are 
finite, yielding 


2/4 
(a — b)? — m?’ (25) 


a —b-+m = odd integer 


(a | sin“) = — 


the largest values corresponding to a = b, m = 1 with a value 27’. 
Likewise, 





mre _1) sn(a+b—™m)x , sin (a — b— m)r 
(a | cos 2 |b) = 3] (atb—m)r (a—b—m)z 
4. sin (a—b+m)nr _ sin(a+ b+ me (26) 
(a—b+m)xr (a+b-+m)xr 
=4 for ja—b|=™m. 


The sine and cosine cases are mutually exclusive. Many modes are 
coupled by the sine variation since for a given value of a and m there 
are multiple values of b satisfying (25). There are only two values of b 
satisfying (26) for the cosine variation. Note that for m = 0, corre- 
sponding to uniform perturbation of the dielectric, the only coupling is 
for a = b. That is, the modulation can occur only within the mode. 
This justifies our earlier statement that the modulation rate is re- 
stricted to the bandwidth of the FPE mode for uniform excitation. 
The cosine distribution yields the largest coupling terms except for 
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the case when a = b and m = 1. The case m = 1 implies that the sine 
variation has a space average part and it is this part which produces 
the coupling for a = b. This is borne out by the fact that for a = |, 
m must be odd and the strength of the coupling decreases with increas- 
ing m. Notice that even though the modes of interest are standing 
waves these results correspond precisely to that predicted by the 8 
condition of (4). The factor of $ for the cosine distribution arises be- 
cause only half the amplitude of the standing wave is effective for a 
given traveling-wave component. The factor 27’ for the sine variation 
can be shown to arise from the fact that although the 6 condition is 
not precisely satisfied, the interaction with each traveling-wave com- 
ponent is down by 7. 

In the next sections, we apply these results to calculate the amplitude 
of the modulation sidebands, their bandwidth, and the intensity modu- 
lation of the light. 


III, AMPLITUDE OF THE MODULATION SIDEBANDS 


Assuming that we have used the appropriate cos mzz/L distribution 
for that component of the microwave field which varies e’, since we 
want large coupling for a ~ b, we can write (22) as 





v : 1 0g , 10” de(wmt) 
EY, €a + Wa Ca + waQuaa eT; 9 0 = eu" [Casas + Caml (27) 
, —1 dep at 2 —1 de(wt) 
1 ta 2d Qo a aaa Say 


which holds for all values of a.f We note that modes é. and éaym are 
coupled together. Since the Q of the modulator is so extremely high, 
it is approximately correct to say that any driving terms at frequencies 
outside the pass band of a given mode have negligible effect on that 
mode. Later we will see that for very strong modulation (a case of 
little practical interest) the high-Q approximation leaves a little to be 
desired. Hence, to good approximation, we may write for any of the 
modes a, 


dy’ Qu tae,/dt ~ Qa 10e,/dt. 
b 


All the neglected terms are at frequencies well outside the passband of 
mode a. With these definitions, the quantity 1/Qar = 1/Q. + 1/Qua 
+ It is also possible to have | a — b| = 1, m = 2 yielding a coupling term 4/37 


which is only slightly less than 3. We will not consider this case which is consid- 
erably more complicated. 
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is properly called the loaded Q of the resonator. The bandpass in the 
absence of mode coupling is defined by Aw = o.u/Qar. The quantity 
2Q. is also properly called the external Q. 

We will assume that the carrier frequency corresponds to the reso- 
nant frequency of the ath mode, w = k,c’ = wa. There is no loss of 
generality in this assumption since the carrier frequency would nor- 
mally be fixed. The microwave modulating frequency, on the other 
hand, must vary over some frequency range. Ideally the microwave 
frequency has the value mc’/2L, since under this circumstance the 
mixing terms, represented by the fourth term of (27), vary at rates 
corresponding to the resonant frequencies of the coupled modes. How- 
ever the microwave frequency must be able to vary over some fre- 
quency range. We define the microwave angular frequency as 


om = (mre'/L) + wad/Qar (28) 


and the parameter 6 measures the excursion from center band frequency. 
In view of the fact that wa/Qar ~~ Awa, the half power frequencies are 
given approximately by 6 = +4. The + nth sideband frequencies are 
given by 


Wa NWm = Wa + ninwe'/L + nbwa/Qar = Watnm ~ NbWa/Qar. (29) 
The pertinent mode amplitudes are defined by 
Casnm = Latnm CXP 1(Wa + Nwm)t + complex conjugate 
de/e = [e/e] EXP tumt + complex conjugate (30) 
e = FE exp tw,t + complex conjugate 
which when substituted into (27) yield 
—[x*Bapm + xHa—m)] + tha = — 31H Qar/Qa (31) 
QndHo-nm — [x*EBo—naym + xBo—a+yml + tHe-am = 0 (82) 
—2nbHa+mn — [x*Bapasym + xLatn— ym) + Batam = 0. (33) 


The coefficients are complex quantities and the asterisk denotes com- 
plex conjugate. The modulus of the coefficients equals half the peak 
value and the argument represents the phase of the related real variable. 
The quantity 


xX = 7Qar[de/€] (34) 


is a “modulation index” whose modulus equals the product of Qa: and 
one quarter of the peak variation of d«/e. The argument of x equals the 
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phase of the microwave field. Note that we have set all the terms 
Qainm equal to Q.. Reference to (19) indicates that the variation in 
Q. from one mode to the next is of order a’ & 10° and thus the sub- 
script is superfluous. Likewise, all other terms of order zc’/2Lw, = a 
have been neglected. 

It is helpful to observe that (33) is redundant. If we define 


ype = (—1)"Ba-nn*™ (35) 


and substitute into (33), we obtain (32). 
Defining the normalized field amplitude 


Gn = (1/x*)"Eanm/(ZLQa1/Qa) (36) 

and substituting into (31) and (32) yields 
2\x PRegi + go = —1 (37) 
(1 — 42n8)gn — Gra + |x Pons = 0 (38) 


which is a three term recursion relation with nonconstant coefficients. 
Placing the question of the microwave bandwidth aside momentarily 
we may assume that the microwave frequency is at its appropriate 
center band value corresponding to 6 = 0 in (38). Under this condition 
solutions for the recursion formula, (38), of the form g, = b”, yield 
as the defining equation for b, | x ’b” + b — 1 = 0. The general solu- 
tion may be written g, = Ab,” + Bb_” with 


olan V law) tel 


We note three things about g,: (i) it is real, since go is real and all the 
coefficients in the recursion formula are real; (ii) it is a function of 
|x |’ and (iii), since forn > 0, (Ha-nm/H) — 0 as the modulation in- 
dex | x | > 0, so must g, | x |" — 0 as follows from (36). 

The solution b_ has the dependence b_” ~ | x |"°" and therefore does 
not have the correct behavior as | x | — 0. Thus, we choose the constant 
B = 0. The constant A is evaluated by substituting g, = Ab” into 
(37), yielding A = —(1 + 4|x |’), so that we may write 


] 
n= —(4)” 
2 2" Tx 





IV1+4[xP-U/Vi+4[xP (40) 
and from (36) 


Eq_nm/$E = —(i12x)~"(Qar/Qa) [M1 + 4x? — 1)” 


ees LD 
+V71+4/x/. 
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The same expression with x — x* applies for Hasnm/H. We note that 
the amplitude of the carrier wave, n = QO, is given by 


E,/3E = —(Qar/Qa)[1 + 4|x[T? (42) 


and is reduced in amplitude from its value of —Q.1/Q, in the no modu- 
lation limit | x | = 0. The value Z, = —4H in the no modulation case, 
with Quar = Qa corresponding to Qaa = ©, is eminently reasonable. 
Remembering that the surface is a “short” the implication is that the 
magnetic field of mode a of the resonator exactly cancels out the half 
of the total carrier magnetic field / on the inner reflecting surface 
which is associated with the reflected wave of amplitude #/2. The inci- 
dent wave also has amplitude //2, yielding a total field H. Likewise 
the magnetic field on the opposite reflecting surface is +H/2. Thus, 
the incident carrier field 1/2, which measures the incident energy, 
correctly accounts for the energy leaving through the opposite face, 
and we have 100 per cent transmission through the etalon and no 
reflection. In microwave parlance, the resonator is matched to the 
incident wave. This always occurs when the external Q equals twice 
the loaded Q. In the presence of internal loss, Qaa < ©, the resonator 
cannot be matched if the etalon plates have equal reflectivity. In the 
presence of the modulation, the carrier mode is loaded by the side- 
bands even without internal loss and the FPE is no longer matched 
to the carrier; consequently carrier energy is reflected as well as trans- 
mitted. This shows up as a reduction in transmitted power as well as 
imperfect cancellation of the reflected wave since | Z,| < | #/2|. 

The relative amplitude of sideband to carrier leaving the opposite 
face of the etalon is given by 


(Ea-nm/Ea) = (12x) "[V/1 + 4 [x — 1)". (43) 


In the limit of small | x | 
(Ha-nm/Ha)|x120 = (x*/t)" (44) 


and the amplitude is down by | x |". The power in the first sideband 
is proportional to | x |’ which is proportional to the microwave modu- 
lating power. For large | x | 


(Ea-nm/Ea)\x1>0 = (|x |/ix)” = exp — in(a/2 + Arg x) (45) 


so that the modulation approaches 100 per cent. As we may note from 
(41), the absolute amplitude of all of the modes decreases with in- 
creasing |x| as (1 + 4|x]?). Thus we might expect maximum 
sideband energy for some finite value of x. Maximizing E._»,/E shows 
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that for the nth sideband the value of | x | for which the energy is a 
maximum is given by the positive real solution of 


16|x — 4n* |x? — nn’ =0. 
For the first sideband the maximum occurs at 


|x| = 4[2(1 + /5)} = 0.636. 


The transmitted light beam contains 23 per cent of its energy in the 
first sideband, and in the absence of internal loss the transmitted carrier 
energy is reduced to 38 per cent of the original carrier energy. As will 
be seen later, this modulation can be accomplished with only a few 
watts of microwave power over a bandwidth of several hundred mega- 
cycles at X-band. 

The microwave bandwidth over which modulation can be accom- 
plished is determined from a solution of (38) with 6 ~ 0. The value of 
6 for which the energy in the first sideband falls to half its value at 
center band is determined by a solution of | gi(5)|’/| gi(0)|? = 3. The 
details of this calculation are too lengthy to be included here. The calcu- 
lation was performed by first noting that when the first sideband mode 
is driven near its band edge the second sideband mode is driven well off 
resonance. As a result, the values of g,(6) decrease rapidly as 6 — 4. 
It follows from (38) that setting g,(6) = 0 is equivalent to neglecting 
terms of order (2n)”” « 1. We assumed that g; = 0 and solved the 
three relations that follow from (37) and (38). The result can best be 
expressed in the form 


6= 43(14+4|x [P+ ---) (46) 


which indicates that 6 is increased over its zero modulation value of +4. 
This result is understandable if we consider that the modulation process 
tends to load the carrier mode and decrease its amplitude. This process 
is more effective when the sidebands are excited strongly. The sideband 
amplitude is proportional to the carrier amplitude; thus the center of 
the sideband modes decreases more rapidly than the sides for a given 
modulation index yielding, in effect, a broader modulation bandwidth. 
From another point of view we may say that the coupling of a given 
mode to adjacent modes increases the rate at which energy can be 
dissipated in the given mode. This decreases the effective loaded Q and 
increases the bandwidth. 

For radiation at 6328 A, f, ~ 5 X 10” eps, fm = 10” eps, it follows that 
a = 5 X 10‘ and for I’ = 0.9 we have Q, ~ 1.5 X 10°. Thus the modu- 
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lation bandwidth is at least 3 X 10° eps. For large values of | x | the 
actual bandwidth may be considerably larger but the modulation will 
not be linear. In any case the high-Q approximation, which is funda- 
mental in this analysis, is probably no longer valid and all conclusions 
for large | x | must be considered as only qualitative. 


IV. EFFECT OF NONUNIFORM DIELECTRIC EXCITATION 


So far we have assumed that the dielectric constant variation induced 
by the microwave field depends only on z and t. However, it is necessary 
to match the microwave phase velocity to the light phase velocity. 
Since the dielectric constant at the microwave frequency is, in general, 
considerably larger than it is at the light frequency, the microwave field 
cannot propagate as a plane wave with no transverse variation. The 
transverse variation could probably be avoided by making the diameter 
of the Fabry-Perot cavity very small so that the microwave fields are 
only partially slowed by the dielectric, but this is not usually possible. 
As a consequence the microwave field must have a variation of ampli- 
tude across the Fabry-Perot cavity. The resultant transverse variation 
of the perturbed part of the dielectric constant will cause coupling of 
the FPE modes. 

In general we must describe the ath FPE mode in terms of a series 
of transversely varying eigenfunctions. Thus 


E. = >. AimRi(Kim-t) sin raz/L (47) 
lm 


in which A;,, is a constant determined by the geometry. The frequency 
of the /mth mode is 


fam = (¢'/2m)[(ma/L)* + kin’|t = foll + Iin’/(na/L)*}) (48) 


so that each of the modes characterized by the number of half wave- 
lengths equal to a has a different resonant frequency. However, kim 
is of order mz/D in which D is some characteristic transverse dimension 
of the FPE. The ratio kin/(aa/L) is of order (m/a) (L/D). Since 
L/D will be of order ten or less and m/a will be order 10° for the lower- 
order transverse variations, we see that the frequencies f..,, differ from 
the fundamental frequency by about five parts in 10°. For 6328 A 
radiation this amounts to about three megacycles, well within the 
passband of the FPE. 

The microwave modulating signal designed to mix adjacent FPE 
modes will not produce mixing of different transverse modes with the 
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same a number but it can produce transverse mixing between adjacent 
modes. The degree of mixing will depend on integrals of the sort 


Ff RiCbeamet) « (See) /e)Ra (erm tS (49) 


evaluated over the cross section of the FPE. The functions Ri(k)m-r) 
are assumed to be suitably normalized and orthogonal. It is apparent 
that when 6e¢/e has only a slow variation over the cross section then 
both U’ and m’ cannot differ much from / and m. Thus the distribution 
of fields over the output face of the PE will not differ much from that 
associated with the incident carrier wave and the angular spread of the 
output carrier and the modulation sidebands will not be seriously 
deteriorated. In addition the amplitude of the modulation index will 
not be drastically reduced. Thus it appears that even though the micro- 
wave field may vary from a maximum value at the center of the FPE 
to a value approaching zero at its transverse boundaries, no serious 
effects will result. Very rapid transverse variations of the microwave 
field, on the other hand, could cause difficulties. 


V. AMPLITUDE MODULATION 


In the traveling-wave modulator proposed by Kaminow' and the 
IPE modulator described here, the fundamental effect is a phase 
modulation. By suitable use of polarizers and birefringent plates, the 
phase modulation of the traveling-wave modulator can convert linearly 
polarized light into elliptic polarization and hence into amplitude or 
intensity modulation. It is therefore of interest to look into the intensity 
modulation properties of the FPE modulator. We require an expression 
for the total electric field associated with the transmitted light L;. 
Using (28), this is given by 


E, = Eq exp twat 1 + >> ((Eo_nm/Ea] exp — inont 
n=1 

(50) 

+ [Hasam/Eal exp tNomt)] 


in which FE, is assumed to be polarized parallel to one of the privileged 
axes which we shall denote as the z-axis of Fig. 1. Combining (385) and 
(41), (50) becomes for the z-component of the FPE field, 


Ey, = Ear XP twat [1 + 2d, (y” + [-y*]")] 


Ear €Xp twat [11 — y) + 1 + y*)7 - 1] (51) 
Ear €XP twat F(y) 
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in which 
Y = [Lam/EHal exp — tomt (52) 
= (12x) (V1 +4] x? — Uexp — tunt, 
Fy) =[Q—-—y)"+d4+y*)7- 1 
and we have used the identity >, y”" = (1 — y)” — 1. Since the 
privileged axes are driven push-pull, the expression for the transmitted 
electric field when the incident light is polarized parallel to the other 


privileged axis is identical, except that x must be replaced by —x, 
and thus y must be replaced by —y. Therefore, we may write 


Ey = Eay exp twat F(y)*. (53) 


Since the incident field has the value E/2 when the light is polarized 
along one of the privileged axes, the relative intensity of the transmitted 
light is given by J; = | 2E,./E |’, and we can write, using Eqs. (42), 
(51) and (52) 


T, = | F(y)P(Qar/Qa)’/fL + 4| x P)- (54) 
It can be shown after some manipulation that 
[1 — 72 |x| cos (wmt + o)] 
[L + 4| x |? cos (ont + ¢)] 
in which ¢ is the argument of x which corresponds to the phase of the 
microwave field. Substituting (55) in (54) yields 
Tr = (Qat/Qz)'/[L + 4| x |’ cos’ (mt + ¢)] (56) 


for the intensity modulation. For comparison with (2), which is the 
intensity modulation derived from the Fabry-Perot equation, we assume 
that the internal loss is negligible so that Q.r = Q. . Using (34) and (19), 
and writing d¢/e = 26u/u, it follows that 

ma (1+T)’ dy 


Fly) = V1 +4|xP (55) 


We note that the equivalent term in (2) is I (madu/u)/(1 — I”), which 
differs by a factor 4T/(1 + IT)’ & 1 from 2 | x |. The negligible differ- 
ence from unity can be attributed to the high-Q approximation. Thus 
the relatively naive considerations which lead to (2) have been shown 
to be precisely correct. 

Other, more efficient, types of intensity modulation can be achieved 
through the use of polarizers. Suppose that the incident light is polar- 
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ized at 45° to the privileged axis, and that the modulator is followed by 
an analyzer which is set parallel to the polarizer. Under these circum- 
stances Haz = Ka, and 


Ey = (Biz + Ety)/W2 = 42, exp twat ([F(y) + F(y)*I. 
Thus we may write 
Ty = (RF (y)(Qar/Qa)*/{l + 41x PI 
= (Qer/Qa)/tL + |x [P cos” (ont + 9) 


which we note differs from (55) in that the denominator is squared. 
As a result the depth of modulation is roughly twice that attained with- 
out the use of polarizer and analyzer. The extra efficiency arises from 
the fact that the cavity also enhances the induced birefringence and 
this effect is brought out by the use of the polarizers. 

When the analyzer is set at 90° to the polarizer 


Ey = (He — Ey)/V2 = $1, exp twat [F(y) — F(y)*) 


(58) 


and 


Ti = Unk’ (y)P(Qar/Qa)/TL + | x PI 
= (Qar/Qa)"4 | x ? cos” (wmt + ¢)/ +4] x ? cos” (mt + 9)? 


and we note that the light is 100 per cent intensity modulated. The 
peak transmitted intensity occurs for | x | cos (wnt + ¢) = 0.5. 

So far, the intensity modulation has been at twice the microwave 
frequency. It is possible to modulate at the microwave frequency by 
introducing a quarter-wave plate between the polarizer and analyzer, 
with the polarizers either parallel or crossed. Under these circumstances, 


Ey = (Bu + iBy)/V/2 = 28a exp twat [P(y) + iF (y)"I, 


(59) 


and 
I, = RF (y) & InF'(y)P(Qar/Q2)?/lL + 4| x FI 
1(Qar/Qa)” [1 F 2| x | cos (amt + ¢)}/ (60) 


[1 +4] x]? cos’ (wmt + ¢)} 


I 


I 


which has a component at the microwave frequency. The introduction 
of a quarter-wave plate provides a bias, or steady retardation, which 
allows linear modulation. 
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VI. MODULATION POWER 


The microwave power required for modulation is a function of the 
percentage modulation. For the sake of discussion we will assume the 
modulation index is 0.636, corresponding to 23 per cent modulation or 
maximum energy in the first sideband as described earlier. In terms of 
the peak variation of the index of refraction u, the modulation index is 
given by x = 3Qa16u/u, as follows from (34). For light polarized per- 
pendicular to the z-axis of KDP’ and parallel to one of the privileged 
axes 


ou = +1yregH (61) 


in which F, is the peak microwave electric field parallel to the z-axis 
and re; is the electrooptic constant. The peak electric field can be re- 
lated to the microwave power absorbed in the KDP by 


P & 1t0cE,;V (62) 


in which P is the absorbed power, V is the volume of KDP and a is the 
microwave conductivity given by 


a = we/Qa (63) 


in which w is the angular microwave frequency, ¢ is the microwave di- 
electric constant and Qa. is the dielectric Q or inverse loss tangent. It 
has been assumed that the energy loss associated with microwave elec- 
tric fields transverse to the z-axis is negligible, a condition which is easily 
satisfied. The factor } arises from a time and space average of the micro- 
wave standing wave pattern. Combining (61) to (63) yields 


P & (we/Qa)(2x/Qaru'res) V. (64) 


At a frequency of 10” cps (X-band) the dielectric constant « © 20c0, 
in which ¢ is the dielectric constant of free space, and Qa & 150.+ For 
fa = 5 X 10” ops, res & 10°" meters/volt, I’ = 0.9, Qa, & 10° and p & 
1.5, (64) yields a power requirement of 240 watts/cc of KDP for | x | = 
0.636. For modulation in the first sideband we must have c’/2L = 10" 
or L = 1 cm. For modulation powers of order ten watts the cross-sec- 
tional area of the KDP must be under 0.04 cm’, corresponding to a di- 
ameter of about 0.2 cm. 

The temperature rise AT from surface to center of a uniformly heated 
rod which is long compared to its diameter can be written 


| Measured at room temperature. 
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AT & LPR?/KV (65) 


in which # is the radius, K is the thermal heat conductivity and P/V 
is defined by (64). For KDP, K = 10” watts/em°K yielding a tem- 
perature rise of about 60°K for a diameter of 0.2 em. Clearly the diame- 
ter of the rod should be made as small as possible. 

Another alternative is that the length of the FPE can be doubled. 
By modulating in the second sideband, the same X-band modulation 
can be accomplished with one-half the power and one-quarter the power 
dissipation per unit volume. The modulation bandwidth would be 
halved. 

The choice of | x | = 0.636 was motivated by the fact that this yields 
maximum energy in the first sideband. I’or superheterodyne systems 
this is a desirable objective. or envelope detection, amplitude modula- 
tion at the modulating frequency is desirable. In this case the combina- 
tion of polarizer analyzer and quarter-wave plate is required yielding 
the modulation described by (60). For this case | x | & 0.1 would be 
quite adequate, and the heat dissipation, temperature rise and required 
power would be reduced by a factor of four. 


VII. OPTICAL CONSIDERATIONS 


It has been supposed that the light beam to be modulated is perfectly 
collimated, thereby allowing the energy to be completely coupled into 
the fundamental IPE mode. Since no light beam is perfectly collimated, 
we would like to examine this assumption and to determine what limi- 
tations are placed on the optical cavity. 

When a diverging beam of monochromatic light is incident on the 
IPE, the transmitted beam will form the typical Fabry-Perot fringes’ 
at infinity. The central fringe corresponds to one of the axial F PE modes. 
Thus we can be assured that the incident energy will excite only the 
axial mode if the angular spread of the beam is less than the angular 
spread of the central fringe. The phase difference between successive 
averaging rays moving at an angle @ within the FPE is given by’ 


20 = 4nfLuc cos 6. 


Substituting into (1) indicates that the transmitted intensity will fall 
to half its @ = 0 value for a bright central fringe at an angle 6, given by 
0, & [e(1 — 1’) / QnfaLply’ ~~ Q.*. Thus the angular spread of the inci- 
dent beam must be less than 2uQ,°. Taking the internal losses of the 
IPE into account, we would find the somewhat larger angle DuQar 


For the FPE cavity discussed earlier with Qa, & 10°, QuQar? 3X 10° 
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corresponding to 10 minutes of arc. For a beam whose angular diver- 
gence is diffraction-limited, the required aperture of the I-PE is 


D & 1.2cQar*/2ufa 


or about 0.25 mm for the example given above. When ‘‘walkoff”’ of light 
rays incident at the angle u@, is considered, the required aperture is in- 
creased by approximately 26,L/(1 — TI’) = 0.15 mm. Since KDP rods 
with a diameter much less than about 1 mm do not seem practical for 
the present, no constraints will be imposed by any diffraction limited 
beam such as that from an optical maser. 

One important constraint imposed on the optical maser is that it 
oscillate in only one longitudinal mode or, if not that, that the modes 
be separated by at least the frequency Af but less than c/2Lu. 

In a manner similar to that by which the half angle 6, was found it 
can be seen that the quantity uL must be controlled such that 


A(uL)/pL < (2Q2)" &3 X 10%. 


Changes in temperature of the cavity will tend to change the product uL. 
Our best guess, in view of the lack of data, is that Au per degree centi- 
grade is of order 2 X 10 ° for one centimeter length. Thus, the tempera- 
ture must be controlled to better than 0.1°C. Likewise, the etalon should 
be flat and parallel to within 3 < 107’ cm. 


VIII. COMPARISON TO TRAVELING-WAVE MODULATORS 


In traveling-wave reactive modulators the light makes one traversal 
of the electrooptic material. The transmitted electric field amplitude, 
for light polarized along one of the privileged axes, is given by 


E, = Ey exp (ty sin wnt) (66) 


in which y = (wL/c’)(6u/u) and Eo is the amplitude of the incident light 
wave at frequency w. The comparison to the FPIX modulator is made by 
relating the amplitudes of the various sidebands. Thus we may write 
for EB t 


$00 


E, = Ey DY Jaly) exp inont (67) 


— 


in which J, is the Bessel function of order n.f The amplitude of the 
nth sideband compared to the carrier is 


-f-00 
+ We have used the identity exp iz sin @ = >> Jn(z) expiné. 
Nm 


178 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


E,/Ey = Jaly)/Jolv) © "7/2"! for yK1 (68) 


which should be compared to (44). We note that for first sideband 
modulation y/2 should be compared to | x |. Thus we compare 


2(wL/c') (du/n) 


and 3{(1 + 1)°/(1 — I°)|(wL/c’)(6u/n). This indicates that the FPE 
modulation is more efficient by the factor 3(1 + T)?/(1 — T°) & 
r(1 — I’)~. For r = 0.9 and equal modulation lengths this is a factor 
of ten in 6u/u, requiring one hundred times the power for equal modula- 
tion. The advantage of multiple reflection can be compensated some- 
what by increasing the interaction length of the traveling-wave modula- 
tor. For example, increasing L by T'(1 — I”) equates | x | and y/2 for 
equal 6u/u but requires (1 — I) times as much power, since for equal 
5u/u the power dissipation per unit length must be the same. Increasing 
L by (1 — I’) & 100 and decreasing 6u/u by (1 — I”) again 
equates y/2 and |x|. In this case the power requirement for cither 
modulator is the same. The potentially greater bandwidth of the travel- 
ing-wave device in this case is probably obviated by the difficulty of 
matching microwave and optical phase velocities closely because of the 
long interaction length. The power dissipation per unit volume is of 
course, lower by I™*(1 — I”)* in the traveling-wave device and this 
is an extremely important advantage. However, the difficulties associ- 
ated with fabricating a long rod of electrooptic material appear to be 
formidable at the moment. In addition, the optical loss of the KDP, 
as has been recently pointed out by Kaminow,’ is not completely negli- 
gible. A rod sufficient to achieve the last mentioned advantages would 
be about one hundred centimeters long at X-band modulation fre- 
quencies. Since the energy decay distance is about 35 centimeters, { the 
light intensity would fall to ~e° of its initial value, i.e., about 13 db 
down. 

In order to determine the transmission loss through the FPE modu- 
lator we may consider the ratio 


(Qar/Qa)” = [1 + Qu/Qaal”. 


Since 
Qa/Qaa & L/t(A — T’) & 2 


the transmitted intensity is decreased by a factor slightly less than three, 
which is comparable to that of the intermediate length modulator. 


jt A rough measurement at 6328 A. 
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Thus a traveling-wave modulator with a length about T(1 — I”’)~ 
times that of the FPE modulator is comparable with respect to heat 
dissipation and transmission loss but requires T'(1 — I’)~* times the 
modulation power. An increase in length by another factor T'(1 — I”)™ 
is somewhat drastic, since the transmission loss increases too rapidly. 


IX. CONCLUSION 


The modulation properties of an electrooptic material placed in a 
Fabry-Perot etalon have been analyzed in detail. The main conclusion 
of the analysis is that the FPE modulator offers considerable advantage 
over the traveling-wave modulator with respect to modulation power 
when the traveling-wave modulator is not extremely long. 

The bandwidth of the modulator is limited by the passband of the 
FPE, which for practical purposes is several hundred megacycles. In 
general this bandwidth is compatible with the bandwidth of the micro- 
wave modulating cavity (without wasting modulation power) since 
the dielectric Q of the material is 150, yielding cavity Q’s of order 75 
and a bandwidth of approximately 130 mc/s. When the electrooptic 
material does not fill the entire cavity or when additional microwave 
bandwidth is required, additional loss must be coupled into the cavity. 
Likewise multimesh broadbanding techniques can be useful. 

The requirement of a stable monochromatic source of light for the 
modulator represents no real disadvantage, since we suppose that the 
modulator would be used in a communication system which needs such 
a source anyway. 

A major requirement for the successful operation of the FPE modu- 
lator is the choice of the correct microwave field configuration which, in 
effect, matches the microwave phase velocity to that of the light in the 
electrooptic medium. 
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Properties of Fast-Decay Cathode-Ray 
Tube Phosphors 


By ARNOLD PFAHNL 
(Manuscript received October 30, 1962) 


Cathodoluminescent efficiency, decay time, aging characteristics, spectral 
distribution, and temperature dependence of luminescent emission have 
been measured for the more intense fast-decay phosphors. Materials studied 
were ZnS(Ag,Ni); ZnO(Zn) = P24; BaPbSO, ; (CaO)2-MgO- (SiOz). - 
(Ce) = P16; and CaO-Al,03-(Si0).-(Ce) = P16. 

The intensities of the P16 and P24 samples were observed, to a good 
approximation, to be proportional to the decay times. Aging due to electron 
bombardment was found to be fastest for BaPbSO, and the P16 CaMg- 
silicate, slower for the P16 CaAl-silicate, and slowest for P24 and ZnS- 
(Ag,N2). The decay was nonexponential for all the phosphors studied and 
the decay time became shorter with increasing beam current density. The 
intensity generally decreased with temperature, reaching 10 per cent of the 
room temperature value at about 200°C, 300°C’, 350°C, 600°C, for the 
ZnS(Ag,N1i), BaPbSO,, P24, and CaAl-silicate, respectively. The CaM g- 
silicate retains 90 per cent of its 20°C intensity at a temperature of 600°C. 

Investigations were also conducted on optimum screen weight, effects of 
particle size distribution, and on screen preparation techniques, all of which 
will be discussed. 

It ts concluded that phosphor characteristics differ significantly between 
batches and manufacturers. A detailed evaluation of the important parame- 
ters should be conducted on individual phosphor batches for applications 
where screen characteristics are critical. 


I. INTRODUCTION 


When selecting a phosphor for cathode-ray tube applications, it is 
important to consider the efficiency, decay time, aging characteristics, 
spectral distribution, and temperature dependence of the luminescent 
emission. The present paper describes the results of measurements of 
these characteristics on numerous samples of fast-decay (r S 1 usec) 
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phosphors of the ZnS(Ag,Ni), ZnO(Zn), BaPbSO,, (CaO),-MgO- 
(SiOz). , and CaO-Al,O3- (SiOs)2 types. The first section deals with ex- 
perimental procedures used in the determination of the phosphor char- 
acteristics. The succeeding sections discuss the specific characteristics 
of the phosphors studied. 

The work described in this paper was undertaken in connection with 
the selection of a phosphor for the screen of the cathode-ray tube of the 
flying spot store, a semipermanent memory developed for use in an ex- 
perimental electronic telephone switching system. 


II. EXPERIMENTAL PROCEDURES 


Many luminescent properties of phosphors depend strongly on the 
mode of excitation because of their dependence on parameters such as 
the absorption coefficient of the exciting radiation, excitation density, 
temperature effects, etc. Phosphor evaluation experiments should 
utilize excitation conditions as similar as possible to those to be used in 
the application. If exact duplication is not attained, it is necessary to 
determine the limits of the validity of the tests. 


2.1 Screen Preparation 


All commercial samples measured were used as received, without 
any grinding or other mechanical processing prior to the settling opera- 
tion. During screen settling, the required phosphor quantity is suspended 
by ultrasonic agitation in a solution of potassium silicate (density 1.07 
g/cm?), This suspension is poured into an aqueous solution of barium 
nitrate (22.5 g/l H.O) which covers the surface to be coated. The set- 
tling time depends on the density of the phosphor and varies between 
one and several hours. The composition of the settling solution is critical 
if optimum screen adherence is to be achieved.! Care must be taken to 
use purified chemicals and water; otherwise screen adherence may be 
poor, or the phosphor may be poisoned by the diffusion of impurities 
during tube baking. 

Screen settling with potassium silicate, K:SiO;, is based on the 
adsorption of silica acid ions (SiO.=) on the surface of the particles and 
the glass. These strong negative charges create a repulsive force which 
prevents coagulation and adherence of the particles. The addition of a 
salt, for example Ba(NOs;)2, decreases the negative potential and en- 
hances the formation of oxygen bridges between the silica acid ions. 
This leads to gelatination and strong adhesion between the particles 
as well as between the particles and the glass.? 34 
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Tor aging measurements, the screens were backed with a layer of 
aluminum 500-600 A thick. This layer acts as an ion trap and also 
reflects that part of the fluorescent light which would otherwise be lost 
in the interior of the tube.® It was found that the thickness of the 
aluminum layer must be greater than 300 A for maximum reflection, in 
agreement with measurements published in the literature.*: 6 


2.2 Measurement of Relative Intensity and Decay Time 


The relative intensity and the absolute value of the decay time 
(1/e point) were measured with a continuously pumped demountable 
tube. This does not approach the environment of a sealed-off tube very 
closely, as the pressure in the system was relatively high (8 to 12 x 
10-* mm Hg). The screens were not aluminized and no processing (heat 
treatment) of the tube was done. Some comparative experiments indi- 
cated, however, that the results were proportional to the values ob- 
tained in sealed-off tubes. 

In agreement with the observation of other authors,’® it was found 
that the use of oil pumps reduces the secondary emission ratio of the 
screens, leading to a rapid decrease of the phosphor surface potential 
and consequently the emission intensity. This effect persisted when a 
liquid nitrogen cold trap was added to the system, but was eliminated 
when the oil diffusion pump was replaced by a mercury type. 

The phosphors were settled on microscope slides. Up to 6 samples 
could be introduced simultaneously into the tube and tested under 
identical conditions. The grid of the electron gun (5W/5Z type) was 
pulsed 1000 times per second with pulses of 5 usec length; the beam 
voltage was 10 kv and the current density in the focussed beam about 
5 ma/cm?. The output light pulse was observed with a photomultiplier 
having an S11 spectral response and displayed on the screen of an oscil- 
loscope. It was possible to expand the display of the decay and to read 
the decay time with high precision directly on the vernier of a variable 
delay (Fig. 1). This method has also been used to determine the form of 
the decay. The shortest decay time which could be measured was of the 
order of 50 nsec. 

Relative intensities were measured using a P16 reference phosphor as 
a standard. The reference slide was changed sufficiently often so that 
deterioration by electron and ion bombardment remained negligible 
during the time necessary for the measurements. 

Measurements of the fluorescent intensity in the demountable station 
were only possible in relative units. To determine the absolute values 
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Fig. 1— Schematic of the circuit for measuring decay times of fast-decay 
cathode-ray tube phosphors. 


of the energy conversion efficiency under operating conditions a few 
small sealed-off cathode-ray tubes were made with screens from three 
different phosphor samples. The emission intensity from these tubes 
operating under standard excitation conditions was then compared 
with a calibrated light source®?° consisting of a P16 phosphor excited 
by the long-life radioactive isotope Sr°°. For technical reasons a light 
flux of 1.1 X 10-* watts from a phosphor screen with P16 spectral 
distribution was called one “‘R.” For green light as emitted by a P24 
phosphor the corresponding figure is 1.5 107 watts. 


2.3 Measurement of the Aging Characteristics 


The decrease of the luminescent efficiency of the screen with life in 
cathode-ray tubes is dependent on the total (time-integrated) deposited 
electron charge per unit area of the screen and the bombardment by 
negative ions produced by electron ionization of the residual gas also 
contribute. 

To minimize ion effects, as well as contaminating deposits from pump 
fluid vapors and rubber gasket joints which are unavoidable in demounta- 
ble stations, it was necessary to study the aging characteristics under 
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clean vacuum conditions. All aging data presented here were therefore 
obtained with sealed-off tubes using aluminized screens. The processing 
of the testers consisted of baking for three hours at 380°C, RI out- 
gassing of the metal parts, and flashing the getters. The seal-off pres- 
sure was about 10-7 mm Hg. 

All tubes were aged with direct current at 10 kv. The aged spots 
ranged in area from 15-25 mm?. A power density of 0.6 watt/cm? was 
found small enough to avoid damage to the screen by overheating. 
At power densities in excess of 1.5 watt/em?, the phosphor screen and 
the aluminum backing volatilize. Each phosphor was tested in at least 
two tubes. The reproducibility of the data was within the limits of +10 
per cent. 

Changes in luminescent efficiency were determined with a photocell 
having $4 spectral response. The spectral distribution of the fluorescence 
has been determined for the aged and unaged phosphors to ensure that 
no changes occurred which would influence the results. 


2.4 Other Characteristics 


The temperature dependence of the luminescent efficiencies was 
determined in the demountable tube. The phosphors were deposited 
on a nickel strip which could be heated directly by passing a current 
through it. The same arrangement was used for determining the tem- 
perature dependence of the decay time. 

The dependence of the decay time and luminescent efficiency upon 
exciting current density was measured in sealed-off tubes as well as 
in the demountable cathode-ray tube. 

Particle size distributions were determined with a Sharpless Micro- 
merograph. 

Spectral distributions were measured with a Bausch and Lomb 
Monochromator “Model 250 mm” and a photomultiplicr tube with an 
S11 cathode. An absolute calibration of the ensemble was made. 


III. RESULTS 


There is only a limited number of phosphor types with decay times 
of one microsecond or less and which have at the same time properties 
suitable for their incorporation into a vacuum tube. They are listed 
in Table I, where organic materials and inorganic ones with vapor 
pressures too high to be used in cathode-ray tubes were not included. 

Some of the general characteristics of these materials are given in 
publications" and a handbook edited by the various manufacturers.!? 
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TABLE I — PHospuors witH Drecay TIMES OF 1 uSEc or LEss, SUITABLE 
FOR INCORPORATION IN VAcuUM TUBES 











Phosphor Type Activator Designation Main Use 
(CaO)2-MgO: (SiOz). Ce P16 TV Scanner 
CaO-AlsO3- (Si02)2 Ce P16 TV Scanner 

ZnO Zn P15, P24 Color TV Scanner 
Zns Ag, Ni Color TV 
BaSO, Pb X-Ray Screens 


The present work includes all samples which were available at the 
time of the survey. Such a detailed study is valuable, since the charac- 
teristics of phosphors presented under one designation may vary con- 
siderably from batch to batch and from manufacturer to manufacturer. 

Tor comparison and control of the measuring equipment some meas- 
urements on CaWO, (P5), a cathode-ray tube and X-ray screen phos- 
phor having relatively well known properties, were included also. 

The phosphor samples were obtained from the following manufac- 
turers: RCA, Sylvania, GE, Du Pont, US Radium Corporation, Levy 
and West, Derby, EMI, Leuchtstoffwerke Heidelberg, Auer, Osram, 
and Philips. Data on a series of ZnO phosphors prepared by the author™ 
are also presented. 


3.1 Light Output vs Screen Weight 


The light output of a screen, observed from the face-plate side of the 
phosphor, goes through a maximum and then falls to zero as the screen 
thickness or weight is increased from zero to infinity (Fig. 2). The shape 
of the light output vs screen weight curve and the position of the maxi- 
mum depend primarily on the scattering and reabsorption of the fluo- 
rescent light in the screen layer. It is also a function of the penetration 
depth of the bombarding electrons. The problem has been treated 
theoretically by Hamaker." 

For phosphors with a single emission band, the maximum light output 
is obtained when the faceplate is just fully covered with phosphor 
crystals and no gaps exist between the grains so that no part of the elec- 
tron beam can hit the glass directly.1> Complete coverage is obtained 
with about one to two layers of particles, depending on the particle 
size distribution, which in turn determines the filling factor of the layer. 
If the phosphor powder contains only small particles of uniform size, 
the coverage and consequently the light output will vary rapidly as a 
function of the screen weight, first increasing steeply, passing through a 
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EMISSION INTENSITY —» 





SCREEN WEIGHT IN MILLIGRAMS PER CM@ 


Fig. 2 — Luminescent intensity versus screen weight for three representative 
phosphors: (a) P16 sample #1, (b) P16 sample #2, (c) P15. Samples used for 
obtaining curves (a) and (b) have different particle size as shown in Fig. 3. 


maximum, and decreasing rapidly. For a phosphor with a broader 
particle size distribution, the dependence of the light output upon the 
screen weight will become less critical. The relations are illustrated by 
Figs. 2 and 3. Phosphor sample *1 of Fig. 3 was used to obtain data of 
Vig. 2(a) and sample *2 for Fig. 2(b). 

The phosphor shown in Fig. 2(c), a ZnO(Zn)-P15 type, is representa- 
tive of a phosphor with two emission bands (blue and green for the 
ease of P15). It is believed that the unusual form of the curve is caused 
by the strong reabsorption of the edge emission band, which enhances 
the green emission band. The phosphor therefore has a normal intensity 
increase with increasing screen weight up to the point of complete 
coverage. With further increase in screen weight, the light output 
increases further until the screen thickness is great enough to absorb the 
blue light completely. For still greater screen thicknesses the light 
output then decreases normally. 

Table II gives a few representative values of optimum screen weight 
for maximum light output at 10 kv. 

There is clearly a wide variation of optimum screen weights for these 
phosphors, all of which are designated P16, which is believed to be due 
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Fig. 3 — Particle size distribution of three representative P16 phosphors. 


mainly to differences in particle size distribution. If reproducibility 
of screen performance is desired (during the settling of the final screens), 
care must be taken to control the particle size distribution because 
resolution and aging characteristics, as well as optimum screen weight, 
depend on it. If, for example, the average particle size is reduced by 


Taste II] — Maximum Licut Ovurrut FoR SCREENS POURED 
IMMEDIATELY AFTER AGITATION OF THE SUSPENSION. 
AutL P16 SAMPLES 





Sample # Composition ae hteaaed in Woe eo, 
1 CaMgSilicate 100 (sharp max.) 1.5 
2 CaMgSilicate 100 (flat max.) 2.8-4.2 
3 CaMgSilicate 74 1.6 
4 CaMgSilicate 90 2.8 
5 CaAlSilicate 31 4.2 
6 CaMgSilicate 103 1.4 
7 CaAlSilicate (?) 49 4.2 
8 CaMgSilicate 43 1.5 
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RELATIVE EMISSION INTENSITY, I 
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Fig. 4 — Effect of mechanical treatment of a P24 phosphor on its emission 
intensity and optimum screen weight. 


permitting the larger particles to settle out of the suspension for a few 
seconds or minutes (presettling) before pouring the screen, the optimum 
screen weight will be less. 

Heavy grinding of the phosphor should be avoided, whenever possible, 
as it may lead to mechanical deformations in the crystal structure and 
a consequent decrease in the emission intensity.“ Gentle grinding may 
sometimes improve the screen characteristics as shown in Iig. 4. Here 
a P24 type phosphor with many large particles was gently crushed in a 
mortar. It is evident that the maximum intensity has not been changed 
by this treatment, but that the screen weight necessary to reach the 
maximum intensity was greatly reduced. 

It is concluded from these data that it is necessary to determine the 
optimum screen weight for every settling process and every batch of 
phosphor, to ensure that optimum conditions are always maintained. 

The particle size distribution and screen thickness could also affect 
the amount of light scatter within the phosphor screen and hence 
influence the screen resolution. Very few systematic investigations of 
the screen resolution could be found in the literature, but to a first 
approximation it has been established that the resolution is degraded 
in proportion to the screen thickness and is independent of the particle 
size (as long as the latter is small compared with the screen thickness). 
The reason for this behavior is that, for increasing screen thickness, 
the light is scattered more frequently. Table III gives some values for 
the resolution of cathode-ray and X-ray screens found in the literature. 
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TABLE IJI — RESOLUTION OF CaTHODE-RAY AND X-RAY SCREENS 

















Screen Type Screen Thickness in» | Particle Size in p Resolution in » Ref. No. 
CR 10 10 1-6 (35) 
CR 35 (estimated) <10 52 (36) 
CR 35 (estimated) <30 49 (36) 
CR = 50 unknown 70 (37) 

X-Ray 300 unknown 500 (38) 








3.2 Hmission Intensity and Decay Time for Unaged Phosphors 


In Fig. 5, the luminescent intensity of the screen element is shown 
schematically as a function of time. The rise and decay follow an expo- 
nential law if the electronic process is confined to a luminescent atom, 
ion, or molecule, depending only on the transition probability of the 
luminescent transition in the center. If the luminescence is the result 
of the recombination of electrons in the conduction band with ionized 
centers, the rise and decay curves follow a second-order (bimolecular) 
law, with perturbations if traps are present. 

The rise time will in general depend on the exciting current density. 
This results from the fact that an equilibrium value between the exciting 
and emitting processes must be reached after the beginning of the 
excitation and that the transition probabilities for the excitation are in 
general independent of those of the emission. 

Under the present experimental conditions and for the phosphors 
studied, it was found that the rise time was always short relative to the 
5-microsecond exciting pulse length, and no systematic measurements 
of this parameter were therefore made. 


EXCITATION EXCITATION 
ENDS 


LUMINESCENT 
INTENSITY 





| RISE Le DECAY _ Le 
TIME TIME 


Fig. 5 — Luminescent intensity of a screen element as a function of the exci- 
tation time. 
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TaBLe TV — Emission Intensity, J, AND Decay TIMs, 7 (1/e POINT), 
FOR UNAGED SAMPLES OF I'ast-DEcay PHospHors. 10 Kv, 5 MA/cM?, 
PuLsE LENGTH 5 wSEC, REPETITION FREQUENCY 1000 TIMES PER SEC 















































Phosphor Type Sample * ge re Decay Time in psec 
P16 1 (Ref.) 100 0.20 
(CaO)2-MgO- (SiO2)e 2 100 0.25 
Ce-activated 3 74 0.25 
4 90 0.25 
6 103 0.20 
8 43 0.30 
P16 5 31 0.10 
CaO- Al.O3- (SiOez)e 7 49 0.15 
Ce-activated 
P15 1 110 1.1 
ZnO 2 100 0.9 
3 81 0.75 
P24 1 101 1.1 
ZnO (Zn) 2 105 0.9 
3 123 1.0 
4 85 0.45 
5 98 0.95 
6 30 0.25 
7 60 0.40 
8 88 0.75 
9 40 0.45 
10 85 0.70 
11 76 0.75 
12 55 0.45 
ZnO (Zn) 1 80 0.60 
Laboratory Preparations 2 71 0.55 
3 65 0.35 
4 36 0.20 
5 31 0.15 
ZnS (Ag,N1) 50 0.15 
BaPbSO, 65 0.50 
CaWO, Control 50 7.5 











Table IV lists the data obtained on the emission intensity, 7, and 
decay time, 7, of the unaged samples from various sources. The intensity 
values are those for optimum screen weight. Some of these data are 
included in Fig. 6, where the decay time is plotted as a function of the 
absolute intensity. For the phosphors of the P16 and P24 groups, there 
is clearly a strong relation between these two parameters, which can 
be expressed to a first approximation in the form J = kr(k = constant). 
A similar behavior can be expected for ZnS(Ag,Ni), but a complete 
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Fig. 6 — Intensity versus decay time for unaged fast-decay phosphors. One R 
corresponds to a light flux of 1.1 X 10-3 watts from a phosphor screen with P16 
spectral distribution. 


investigation was not possible becausc of the limited availability of 
samples. Further work on ZnO(Zn) phosphors prepared in the labora- 
tory showed that below an upper limit, which under the given excitation 
conditions is of the order of 1 microsecond, any desired 7 value can be 
obtained by using appropriate preparation conditions." 

The proportionality between J and 7 for ZnO(Zn) is believed to be 
due to the fact that crystals with shorter decay time have a greater 
number of lattice defects introduced by firing at higher temperatures." 
Radiationless transitions are induced by these defects and compete 
with the luminescent transitions in the excess Zn ions. Quenching effects 
can also be obtained by the introduction of special impurities, the so- 
called “killers,” such as Ni, Co, or Fe in ZnS.'® Quantities as low as 
2 X 10-5 per cent of Ni in ZnS(Cu) are sufficient to have an appreciable 
action on deep-lying traps, thereby quenching the long-lasting after- 
glow without materially affecting the luminescence efficiency during 
excitation.” A reduction of the emission intensity at a higher Ni con- 
centration appears with complete extinction at 10-2 per cent Ni. The 
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ZnS(Ag,Ni) sample studied here is representative of the application of 
this quenching principle to obtain an appropriate combination of decay 
time and intensity. 


3.3 Aging of Efficiency and Decay Time 


Data on the changes of efficiency and decay time (aging) under 
electron bombardment are especially important, as published data on 
quantitative measurements are very limited. 

The aging of the luminescence efficiency is shown in Fig. 7 for several 
examples of P16 phosphors, in Fig. 8 for one sample each of ZnO(Zn), 
ZnS(Ag,Ni), and BaPbSO,. Differences appear in the aging behavior 
of the P16 samples. This is due to differences in the manufacturing 
processes of samples of similar chemical composition as well as to the 
fact that some of the samples (*1, *2, *4, and 6) are of the (CaQ),- 
Mg0O-(SiO2)2 type, and others (#5 and 7), of lower initial intensity, 
are of the CaO- Al.O3- (SiO2)2 variety. All P16 samples, especially those 
of the Mg-silicate type, are characterized by a sharp percentage de- 
crease of the intensity during the first 0.1 coulomb/cm?. This is not too 
apparent when plotted on a logarithmic scale, but the rapid initial 
intensity decrease makes screen preaging useful where a high degree of 
screen uniformity is needed. 

For two P16 samples (#1 and *5), the aging of the efficiency was 
also determined at a lower beam acceleration voltage of 7.5 kv, and no 
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Fig. 7 — Aging of emission intensity of various P16 samples. 
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Fig. 8 — Aging of emission intensity of ZnO(Zn), BaPbSO,, and ZnS(Ag,Ni) 
phosphors. 


difference of the aging curve as compared to 10 kv was found. Also, 
even though the range from 7.5 to 10.0 kv is small, it confirms previous 
evidence of the independence of the aging with accelerating voltage.!8 

Some variations were found in the aging behavior of ZnO(Zn) = 
P24 samples. Phosphors prepared at firing temperatures over 1000°C 
aged slightly faster than those prepared below this temperature.” 
The differences are not great enough to be of importance for the present 
evaluations. Otherwise, the form of the normalized aging curve remains 
the same, regardless of the initial intensity of the P24 sample. 

Rapid aging occurs at the beginning of the bombardment for Ba: 
PbSO,, similar to P16. The ZnS sample, on the other hand, ages more 
slowly, approximating the P24 type phosphors. 

In the case of Ce-activated P16, it was proposed!® that the electron 
bombardment might act as a reducing agent, transforming the lumines- 
cent Ce*+ form of the activator to the nonluminescent Ce*+. As a con- 
sequence, it was proposed to compensate for this reducing action by 
introduction of an oxidizer into the screen of the cathode-ray tube.?° 
Experiments by R. C. Miller*! and the author on several P16 samples 
were unsuccessful in achieving the desired result. 

No difference in the aging curves for a P16 and a P24 sample was 
found between pulsed or direct current excitation. The result is impor- 
tant as it shows that accelerated aging tests can be performed with 
validity, at least within reasonable limits. 

The aging data given in this survey are valid only to the first approxi- 
mation, as phosphors with the same chemical composition and the 
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same activators may show variations in behavior depending both upon 
ditferences in physical characteristics of the material and operating 
conditions. Particle size, form of crystallization, operating temperature, 
vacuum, and mode of excitation are all important parameters that must 
be known and maintained under control. 

The decay time of all phosphors studied decreased with increasing 
electron bombardment. The greatest change occurred with P16 type 
samples (CaAl-silicates) for which a deposit of 0.04 coulombs/cm? re- 
duced the decay time from about 200 nsec to 50 nsec or less. After this 
initial decrease, further reduction of 7 was small. The decay time of the 
P24 type phosphors decreased by about 20 per cent after 40 coulombs/ 
cm2, independent of the initial value. Similar behavior was found for 
aging of ZnS(Ag,Ni) and BaPbSO,. 

The principal aspects of aging behavior under electron bombardment 
can be summarized as follows: The aging data can be represented, to a 
first approximation, by the formula J = Jo/(1 + CN), where J» = ini- 
tial intensity, J = aged intensity, C = burn parameter (a constant for a 
given phosphor), and N = number of electrons deposited per cm?.22:%3 
Phosphors with larger particles age more slowly than those with smaller 
particles.2* The form of crystallization has some influence on the aging, 
with phosphors quenched from higher firing temperatures aging faster 
in general.” Aging rates increase at higher sample temperatures, but 
the acceleration voltage has practically no influence.’ 

The mechanism of the aging process depends on the specific material 
and the type of exciting radiation. In general, there are two possibilities: 
(a) the bombardment destroys the luminescent center, e.g., changes the 
valence of the activator ion, or (b) it introduces defects which act as re- 
combination centers for radiationless transitions.2”:3 The first process 
probably applies to all materials where the luminescent process takes 
place in one ion or group of ions, such as organic molecules, CaWO, , 
etc. Phosphors with second-order kinetics, such as ZnS, ZnO, etc., are 
probably in the second group. If the phosphor is exposed to exciting 
radiation in an atmosphere which reacts with the material, further 
changes may take place. A well known phenomenon illustrating this 
effect is the photolysis of ZnS.?9 


3.4 Other Important Parameters 


3.4.1 Form of Decay 


Tig. 9 shows the decay curves (form of the decay) of the various phos- 
phors studied. None of them showed an exponential decay, except 
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CaW0O,, which previously has been measured*° and is shown here only 
for comparison. As would be expected for phosphors with recombination 
type of luminescence, at least a portion of the decay is due to the pres- 
ence of traps. This explains the long tail of the afterglow at low levels 
of intensity.*! In ZnO(Zn) the presence of traps was shown by glow 
curves” and conductivity measurements.® 

During the excitation and at the beginning of the decay, the concen- 
tration of excited electrons and ionized centers is large. Thus, the rate 
of the second-order processes, e.g. transition of electrons from the con- 
duction band to the excited or upper state of the center, can become 
larger than the transition within the center which is responsible for the 
luminescence. Consequently an exponential decay could be observed ini- 
tially with a decay constant characteristic for the transition within the 
center (activator). This effect is clearly present at the beginning of the 
decay of P16 and BaPbSO, phosphors (I'ig. 9). 

The increase in the recombination rate with increasing excitation 
density results in a decrease of the decay time. Fig. 10 shows this effect 
for P16 and P24 phosphors. Similar behavior was found for BaPbSO, 
and ZnS(Ag,Ni). 






100 

AND re eae aca (I aR PR He | 
80 
56 [eh tcawo, | | 












NC OELES 
RCEECEEEL EE 
ANDSEECEC EL 


20 








PER CENT OF INITIAL INTENSITY 


TIME IN MICROSECONDS AFTER END OF EXCITATION 


Fig. 9 — Form of decay of various fast-decay phosphors. 


FAST-DECAY PHOSPHORS 197 


DECAY TIME IN MICROSECONDS 





fo ae aaa (0 i ee a | | 
CURRENT DENSITY IN AMPERES PER CM2 


Fig. 10 — Decay time of a representative P16 and P24 sample as a function of 
the excitation level (as determined by the exciting current density). 


3.4.2 Temperature Dependence of Intensity and Decay Time 


An increase in temperature enhances, in general, the probability of 
radiationless transitions and should therefore result in a reduction of J 
and 7. This is shown in Figs. 11 and 12 for the fast-decay phosphor types 
studied here. ZnO(Zn), ZnS(Ag,Ni), and BaPbSO, exhibit a decrease 
of both parameters with increasing temperature. P16 is an exception, 
as its intensity increases slightly above room temperature and decreases 
for temperatures exceeding 150°C and 400°C for the aluminum and 
magnesium silicates respectively. This is similar to the behavior of Mn- 
activated silicate phosphors* which are frequently preferred for projec- 
tion kinescopes despite a lower intensity at room temperature as com- 
pared with ZnS-type phosphors. 


3.4.3 Spectra 


The spectra of the various fast-decay phosphors are shown in Fig. 18. 
No variation of the spectral distribution was observed during the aging 
for any of the fast-decay samples. 


IV. CONCLUDING REMARKS 


The data obtained in this study show that phosphors obtained under 
one designation and of the same chemical composition show wide varia- 
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Fig. 11 — Temperature dependence of emission intensity for various fast- 
decay phosphors. 
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Fig. 12 — Temperature dependence of the decay time of fast-decay phosphors. 
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Fig. 13 — Distribution of the spectral emission of various fast-decay phosphors. 


tions of their properties. This, of course, is to be expected as the lumi- 
nescent properties depend on the incorporation of small amounts of ac- 
tivators as well as on crystal defects which might be introduced during 
the preparation. 

Several samples of a given phosphor must therefore be evaluated to 
obtain the range of variations of its intensity, decay time, and optimum 
screen weight. In this manner the optimum phosphor can be selected 
for a specific application. Such a survey will be especially necessary if 
the screen performance is near the limit of the system requirements. 


Vv. ACKNOWLEDGMENTS 


The encouragement of R. W. Sears throughout the study has been 
most valuable. R. S. Ehle has performed many measurements and his 
help is greatly appreciated. 


200 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


REFERENCES 

1. DeBoer, F., and Emmens, H., Sedimentation of Fluorescent Screens in 
Cathode-Ray Tubes, Philips Techn. Rev., 16, 1955, p. 237. 

2. Vail, J. G., Soluble Silicates, Vol. 2, Their Properties and Uses, New York, 
Reinhold, 1952, p. 235. 

3. Sadowski, M., The Preparation of Luminescent Screens, J. Electrochem. Soc., 
95, 1949, p. 112. 

4. Hangos, I., Die Adhision von microkristallinischen Phosphoren an Glasfli- 
chen, Exp. Techn. d. Phys., 7, 1959, p. 71; 6, 1958, p. 228. 

5. Bril, A., and Klasens, H. A., Intrinsic Efficiencies of Phosphors under Cath- 
ode-Ray Excitation, Philips Res. Rep., 7, 1952, p. 401. 

6. Epstein, D. W., and Pensak, L., Improved Cathode-Ray Tubes with Metal- 
Backed Luminescent Screen, RCA Rev., 7, 1946, p. 5. 

7. Dowling, P. H., and Sewell, J. R.; The Effect of Organic Vapor on the Second- 
ary mission of Phosphors, J. Appl. Phys., 25, 1954, p. 228. 

8. Rottgardt, K. H. J., and Berthold, W., Sticking Potential of CR Screens, 
J. Appl. Phys., 26, 1955, p. 1180. 

9. Hoover, C. W., Jr., and Raag, H., Method for the Measurement of Radiant 
Flux, in Solid State Physics, New York, Academic Press, 1960, p. 688. 

10. Fuls, i. N., and Wittwer, N., private communication. 

11. Bril, A., and Klasens, M. A., New Phosphor for Flying Spot Cathode-Ray 
Tubes, Philips Res. Rep., 7, 1952, p. 421. 

Smith, A. L., Some New Complex Silicate Phosphors Containing Ce, Mg and 
Be, Trans. Am. Electrochem. Soc., 96, 1949, p. 287. 

Monnot, R., TV Transmission of Pictures of Variable Transparencies by 
Systems Which Make No Use of Storage, L’Onde Electrique, 30, 1950, p. 
377. 

12. Optical Characteristics of Cathode-Ray Tube Screens, formulated by JEDEC 
Electron Tube Council N. Y., June 1960. 

13. Pfahnl, A., Preparation and Properties of ZnO Phosphors, J. Electrochem. 
Soc., 109, 1962, p. 502. 

14. Hamaker, H. C., Radiation and Heat Conduction in Light Scattering Mate- 
rial, Philips Res. Rep., 2, 1947, pp. 5, 103, 112, 420. 

15. Leverenz, Luminescence of Solids, New York, John Wiley, 1950, p. 381. 

16. Byler, W. H., Studies on Phosphorescent ZnS, J. Am. Chem. Soc., 60, 1938, 
p. 632. 

17. Levy, L. A., and West, D. W., New and Rapid Intensifying Screens, Brit. J. 
Radiology, 6, 1933, p. 85. 

18. Rottgardt, K. H. J., Destruction of the Luminescence of Cathode-Ray Tube 
Screens by Electrons, Z. angew. Phys., 6, 1954, p. 160. 

19. Koontz, D., private communication. 

20. Cook, J. C., and Calbick, C. J., Luminescent Coating, US Patent 2,252,552. 
August 12, 1941. Sadowsky, M., and Groner, M. G.; Burn Resistant Lumi- 
nescent Products and Method of Making Same, US Patent 2,689,830. Sep- 
tember 1, 1954. 

21. Miller, R. C., private communication. 

22. Hanle, W., and Rau, K. H., Light Efficiency and Burn of Phosphors under 
Electron and Ion Excitation, Z. Phys., 133, 1952, p. 297. 

23. Rottgardt, K. H. J., Destruction of the Luminescence as a Consequence of 
Electron Bombardment, Naturwiss., 40, 1953, p. 315. 

24. See Ref. (15), p. 449, also Berthold, W., Reversible Phosphor Damage by 
Positive Ions, Naturwiss., 42, 1955, p. 436. 

25. Martin, W., Change of Luminescent Properties of Phosphors by H:t Ions, 
Z. Phys., 147, 1957, p. 583. 

26. Short, M. A., and Steward, E. G., The Effect of Grinding on the Structure 
and Luminescence of ZnS and ZnCd§, Z. Phys. Chem., 18, 1957, p. 298. 

27. Bandow, F., Lumineszenz, Stuttgart, Wiss. Verlagsges, 1950, p. 166. 

28. Broser, I., and Reichard, W., On the Relation Between Light Output and 


Particle Size in Polycrystalline Crystal-phosphors, Z. Naturf., 6a, 1951, 
p. 466. 


29. 
30. 


ol. 
32. 
33. 
34. 
35. 
36. 


37. 
38. 
39. 


FAST-DECAY PHOSPHORS 201 


Schleede, A., On the Darkening of ZnS by Light, Z. Phys. Chem., 106, 1923, 
p. 386. 

Pfahnl, A., Study of an X-Ray Tube with Control Grid, Application to the 
Study of Phosphorescence, Rev. Gen. Electricité, 66, 1957, p. 159. 

Coltmann, J. W., Ebbinghausen, E. G., and Altar, W., Physical Properties 
of Calcium Tungstate X-Ray Screens. J. Appl. Phys., 18, 1947, p. 580. 

See for example Kréger, F. A., Inorganic Crystal Phosphors, in Ergeb. d. 
exakten Naturwiss., Berlin, Springer, 1956, p. 80. 

Garlik, G. F. J., and Wilkins, M. H. F., Short Period Phosphorescence and 
Electron Traps, Proc. Royal Soc., (London), A184, 1945, p. 408. 

Weiss, H., The Photoelectric Conductivity of ZnO, Z. Phys., 182, 1952, p. 
335. 

Kroger, F. A., Some Aspects of the Luminescence of Solids, Amsterdam, E]se- 
vier Pub. Co., 1948. 

Cooper, H. G., Design of a High-Resolution Electrostatic Cathode-Ray Tube 
for the Flying Spot Store, B.S.T.J., 40, 1961, p. 723. 

Deubner, B., and Hieber, F., Vergleichende Messungen an Aufsicht- und 
Durchsicht Leuchtschirmen bei Elektronenanregung, Z. angew. Phys., 6, 
1954, p. 112. 

Borries, V. B., Die energetischen Daten und Grenzen der Ubermikroskopie, 
Optik, 3, 1948, p. 321. 

Frost, D., Uber die Bildunschirfe bei Verstarkerfolien und Leuchtschirmen, 
Roéntgen-Blitter, 10, 1957, p. 1. 

Hoover, C. W., Jr., Staehler, R. E., and Ketchledge, R. W., Fundamental Con- 
cepts in the Design of the Flying Spot Store, B. 8S. T. J., 37, September, 
1958, p. 1161. 


Contributors to this Issue 


Cuaruss D. Bracu, B.S. in 8.E., 1959, University of Colorado; M.S. 
in E.E., 1961, New York University; Bell Telephone Laboratories 1959- 
1962. His work at the Laboratories involved theoretical analyses of 
interchannel modulation distortion in tropospheric scatter radio sys- 
tems. He also participated in an experimental study of intermodulation 
performance of such systems. Since February, 1962, he has been on 
leave of absence from the Laboratories while teaching and pursuing 
doctoral studies at the University of Colorado. Member Tau Beta Pi, 
Eta Kappa Nu, Sigma Tau. 


Rover R. Buarr, B.S. in E.E., 1930, Rose Polytechnic Inst.; M.A. in 
Physics, 1939, Columbia University; Bell Telephone Laboratories, 
1930—. Mr. Blair’s initial work at the Laboratories in transmission re- 
search led to the design of circuits for broadband amplifiers, regulators 
and test equipment for coaxial systems. During World War II he was 
concerned with circuits for proximity fuses, radar indicators and other 
security-classified devices. There followed a period in apparatus develop- 
ment on germanium and silicon varistors and on processing and test 
equipment used in that field. He later developed high-speed transistor 
switching circuitry for digital computer applications for the military. 
Since 1958, he has been investigating and studying the effects of nuclear 
radiation on various semiconductor devices and on their performance in 
a wide variety of applications. He is a senior member of the Institute of 
Radio Engineers and a member of Tau Beta Pi. 


W. F. Boprmanwn joined Bell Telephone Laboratories in 1941 and 
during 1942-1945 served in the U.S. Air Corps. After returning to the 
Laboratories in 1945, he worked on various circuits and components for 
the TD-2, TJ and TL radio relay systems, and on transistor circuits for 
wideband amplifiers. Later he was concerned with the FM feedback 
receivers for Project Echo, and more recently has been working on the 
frequency compression demodulators for Project Telstar. 


W. L. Brown, B.S., 1945, Duke University; A.M., 1947, Ph.D., 1951, 
Harvard University; Bell Telephone Laboratories 1950—. Mr. Brown 


203 


204 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


has been engaged in research on the physical properties of semiconductor 
surfaces and the nature of defects produced in semiconductors by high 
energy radiation. He has recently been responsible for radiation damage 
experiments aboard the Telstar satellite. Fellow of American Physical 
Society, Sigma Xi and Phi Beta Kappa. 


James C. Coyne, A.B., 1953, B.S., 1954, and M.S., 1958, Columbia 
University; Bell Telephone Laboratories, 1958—. Mr. Coyne partici- 
pated in the planning for the Columbus Branch Laboratory and was en- 
gaged in development work on welded connections. At present he is en- 
gaged in reliability and application studies of solderless wrapped and 
percussive welded connections. Tau Beta Pi. 


IiucEne I. Gorpvon, B.S8., 1952, City College of New York; Ph.D., 
1957, Massachusetts Institute of Technology; Bell Telephone Labora- 
tories, 1957—. A member of the Electron Device Department, he is en- 
gaged in research in optical masers. Member American Physical Society, 
Phi Beta Kappa, Sigma Xi. 


H. JAMES GraFF, B.S.M.I., 1955, University of Notre Dame; Bell 
Telephone Laboratories 1956—. Mr. Graff has been concerned primarily 
with the development of a mechanical splicing method for polyethylene 
insulated conductor cable. At present he is in charge of a group responsi- 
ble for the development of load coils, coil cases, and other outside plant 
items. 


Jessip MacWiuuiaMs, B.A., 1939, M.A., 1941, Cambridge (Eng- 
land); Ph.D., 1962, Harvard; Bell Telephone Laboratories 1956—. Mrs. 
MacWilliams has been concerned with writing computer programs for 
the analysis and synthesis of transmission networks. She is now engaged 
in data systems studies, particularly the study of algorithms for decoding 
systematic error-correcting codes. Member Mathematical Association of 
America. 


Joun M. Peacock, B.S.E., 1947, Princeton University; Western Elec- 
tric Co., 1947-50; Bell Telephone Laboratories 1953—. At Bell Labora- 
tories Mr. Peacock first worked on the design of antennas, radio towers 
and other structures. Later he was engaged in fundamental development 


CONTRIBUTORS TO THIS ISSUE 205 


work on expanded dielectrics for cable conductor insulation. At present 
he is supervisor in charge of outside plant wire joining and cable splicing 
activities. Member, American Wire Association. 


D. Stewart Peck, B.8.E., 1939, M.S., 1940, University of Michigan; 
Bell Telephone Laboratories 1947—. Mr. Peck has been concerned in the 
design for production of gas-filled electron tubes such as rectifiers, thyra- 
trons and cold-cathode tubes. More recently, he has been in charge of 
work on reliability studies, applications engineering, and specifications 
for electron devices. On Project Telstar Mr. Peck’s department has been 
generally responsible for the specification, aging and selection programs 
for semiconductor components used in the satellite. Member A.I.E.E., 
American Standards Association, Tau Beta Pi, Eta Kappa Nu, Sigma 
Xi and Phi Kappa Phi. 


ARNOLD PFAHNL, Dr. Phil., University of Graz, Austria, 1948, D.Sc., 
University of Paris, France, 1956; G. Massiot & Cie. and Massiot-Fluor, 
Courbevoie, I'rance, 1950-57, Bell Telephone Laboratories, 1957—. At 
the Laboratories, Mr. Pfahnl has worked on the study of the properties 
of fast-decay cathode-ray tube phosphors, the development of ZnO phos- 
phors with application to cathode-ray tubes used in electronic switching, 
and the study of injection electroluminescence of gallium phosphide. 


J. DANE RiapEN, B.Sc., 1955, and Ph.D., 1958, Reading University, 
Reading, England; National Research Council of Canada, 1958-60; Bell 
Telephone Laboratories 1960—. Mr. Rigden has specialized in the de- 
velopment of gaseous optical masers and is the co-developer of the visible 
helium-neon maser. 


C. L. Rururorr, B.S., in E.E., 1950, and M.A., 1952, University of 
Nebraska; Long Lines, 1946-52; Bell Telephone Laboratories, 1952—. 
Since transferring from Long Lines in Lincoln, Nebraska, where he was 
a central office maintenance man, Mr. Ruthroff has concentrated on 
radio research with emphasis on FM problems and transistor circuits for 
microwave components and systems. Member, I.R.E. 


Frrepour M. Smits, Dipl. Phys., 1950, Dr. rer. nat., 1950, University 
of Freiburg, Germany; research associate, Physikalisches Institut, Uni- 


206 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1963 


versity of Freiburg, 1950-54; Bell Telephone Laboratories, 1954-62. Mr. 
Smits went to the Sandia Corporation in May 1962. His work at Bell 
Telephone Laboratories included studies of solid-state diffusion in ger- 
manium and silicon, device feasibility, and process studies, as well as 
the development of UHI semiconductor devices. He supervised a group 
that conducted radiation damage studies on components, particularly 
solar cells, used in the Telstar experimental satellite. Member of the 
American Physical Society and the German Physical Society. 


JosEPH M. Trecxer, B.S. in E.E., 1949, Iowa State College; Bell 
Telephone Laboratories, 1949—. His early assignments involved studies 
of information requirements and planning of communication systems 
for military applications. More recently, he was engaged in studies of 
transmission characteristics of tropospheric scatter radio systems. He is 
now in charge of a group studying information flow and ground com- 
munication network arrangements for the Apollo space mission. Member 
Tau Beta Pi, Eta Kappa Nu. 


J. JoHN ZALMANS, B.S.M.E., 1953, University of Arizona; Bell Tele- 
phone Laboratories, 1953—. Mr. Zalmans is engaged in the development 
of connecting devices. American Society of Mechanical Engineers, Tau 
Beta Pi, Pi Mu Epsilon. 


